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The  work  presented  in  this  report  was  performed  by  Bolt  Beranek  and 
Newman  lne.,  Cambridge,  Massachusetts,  under  USAF  Contract  No.  F33615-70-  fW 
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Problems  in  Air  Force  Flight  Vehicles",  Task  147102,  "Prediction  and 
Control  of  Noise  Associated  with  USAF  Flight  Vehicles",  and  administered  by 
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Smith  served  as  Project  Engineer;  his  assistance  and  cooperation  are  grate¬ 
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Inc.  Report  No.  2104. 
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appendix  on  "Atmospheric  Propagation  Effects"  and  assembled  that  part  of  Lne 
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Dr.  Ulrich  J.  Kurze  authored  the  appendix  on  "Sound  Attenuation  in  Circular 
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combustion ,  and  turboshaft  engines.  Mr.  Davey  L.  Smith  and  Gapt,  Rodney  P. 
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Determining  Uncorrected  Detection  Range".  Dr.  Eric  h .  iJngar  was  responsible 
for  tiie  coordination  and  integration  of  all  efforts  into  Luis  report,  as  well 
as  for  die  rest  of  the  technical  work  presented  Here. 
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ABSTRACT 


The  concepts  which  underlie  the  detection  of  aircraft  by 
the  human  ear  are  described.  A  scheme  Is  delineated  for  predic¬ 
ting  the  range  at  which  a  given  aircraft  will  first  be  heard  by 
an  "average"  listener;  this  scheme  is  also  applicable  to  compar¬ 
ing  the  aural  detectabilities  of  alternate  aircraft  configura¬ 
tions  and  to  identifying  those  components  of  the  noise  of  a  given 
configuration  which  bear  prime  responsibility  for  the  aircraft's 
detectability.  Means  are  presented  for  predicting  the  noise  due 
to  all  sources  likely  to  be  significant  for  light  aircraft,  the 
attenuation  of  acoustic  signals  propagating  from  an  aircraft  to 
a  listener  on  the  ground,  and  the  ability  of  a  listener  to  de¬ 
tect  an  acoustic  signal  In  the  presence  of  background  noise. 
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INTRODUCTION 


PURPOSE  AND  SCOPE 

It  is  desirable  for  all  military  aircraft  to  be  impossible 
to  detect  by  the  enemy  until  it  is  too  late  for  him  to  employ 
countermeasures  or  take  evasive  action.  Modern  high-speed  strike 
or  reconnaissance  aircraft  approach  this  ideal  by  flying  at  tree- 
top  level,  but  their  high  speeds  limit  their  utility  for  certain 
missions.  Consequently,  there  exists  a  continuing  interest  in 
slow-flying  aircraft  of  minimum  detectability. 

Much  is  known  about  designing  aircraft  and  countermeasures 
to  reduce  detectability  by  radar,  optical  and  infrared  means; 
also,  all  of  these  means  tend  to  be  useful  only  for  very  limited 
ranges  against  aircraft  flying  at  treetop  level.  In  contrast, 
the  acoustic  noise  emanating  from  aircraft  may  be  detected  -  and 
localized  -  at  very  considerable  ranges,  even  by  the  unaided 
human  ear.  Thus,  low-speed  aircraft  tend  to  be  extremely  vulner¬ 
able  to  aural  detection  by  personnel  on  the  ground. 

Although  there  exists  much  information  relevant  to  the  de¬ 
sign  of  "quiet"  aircraft,  this  information  is  widely  scattered 
throughout  the  literature  and  is  generally  inaccessible  to  the 
nonspecialist  in  acoustics.  It  is  the  purpose  of  the  present 
report  to  provide  a  collection  of  information  pertinent  to  the 
design  of  aircraft  from  the  standpoint  of  minimum  aural  detection, 
and  to  present  this  information  in  a  form  in  which  it  may  be  read¬ 
ily  used  without  reference  to  a  multitude  of  other  documents. 

A  "quiet"  aircraft  must  satisfy  aural  detectability  specifi¬ 
cations,  in  addition  to  the  usual  performance  specifications  (e.g., 
payload,  range,  speed).  Generally,  many  different  alternative 
configurations  will  satisfy  the  various  requirements  to  various 
dppr>pest  3nd  the  aircraft  designer  usually  is  faced  with  selecting 
the  best  compromises.  The  designer  of  an  aircraft  that  must  meet 
stringent  aural  detectability  requirements  will  generally  need  to 
begin  with  a  first-cut  configuration  designed  on  the  basis  of  both 
performance  and  acoustic  requirements;  he  must  then  evaluate  the 
aural  detectability  of  his  initial  design,  identify  what  acoustic 
changes  are  required,  modify  or  redesign  the  aircraft  accordingly 
-  and  repeat  the  evaluation-modification  cycle  until  he  reaches 
an  acceptable  result. 
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The  information  provided  ir.  the  present  report  may  be  used: 
(1)  to  compare  the  aural  detectabilities  of  alternate  components 
(in  order  to  enable  one  to  make  favorable  initial  choices),  (2) 
to  evaluate  and  compare  the  aural  detectabilities  of  aircraft 
designs,  (3)  to  estimate  the  ranges  at  which  aircraft  may  be  de¬ 
tected  by  ear,  and  (4)  to  identify  what  components  of  noise  con¬ 
tribute  most  to  the  aural  detectability  of  a  given  design  (and 
thus  what  components  of  noise  must  be  reduced  in  order  to  reduce 
detectability).  Although  some  hints  are  provided  concerning  tac¬ 
tics  (e.g.,  selection  of  favorable  terrain  and  meteorological  con¬ 
ditions)  for  reducing  the  probability  of  aural  detection  and  con¬ 
cerning  how  one  might  obtain  more  precise  detection-range  predic¬ 
tions  for  given  situations,  these  topics  are  considered  generally 
beyond  the  scope  of  this  report. 

ORGANIZATION  OF  REPORT 

The  first  of  the  following  sections  summarizes  the  basic 
concepts  involved  in  aural  detectability  analysis,  presents  com¬ 
monly  used  acoustical  quantities  and  notation,  and  indicates  the 
most  important  characteristics  of  these  quantities.  The  second 
section  delineates  a  scheme  for  determining  the  range  at  which 
an  aircraft  can  first  be  heard,  provides  guidelines  for  obtaining 
the  information  needed  to  apply  this  scheme,  and  illustrates  its 
use . 

In  the  fourteen  appendixes  of  this  report  are  provided  the 
tools  for  implementing  the  suggested  detection  range  estimation 
scheme.  These  appendixes  deal  with  the  various  noise  sources 
likely  to  be  significant  for  light  "quiet"  aircraft,  with  the 
ability  of  human  listeners  to  detect  an  acoustic  signal  in  the 
presence  of  background  noise,  and  with  the  attenuation  acoustic 
signals  experience  as  they  propagate  through  the  atmosphere. 

Each  appendix  attempts  to  provide  the  reader  with  some  .insight 
into  the  basic  processes;  some  appendixes  -  in  particular  those 
which  present  newly  developed  information  -  even  contain  complete 
derivations.  Each  appendix  also  includes,  usually  at  its  end, 
specific  directions  for  application  of  the  given  information. 


SOME  BASIC  CONCEPTS 


AURAL  DETECTABILITY 


Aural  detectability  involves  three  major  considerations; 

1.  Characteristics  of  acoustic  noise  generated  at  the  source, 

2.  Modifications  experienced  by  acoustic  signals  as  they  propa¬ 
gate  from  the  aircraft  to  listener  locations  on  the  ground, 

3.  Ability  of  human  listeners  to  detect  aircraft  noise  in  their 
background  noise  environment. 

A  noise  source  may  be  described  most  conveniently  in  terms 
of  its  acoustic  power  spectrum  (i.e.,  the  frequency  distribution 
of  the  acoustic  power  emitted  by  the  source),  its  directivity 
(i.e.,  the  variation  of  the  acoustic  power  with  direction  as  mea¬ 
sured  from  reference  axes  attached  to  the  source),  and  any  time  - 
variations  of  the  spectrum  and  directivity.  Because  the  noise 
sources  associated  with  aircraft  in  steady  flight  are  not  likely 
to  change  with  time,  time-variations  need  not  be  considered  here; 
rather  slow  variations,  such  as  may  result  from  maneuvers  or 
changes  in  engine  or  propeller  settings  may  be  analyzed  simply  in 
terms  of  sequences  of  quasi-steady  conditions. 

In  most  military  situations  in  which  one  is  concerned  with 
aural  detection,  enemy  personnel  may  be  located  in  any  direction 
from  the  aircraft.  One  therefore  would  prefer  to  design  such 
an  aircraft  so  that  there  exist  no  directions  in  which  the  radi¬ 
ated  sound  Is  more  intense  than  in  others  -  except  that  one  may 
wish  to  direct  much  of  the  noise  upward,  where  there  are  no  lis¬ 
teners,  or  possibly  backward,  where  the  noise  heard  at  a  fixed 
location  on  the  ground  after  the  aircraft  has  passed  is  likely  to 
be  less  intense  than  when  the  aircraft  is  more  nearly  overhead. 
However,  for  most  noise  sources  of  interest  here  one  can  do  little 
to  change  their  directivity  significantly.  Although  the  various 
noise  sources  discussed  later  do  have  some  inherent  directivity 
characteristics,  the  directivity  effects  are  generally  not  very 
pronounced.  For  this  reason,  and  because  consideration  of  direc¬ 
tivity  effects  complicates  the  analysis  and  greatly  increases 
the  computational  burden,  directivity  considerations  are  not  em¬ 
phasized  here,  (But  directivity  information  has  been  included 
in  the  appendixes,  as  far  as  it  is  readily  at  hand,  so  that  it 
may  be  available  for  more  detailed  analyses  than  those  generally 
advocated  here.) 


As  an  acoustic  signal  travels  away  from  a  source,  its  energy 
is  spread  out  over  a  larger  area,  so  that  less  energy  reaches  a 
sensor  (ear  or  microphone)  of  a  given  size.  In  addition,  the 
atmosphere  through  which  a  signal  propagates  absorbs  energy  from 
it,  again  reducing  the  signal  that  reaches  a  given  sensor.  Atmos¬ 
pheric  gradients  also  refract  (deflect)  signals,  and  n.ay  either 
reduce  or  Increase  the  amount  of  energy  reaching  a  sensor;  the 
same  is  true  of  terrain  and  vegetation,  which  may  reflect  and/or 
absorb  acoustic  energy.  An  acoustic  signal  generally  is  described 
by  its  spectrum  (frequency-distribution)  of  acoustic  pressure; 
only  spreading  attenuates  all  frequency-components  equally,  the 
other  attenuation  effects  generally  are  frequency-dependent. 

A  human  listener  can  detect  an  acoustic  signal  only  if  it 
is  strong  enough  to  be  above  the  threshold  of  sensitivity  of  his 
hearing  system  and  if  it  is  sufficiently  intense  so  that  he  can 
distinguish  it  from  the  background  noise  that  always  reaches  his 
ear  (e.g.,  due  to  wind,  rustling  leaves,  machinery,  human  activ¬ 
ity).  The  hearing  sensitivities  of  different  Individuals  differ 
widely  and  are  also  affected  by  exposure  to  noise.  Of  course, 
any  ear  covering  tends  to  reduce  this  sensitivity,  whereas  lis¬ 
tening  aids  (even  cupped  hands)  may  enhance  it.  An  individual's 
ability  to  detect  an  acoustic  signal  also  depends  on  many  psycho¬ 
logical  factors,  such  as  fear  and  fatigue,  whose  effects  cannot 
be  readily  predicted.  The  "masking"  effect  of  background  noise 
depends  on  the  spectral  (and  time-dependent)  character  of  both 
the  noise  and  the  signal.  Prediction  of  whether  a  given  signal 
will  be  detected  (or  better,  of  the  probability  that  a  given  sig¬ 
nal  will'  be  detected)  under  given  conditions  clearly  is  a  very 
complex  problem.  Fortunately,  this  problem  and  most  of  its  com¬ 
plexities  may  be  by-passed  here,  since  for  purposes  of  evaluating 
alternate  aircraft  components  and  designs  it  suffices  to  consider 
an  average  listener  under  representative  conditions. 

ACOUSTICAL  QUANTITIES 

Sound  Power  Level  and  Sound  Pressure  Level 

Because  the  values  of  sound  power  and  sound  pressure  encoun¬ 
tered  in  most  acoustics  problems  extend  over  several  orders  of 
magnitude,  it  has  become  customary  to  represent  these  quantities 
in  logarithmic  terms.  Thus,  the  sound  power  level  Lw  of  a  source 
which  radiates  acoustic  power  W  is  defined  as* 

Lw  =  10  log  (W/Wref)  (1) 


*A11  logarithms  in  this  report  refer  to  base  10. 


where  Wref  is  a  reference  power.  The  internationally  accented 
standard  value  Wref  =  in-1-  watt  *  7.38  *  i0“13  ft  lb/sec  is  used 
throughout  this  report.  The  power  level  is  a  dimensionless 
quantity ,  but  to  indicate  its  logarithmic  character  and  the  fact 
that  it  is  defined  with  a  factor  of  10,  it  is  given  the  "units" 
of  decibels  (dB).  Because  different  reference  values  have  been 
used,  it  is  also  customary  to  indicate  the  reference  value  one 
uses.  For  example,  the  power  level  corresponding  to  an  acoustic 
power  of  2.0  watt  would  be  written  as  4,  =  123  dB,  re  10"12  watts. 

The  sound  pressure  level  Lp  corresponding  to  a  mean-square 
acoustic  pressure  p2  is  defined  in  a  manner  similar  to  Eq.  (1)  as 

Lp  -  10  log  (p2/pj,sf)  -  20  log  (p/pref)  (2) 

where  p£ef  represents  a  reference  mean-square  pressure,  and  where- 
p  and  Pj-ef  represent  the  root-mean-square  values  that  correspond 
to  p2  and  Pref*  The  internationally  accepted  standard  reference 
value  pref  =  0.0002  microbar  =  0.0002  dyne/cm2  =  2  >  10“ snewton/m2 

*  2  x  I0~10atm  *2.86  x  I0"spsi  is  used  throughout  this  report. 
Again,  it  is  customary  to  append  the  "dB"  designation  to  the 
sound  pressure  level  and  to  indicate  the  reference  quantity;  for 
example,  the  pressure  level  corresponding  to  an  acoustic  pressure 
of  p  *  10  microbar  would  be  written  as  L  =  9^  dB,  re  0.0002 
microbar.  p 

In  the  foregoing  definitions  W  and  p  may  represent  either 
the  total  power  and  mean-square  pressure  measured  at  all  frequen¬ 
cies,  or  those  same  quantities  measured  in  specified  frequency 
bands.  If  these  quantities  are  measured  in  bands  and  the  result¬ 
ing  Lw  and  Lp  values  are  plotted  as  a  function  of  frequency  (usu¬ 
ally,  the  center  frequency  of  the  bands  used),  one  obtains  sound 
power  and  sound  pressure  spectra.  Spectra  are  most  often  given 
in  octave  bands  (which  span  a  factor  of  2  in  frequency)  and  1/3 
octave  bands  (which  span  a  factor  of  21/*  in  frequency). 

Dependence  of  Sound  Pressure  on  Sound  Power  Level 

The  mean-square  acoustic  pressure(measured  at  a  given  fre¬ 
quency  or  in  a  given  frequency  band)  at  a  given  location  with 
respect  to  a  sound  source  is  proportional  to  the  acoustic  power 
produced  by  that  source  (at  the  same  frequency  or  in  the  same 
frequency  band).  In  a  free  field,  where  no  reflections  occur, 
the  ac-iustic  pressure  decreases  with  increasing  distance  from 
the  source;  at  several  wavelengths'  distance  from  the  source,  the 
mean-square  pressure  varies  inversely  as  the  square  of  the  dis¬ 
tance  from  the  source  (because  in  essence  the  acoustic  power 
emanating  from  the  source  is  spread  over  a  larger  surface  area 
and  a  sensor  of  a  given  area  intercepts  less  of  the  power) . 


The  aforementioned  variation,  when  expressed  logarithmically 
and  with  the  appropriate  constants,  results  in  a  simple  expression 
for  the  sound  pressure  level  Lp  measured  at  a  distance  R  from  a 
nondirective  source  (i.e.,  a  source  that  radiates  uniformly  in 
all  directions)  emitting  sound  with  a  sound  power  level  L  : 

L  (dB,re  0.0002  microbar)  -  L  (dB,re  10”12watt)  -  20  log  R  -  0,5  (3) 

IP  w 

where  R  denotes  the  distance  in  feet. 

Dependence  of  Levels  on  Bandwidth 

It  is  often  necessary  to  express  levels  given  in  bands  of 
one  width  in  terns  of  levels  in  bands  of  other  widths.  For  sig¬ 
nals  whose  narrow-band  spectra  are  relatively  flat  within  the 
broader  (overlapping)  bands  Afj  and  &f2 ,  one  may  assume  the  power 
to  be  approximately  uniformly  distributed  over  frequency  in  the 
vicinity  of  the  bands  of  concern.  If  W,  represents  the  power  in 
the  band  Afj  and  W2  that  in  Af2 ,  then 

W,/Af,  =  V/2  /  Af2  .  (4) 

The  power  level  Lw, in  the  band  Afj  then  is  ^elated  to  the  level 
Lw  2  ^  s 

LW,  a  I*2  +  10  lo?  (A  f  j  /A  f2  )  (5) 

If  Af,  is  an  octave  band  and  Af2  a  1/3-octave  band  of  the  same 
center  frequency,  then  Af,/Af2  =  3.0  and  Lwl  -  LW2  »  5  dB. 

Multiple  Sources 

The  total  acoustic  power  radiated  by  several  sources  is 
equal  to  the  sum  of  the  individual  contributions,  This  is  true 
of  the  overall  power  (including  all  frequencies),  as  well  as  in 
any  frequency  band. 

The  mean-square  pressure  observed  at  a  given  location  as  the 
result  of  several  sources  similarly  is  equal  to  the  sum  of  the 
individual  mean-square  contributions  (overall,  or  in  specified 
bands),  provided  that  the  various  individual  contributions  are 
uncorrelated  (i.e,,  not  substantially  in  phase  with  each  other)? 

*Fcr  correlated  signals,  the  total  roof -mean-square  pressure  is 
equal  to  the  sum  of  the  roof-mean-square  components. 


Since  it  is  unlikely  that  the  noise  components  produced  by  dif¬ 
ferent  sources  on  an  aircraft  will  be  correlated,  one  may  consid¬ 
er  them  generally  uncorrelated  for  all  practical  purposes. 


It  is  usually  convenient  to  work  only  in  terms  of  the  levels, 
and  to  avoid  converting  from  levels  back  to  the  basic  quantities, 
combining  these  quantities  as  outlined  above,  and  then  converting 
the  results  again  into  levels.  One  may  obtain  the  (power  or  pres¬ 
sure)  level  Lcomb  which  corresponds  to  the  combination  of  two 
components,  one  at  level  L,  and  one  at  L2  ,  by  using  the  chart  of 
Fig.  1.  For  example,  by  use  of  this  chart  one  may  determine  very 
easily  that  a  signal  of  90  dB  and  one  of  96  dB  combine  to  produce 
a  signal  of  97  dB  ( not  186  dB!). 


One  may  obtain  the  level  corresponding  to  a  combination  of 
several  signals  by  using  the  chart  repetitively  for  two  components 
at  a  time;  the  sequence  of  forming  the  two-at-a-time  combinations 
is  immaterial.  For  example,  consider  four  levels:  Li  =  60  dB, 

L2  =  65  dB,  L3  =  66  dB,  L4  =  68  dB.  Then  one  finds  L1+2  =  66.2  dB, 
L3  +  4  =  70.1  dB,  and  L(  !  +  2 )  +  (  j+i. )  *  71.6  dB;  or 

L(i+2)+3  =  69*°  dB*  L[(i+2)+3]+4  =  71*6  dB‘ 


DETECTION  RANGE  CALCULATION 


GENERAL  APPROACH 

The  "aural  detection  range"  of  a  given  aircraft  operating  at 
specified  conditions  may  be  defined  as  the  greatest  distance  at 
which  an  "average"  listener  on  the  ground  can  distinguish  the  air 
craft  sound  from  his  ambient  background  noise  environment.  In 
order  to  estimate  this  detection  range,  one  may  proceed  as  follow 

1.  Determine  the  frequency  distribution  of  the  acoustic  power 
emitted  by  the  aircraft,  due  to  all  sources; 

2.  Determine  the  "detection  level  spectrum",  i.e.,  the  sound 
pressure  level  spectrum  of  signals  that  the  average  listen¬ 
er  can  just  detect  in  the  presence  of  the  ambient  noise;  • 

3.  Calculate  the  distance  within  which  the  aircraft  sound  (at 
each  frequency)  reaches  the  level  at  which  it  can  be  de¬ 
tected  by  the  listener.  This  process  results  in  a  frequency 
distribution  of  ranges,  or  a  "detection  range  spectrum".  The 
maximum  of  this  spectrum  then  is  the  detection  range. 

This  approach  is  illustrated  graphically  in  Fir.  2,  to  which  fur¬ 
ther  reference  will  be  made  below. 

ACOUSTIC  POWER  EMITTED  BY  AIRCRAFT 

In  order  to  determine  the  spectrum  of  the  acoustic  rower 
generated  by  a  given  aircraft  under  specified  operating  condi¬ 
tions,  one  generally  may  proceed  best  by  considering  each  noise 
source  separately,  and  then  combining  the  effects  of  all  such 
sources,  Where  acoustic  data  directly  pertinent  to  specific 
components  (e.g.,  engines,  propellers,  or  engine-prorel ler  combi¬ 
nations)  are  available,  this  data  should  be  used.  Jf  such  data 
are  not  in  hand,  one  may  use  the  estimation  technique?  presented 
in'  the  appendixes  to  this  report. 

"’or  propeller- drive  aircraft ,  one  should  consider  the  fol¬ 
lowing  sources  (listed  in  the  rrobable  order  of  import ance ) * : 


#Note:  The  capital  letters  preceding  each  noise  component  listed 
indicate  the  appendix  which  may  be  used  to  estimate  the  noise 
from  that  component. 


1.  Proneller 

A  Rotational  (pure-tone)  noise 
B  Vortex  (broad-band)  noise 

2.  Engine  exhaust,  intake,  ana  casing 
C  Piston  engines 

D  Rotating  combustion  engines 
E  Turboshaft  engines 

3.  Power  transmission  systems 
J  Gears 

Aerodynamic  noise 
K  Boundary  layer  turbulence 
K  Movement  through  atmospheric  turbulence 
K  Turbulent  wakes 

5.  Vibrating  internal  components 

L  Radiation  from  point-excited  surface  structures 

For  jet-engine  drive  aircraft ,  one  should  consider  the  fol¬ 
lowing  sources  (again  listed  in  the  probable  order  of  importance)*: 

1.  Propulsion  system 

E  Turbine  exhaust,  inlet,  casing 
F  Pan  and  compressor  (nure-tone)  noise 
0  Jet  exhaust  flow 

H  Effects  of  inlet  and  exhaust  ducts 

2.  Aerodynamic  noise 

K  Roundarv  layer  turbulence 
K  Movement  through  atmospheric  turbulence 
K  Turbulent  wakes 

3.  Vibrating  internal  components 

L  Radiation  from  point-excited  surface  structures 

The  upper  part  of  Fig.  2  illustrates  the  acoustic  power 
spectra  for  an  hypothetical  prone ller-dri ven  aircraft.  For  the 
sake  of  clarity,  only  two  broad-band  source-spectra  (propeller 
vortex  noise  and  engine  casing  noise**)  and  pure  tones  from  only 
one  source  (propeller  rotational  noise)  are  shown.  The  figure 
also  shows  the  total  broad-band  noise  spectrum,  obtained  by  com¬ 
bining  the  two  contributions  according  to  Fig.  1. 


:iNote :  The  capital  letters  preceding  each  noise  component  listed 
indicate  the  appendix  which  may  be  used  to  estimate  the  noise 
from  that  component. 

**For  tills  hypothetical  aircraft,  the  exhaust  was  assumed  to  be  well 
muffled,  so  that  exhaust  noise  is  effectively  absent. 
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DETECTION  LEVEL  SPECTRUM 


As  discussed  in  Appendix  N,  the  sound  pressure  level  of  an 
acoustic  signal  that  a  listener  can  just  detect  depends  on  the 
listener’s  threshold  of  hearing  and  on  the  ambient  noise  that 
surrounds  him.  The  threshold  of  hearing  for  pure  tones  is 
slightly  different  from  that  for  bands  of  noise,  and  a  given 
(broad-band)  background  noise  generally  can  mask  a  broad-band 
signal  better  than  it  can  a  pure  tone.  Therefore,  one  must  treat 
pure-tone  and  broad-band  signals  separately. 

Appendix  N  presents  standard  hearing  threshold  curves  and 
indicates  how  one  may  determine  the  detection  levels  which  corre¬ 
spond  to  given  ambient  noise  spectra.  (The  detection  level  repre¬ 
sents  the  sound  pressure  level  of  a  signal  which  is  just  detect¬ 
able  against  the  background  noise.)  A  signal  is  detectable  if 
it  exceeds  at  any  frequency  both  the  hearing  threshold  level  and 
the  detection  level  -  thus,  if  it  exceeds  the  higher  of  these  two 
levels.  A  curve  which  at  any  frequency  follows  either  the  hear¬ 
ing  threshold  level  curve  or  the  detection  level  curve,  whichever 
is  higher,  thus  constitutes  a  "detection  level  spectrum",  which 
a  signal  must  exceed  to  be  detectable. 

Information  on  the  ambient  noise  in  some  environments 
likely  to  be  of  interest  in  relation  to  quiet  aircraft  design 
is  also  presented  in  Appendix  N,  in  terms  of  detection  levels  for 
pure-tone  signals. 

Figure  2  illustrates  how  one  may  derive  pure-tone  and  octave- 
band  detection  level  spectra  from  the  hearing  threshold  curves 
and  detection  levels.  The  hearing  threshold  curves  in  Fig.  2 
have  been  taken  from  Fig.  Nl.  The  pure  tone  detection  level  cor¬ 
responds  to  nighttime  data  for  a  dense  Panama  jungle  vfrom  Fig. 

N4);  the  octave-band  detection  data  has  been  derived  from  the  pure- 
tone  data  with  the  aid  of  Fig.  N2.  The  heavy  solid  and  dashed 
curves  then  correspond  to  the  detection  level  spectra,  for  pure- 
tones  and  octave-band  noise,  respectively. 

DETECTION  RANGE  SPECTRUM 

If  one  knows  the  spectrum  of  acoustic  power  levels  Lw  of 
the  noise  of  an  aircraft  and  also  the  spectrum  of  detection  levels 
Ld  applicable  to  given  conditions,  one  may  easily  determine  the 
distance  over  which  each  frequency-component  of  the  aircraft 
noise  must  travel  (in  an  ideal  non-absorbing  and  non-scattering 
atmosphere)  so  that  the  sound  pressure  level  of  the  signal  is 
eaual  to  the  detection  level.  In  view  of  Eq.  (3),  this  distance 
Ru  (ft)  may  be  found  from 


1  0 


(6) 


20  log  Ru  »  Lw  -  Ld  -  0.5  *  Lw  -  Ld 

Lw  and  L(j  are  expressed  in  decibels,  referred  to  the  standard 
reference  power  and  pressure,  respectively.  Because  of  the  un¬ 
certainties  usually  Involved  in  estimating  Lw  and  L^,  one  may 
usually  neglect  the  0.5  dB  correction. 

By  applying  Eq.  (6)  to  each  frequency  component,  one  may 
arrive  at  a  spectrum  of  "uncorrected  detection  ranges"  Ru.  The 
designation  "uncorrected"  implies  that  corrections  for  atmospher¬ 
ic  and  terrain  attenuation  effects  have  not  yet  been  included. 

If  Lw  and  are  plotted  to  the  same  scale,  as  in  Fig.  2, 
one  may  determine  20  log  Ru  graphically,  as  shown  in  the  figure. 
Again,  one  must  carry  out  separate  calculations  for  the  pure-tone 
and  for  the  broad-band  noise  components.  The  uncorrected  detec¬ 
tion  range  spectrum  one  obtains  from  Fig.  2  is  shown  as  the  dot¬ 
ted  curve  of  Fig.  3  for  broad-band  noise  and  as  the  open  circles 
for  pure-tone  noise. * 

From  the  information  presented  in  Figs.  Ml  though  M3  (of 
Appendix  M)  one  may  estimate  the  total  absorption  coefficient 
atot(dB/ft)  due  to  atmospheric  and  terrain  effects  at  each  fre¬ 
quency.  From  Fig.  MA  or  the  corresponding  relation  (see  last 
footnote  of  Appendix  M)  one  may  then  correct  the  previously  found 
"uncorrected  detection  range"  spectrum,  to  obtain  a  spectrum  of 
actual  (corrected)  detection  range.  The  maximum  of  this  spectrum 
determines  the  range  at  which  the  aircraft  wil]  first  be  heard 
under  the  specified  conditions;  by  noting  the  frequency  at  which 
this  maximum  is  found  and  ref'eriing  back  to  the  noise  source 
spectra,  one  may  then  determine  which  noise  component  is  primarily 
responsible  for  earliest  aural  detection  of  the  aircraft  -  i.e., 
which  noise  component  must  be  reduced  to  reduce  the  detection 
range  significantly. 

The  solid  curve  and  filled-in  circles  of  Fig.  3  represent 
the  corrected  detection  range  spectrum  for  the  octave-band  and 
pure-tone  noise  components,  respectively,  as  obtained  from  the 
uncorrected  spectrum  for  certain  conditions  (atot  a  0.02  dB/ft). 

One  finds  that  the  maximum  of  the  corrected  range  spectrum  cor¬ 
responds  to  20  log  Rc*  70  indicating  that  the  aircraft  will  first 
be  detected  aurally  at  a  range  of  about  3000  ft.  One  finds  that 
the  160  Hz  pure  tone  component  of  propeller  rotational  noise  will 
be  heard  first.  If  this  component  can  be  reduced  by  a  little 
more  than  3  dB,  then  the  80  Hz  tone  becomes  predominant,  and  the 
detection  range  is  reduced  to  about  2000  ft.  If  all  the  propeller 
pure-tones  could  be  eliminated  (or  reduced  by  10  dB  or  more),  the 
detection  range  would  be  about  1000  ft;  with  the  noise  components 

»Aopendix  P  presents  an  alternate  method  for  determining  Ru.  This 
alternate  method  is  oreferable  if  recorded  flyby  noise  data  are 
available,  or  if  there  exists  considerable  uncertainty  about  whether 
a  noise  component  should  be  considered  as  a  pure-tone  or  broadband. 
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between  125  and  160  Hz  being  heard  first.  One  notes  from  Pig. 
that  at  these  frequencies  the  noise  is  primarily  engine  casing 
noise;  therefore,  one  would  need  to  reduce  the  noise  from  this 
source  in  order  to  achieve  any  further  reduction  in  detection 
range . 
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FIGURE  1  CHART  FOR  COMBINING  TWO  UNCORRELATED  ACOUSTIC  LEVELS, 
L]  and  Lg.  LEVELS  MAY  BE  POWER  LEVELS  OR  SOUND- 

PRESSURE  LEVELS 


i 
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FIGURE  2  ILLUSTRATION  OF  DETECTION  RANGE  CALCULATION  APPROACH 


DETECTION  RANGE  SPECTRA  CORRESPONDING  TO  FIGURE 


APPENDIX  A 


ROTATIONAL  NOISE  OF  PROPELLERS 


DIPOLE  SOUND 

Far  ly  investigations  of  propeller  noise  (Ref.  L )  showed  that 
the  lift  and  drag  forces  acting  on  the  propeller  blades  give  rise 
to  equal  and  opposite  inertia  forces  on  the  fluid  medium,  and 
these  forces  constitute  arrays  of  acoustic  dipole  sources  distri¬ 
buted  over  the  propeller  disk  and  rotating  in  the  propeller  plane 
at  the  propeller's  angular  velocity. 

In  order  to  analyze  the  acoustic  effects  of  these  dipole 
sources,  it  is  convenient  to  treat  first  the  case  of  a  propeller 
having  B  blades,  rotating  at  angular  velocity  ft,  (but  not  advanc¬ 
ing  axially)  with  uniform  inflow  to  the  propeller.  An  element 
of  the  fluid  medium  at  polar  coordinate  position  r^pj  in  the 
plane  of  the  propeller  is  subjected  to  forces  Fg  and  as  the 
propeller  passes  the  element's  location,  which  correspond  to 
infinitesimal  increments  of  the  thrust  and  the  drag  on  the  pro¬ 
peller  blade.  When  there  is  no  propeller  blade  present  at  the 
position  Tjjifj,  there  is  no  force  exerted  on  the  fluid. 

The  duration  of  each  force  pulse  is  determined  by  the  local 
width  of  the  propeller  blade,  projected  onto  the  propeller  disk 
plane,  and  by  the  propeller'  rotational  speed.  The  force  on 
the  rluid  at  a  typical  location  then  varies  with  time  as  shown 
in  Fig.  Al.  If  one  expands  such  a  periodic  function  in  a  Fourier 
series,  one  obtains  a  sum  of  harmonic  functions  whose  frequencies 
are  multiples  of  the  blade  passage  frequency  Wg  =  ftB. 

One  may  express  (see  Ref.  2)  the  axial  force  distribution  as 

Fp,  ( t ;  r  j  ,  j  )  =  £  Asexp[-iswB(t-4»1/n)] 

s 

oo 

=  Frj^dgexpdsB^jJexpfisBftt)  (1) 

_00 

where  As  and  cij  are  Fourier  coefficients;  s  is  the  harmonic  num¬ 
ber  index;  F8  is  the  thrust  force  amplitude. 


Similarly,  the  drag  force  distribution  may  be  expressed  as 


exp(isB(J>,  )exp(-isBfit) 

^  —  CD 


(2) 


As  shown  in  Ref.  3,  application  of  a  concentrated  force  at  a 
point  in  a  fluid  is  equivalent  to  locating  an  acoustic  dipole 
source  at  that  point,  with  the  dipole  axis  in  the  direction  of 
the  force.  Thus,  the  double  array  of  forces  over  the  propeller 
disk  due  to  thrust  and  drag  gives  rise  to  acoustic  radiation  from 
a  double  array  of  dipole  sources,  the  thrust-related  array  having 
all  sources  aligned  parallel  to  the  propeller  shaft  and  the  drag- 
induced  sources  aligned  parallel  to  the  propeller  plane. 

To  find  the  sound  field  produced  by  the  previously  given 
force  distributions,  one  may  use  the  Green’s  function  for  a  di¬ 
pole  source  of  unit  strength  in  the  z-direction  and  of  sinusoidal 
frequency  to,  (Ref.  2) 


G 

i  z 


ikr  +»  im(  $-<}>.  )  .  „ 

ikcosv  'jj—  £  1  J  (kr1sinv)e  e“la,t 


(3) 


where  the  distances  r  and  r,  and  angles  v,<J>  and  <Ji  are  defined 
In  ^ig.  A2,  and  where  Jm  renresents  the  m'th  order  Bessel  function, 
Also,  k  =  w/c  is  the  wavenumber  and  c  represents  the  wave  propa¬ 
gation  speed,  i.e.  the  speed  of  sound.  The  Green's  function  for 
a  unit  dipole  in  the  <!> t -direction  is  similarly  given  by 


1  ^  ~  in ( )  .  .. 

n  -  1 _  _ _ V  n,v,  i  -lwt 


Enl'j  ( kr  sinv)e 
„  m  i  • 


The  foregoing  forms  of  the  Green's  functions  involve  approxi¬ 
mations  which  are  valid  only  if  the  point  of  observation  is  in  the 
geometric,  as  well  as  in  the  acoustic  far-field  of  the  source,  - 
i.e.*  if  the  observation  distance  is  much  greater  than  both  the 
physical  dimensions  of  propeller  disk  and  the  acoustic  wavelength 
in  the  fluid  medium. 

The  sound  pressure  due  to  the  thrust  force  then  may  be  de¬ 
termined  from 

Pz(r,*,v,t)  =//Fz(t;rl,4>1)Glz(r1,<t)1;r,(j),v)rldr1d<()1  (5) 
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and  that  due  to  the  drag  force  from 


P<p(r,<l>,v,t)  =/fP<fl(t;rl,<l>l)ai(l)(rl)4>1;r,<t>,v)r1drld<p1  (6) 

By  substituting  Kqs.  (1)  and  (3)  into  (5)  and  noting  that  the 
frequency  in  the  Green's  function  must  be  the  same  as  that  of  the 

r’  RQ 

source  component  ( co  =  sBfi,  k  =  =  sk  t ) ,  one  obtains 


ik  cosv 

Pz<r,«,v,t)  -//  -  -  I-,- 


«<s.imJn(sk1  rt  siny  )e^sl<  i r 


im(q-4i1)  -isB»2t  is B<f> t 
e  e  e 


P°m,.2nr  jdrj  d(()  j 


Because  of  the  factors  e 


in( <!>-<(>  )  isBtj), 

?  and  e  ,  the  integration 


over  <f>j  results  in  zero  terms,  except  for  m  =  sB,  and  where  the 
integration  introduces  a  factor  of  2ir.  (since  the  radiation  func¬ 
tion  GlZ  is  built  up  of  components  periodic  in  ((>,  ,  as  measured 
from  some  reference  angle,  and  since  the  source  distribution  is 
also  composed  of  terms  periodic  in  <j>j  ,  only  those  radiation  func¬ 
tion  harmonics  which  are  spatially  in  phase  with  the  source  dis¬ 
tribution  terms  contribute  to  the  sound  pressure  disturbance 
arriving  at  the  point  of  observation.)  The  summation  with  re¬ 
spect  to  m  and  the  integration  over  the  angle  (}>  can  thus  be 
handled  quite  simnly,  and  F,q.  (7)  nay  be  reduced  to 


-i k i cosv  r  isk.r 

7, nr  51  e  e  8 


J*P”nto  J,.^(  sk ,  r  j  sinv  )2mr  jdr , 


Similarly,  the  sound  pressure  due  to  the  torque-related  forces 
is.  ^ound  to  be 


z  .  4  \  Bi  .  s  13  ^  *  isB(f>  —  isB12tf  1  nio  , 

P/r,-q,v,t)  =  777  Xa  1  c  p  G  J  F.^VsB 


( sk  j r  sinv )2ur 


idri  | 
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This  elementary  theory  based  on  the  rotating  forces  exerted 
on  the  fluid  medium  by  the  propeller  blades  predicts  two  series 
of  discrete-frequency  pressure  disturbances;  each  series  is  com¬ 
posed  of  multiples  of  the  blade  passage  frequency,  and  each  pres¬ 
sure  term  varies  Inversely  with  the  first  power  of  the  distance 
from  the  propeller  hub  to  the  observation  point.  Each  series 
term  i  .volves  a  directivity  factor,  and  the  source  level  of  each 
term  is  related  to  the  aerodynamic  harmonic  forcing  function 
which  produces  the  dipole  source.  From  the  exponential  factors 
in  each  series  representation,  one  may  note  that  at  a  fixed  value 
of  r,  the  angular  position  4>  of  constant  phase  rotates  with  the 
shaft  velocity  ft. 

Since  the  factors  F°,  F^,  and  as  in  general  are  functions 
of  radial  position  rj ,  evaluation  of  the  integrals  in  Eqs.  (8) 
and  (9)  presents  some  difficulties.  Some  approximations  are  pos¬ 
sible,  as  presented  below.  By  use  of  the  mean-value  theorem,  one 
may  write 


J  P°VsB(skirisinv)27rridri  = 

asJsB(skir\sinv)lr  =a  J  Fl27tridri  =  asJsB(skirisinv)  r=a  ,Ft 

1  Z  *  Z 


where  a  denotes  the  propeller  radius,  az  is  an  "effective  radius" 
at  which  asJsR  is  to  be  evaluated  in  order  to  produce  the  correct 
value  of  the  integral  ,  and  Ft  represents  the  total  thrust  force 
exerted  on  the  propeller.  Equation  (8)  then  may  be  reduced  to 


PzO,4>,v,t) 


and,  similarly,  Eq.  (9)  may  be  rewritten  as 
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whore  a*  is  the  effective  radius  for  drag,  analogous  to  az  for 
thrust,  and  T  represents  the  total  drag  torque.  The  values  of 
az  and  a^  depend  upon  both  the  radial  variation  in  loading  repre¬ 
sented  by  F°  and  and  the  angular  variation  of  blade  forces 
represented^by  the^factor  ;  since  the  value  of  a§  is  different 
for  each  frequency  contribution,  the  effective  radii  az  and  a^, 
generally  also  depend  on  the  harmonic  index  s. 

Morse  and  Ingard  (Ref.  2)  introduce  the  approximation  that 
az  =  =  a^  so  that  the  total  sound  pressure  due  to  the  two  di¬ 

pole  source  mechanisms  is  given  by  the  sum 
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From  this  r  suit  one  finds  that  the  amplitude  of  the  sth  harmonic 
frequency  ;  mponent  is  given  by 


(  a  )  r.Bft 

„  =  _  3 _ 

f  s  Tirrc 
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t  ccn;v  +  Wj  JsB(5BMd  sinv) 
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where  repi  esents  the  rotational  ?lach  number  at  the  effective 
radius  a^, 

a 

M  =  -1- 
d  c 


and  where  =  T/a^  represents  the  drag  force. 

liquation  (13)  gives  the  amplitude  of  propeller  noise  due  to 
dipole  source  mechanisms  associated  with  the  aerodynamic  force 
loadings  on  the  propeller  blades  which  give  rise  to  thrust  and 
rotational  drag.  By  virtue  of  the  cosv  term,  the  thrust  contri¬ 
bution  vanishes  ,J.n  the  plane  of  the  propeller;  the  sum  of  the 
two  contributions  vanishes  on  the  propeller  axis,  in  view  of  the 
Bess, el  function  directivity  factor.  [Equation  (13)  may  be  seen 
to  be  a  generalization  of  Tutin's  formulation.] 
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MONOPOLE  SOUND 


In  addition  to  providing  net  lift  and  drag,  a  propeller  blade 
element  displaces  some  fluid  by  virtue  of  its  thickness,  and  thus 
acts  as  a  monopole  sound  source.  In  order  to  analyze  this  sound 
contribution,  one  may  describe  the  propeller  in  terms  of  a  blade 
thickness  function  Sj(rJ,$}),  which  defines  the  blade  thickness 
distribution  in  polar  coordinates  r  J  ,4> |  rotating  with  the  propeller 
For  a  fixed  radial  distance  r*  in  the  propeller  plane,  the  thick¬ 
ness  function  may  be  expressed  as 


vr; 


)  = 


D  ( r '  )  S 
i  1  1 


Up  = 


°i(rp 


+  CO 

* 


UB*{ 


(14) 


where  D,(rJ)  gives  the  maximum  blade  thickness  as  a  function  of 
radial  distance,  and  S  Up  and  its  Fourier  expansion  describes 
the  blade  cross-section  profile.  In  terms  of  the  nonrotating  co¬ 
ordinates  Tj  =  rj  and  <J>  t  =  4>  J  +  fit,  the  thickness  function  becomes 
3  ( r  j  ,<f)  j-flt )  and  has  the  form  sketched  in  Fig.  A3. 

The  monopole  source  term  in  the  general  wave  equation  is  the 
time  derivative  of  the  rate  Q  at  which  mass  is  displaced  per  unit 
area  of  the  propeller  plane.  Here  one  finds  that  Q  is  given  by 

=  o(H)  =  -»«s'  ■  s  <15) 

where  p  denotes  the  density  of  air,  so  that  the  acoustic  source 
term  is 
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Using  Eq .  (14),  one  may  rewrite  this  term  as 


f?  (VV  =  _pn2D 


+  00 

,(v  .£ 


t=-<» 


‘  B 


UB*| 


(l/) 


To  obtain  the  acoustic  pressure  at  the  observation  point 
(r,4>,v),  one  may  form  the  product  of  the  source  function  and 


Green's  function  for  a  monopole  source,  and  integrate  this  pro¬ 
duct  over  the  source  plane.  After  some  manipulation,  one  arrives 
at 
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As  before,  one  is  faced  with  evaluating  an  integral  containing 
the  product  of  a  Bessel  function  and  of  two  possibly  complicated 
functions  of  the  variable  of  integration.  One  may  again  apply 
the  mean-value  theorem  and  express  the  integral  as 


/  Dj  (r,  )6JlJiB(Jlk1r1  sinv  )2*r  x  dr  , 
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where  VD  denotes  the  volume  of  a  solid  disk  having  the  radius  a 
of  the  propeller  and  the  same  radially-varying  thickness  as  a 
propeller  blade,  and  am  represents  an  appropriate  mean  radius. 

Thus,  the  sound  pressure  amplitude  of  the  sth  harmonic  com¬ 
ponent  due  to  the  blade  thickness  (monopole  source  effect)  is 
given  by 
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One  may  observe  that  here  the  directivity  (the  dependence  on  v) 

1.;  the  same  as  for  combined  dipole  sources,  see  Eq.  (13)*  Also, 
the  monopole  term  of  Eq.  (20)  varies  inversely  with  distance,  just 
like  the  dipole  term  of  Eq .  (13).  Comparing  the  monopole  and  di~ 
pc> ]  a  source  terms  further,  one  find.;  that  the  monopole  term  varies 
a.;  the  square  of  sBft,  whereas  the  thrust-dipole  term  involves  the 
first  power  of  sBfi  and  the  drag-dipole  term  involves  the  first 


power  of  sBfi,  modified  by  the  rotational  Mach  number  M(j,  AIbo, 
one  may  observe  that  the  dipole  and  monopole  termB  are  90°  out 
of  phase. 

Effects  of  Disturbed  Inflow  and  Forward  Motion 

The  previous  analysis  has  dealt  only  with  the  simple  ease  of 
a  rotating  propeller  with  no  forward  motion  and  no  nonuniformities 
in  the  inflow. 

The  effects  of  nonuniform  flow  are  considered,  for  example, 
in  the  theoretical  analysis  of  Morse  and  Ingard  (Ref.  2)  and  have 
been  given  extensive  attention  in  studies  of  compressor  noise, 
(e.g.,  Refs.  4,5,6)  where  the  presence  of  struts,  stators,  etc. 
ahead  of  the  rotating  blades  serves  to  disturb  inflow  and  thus 
produces  perturbations  of  the  aerodynamic  loadings  on  the  propel¬ 
ler  blades.  For  small  propeller-driven  aircraft  designed  for 
quiet  operation,  this  source  of  noise  is  likely  to  be  unimportant 
and  is  not  discussed  here. 

The  effects  of  motion  of  the  acoustic , source (s )  -  that  is, 
of  forward  motion  of  the  aircraft  -  have  also  been  studied  (Refs. 
2,7, 8) . Primarily ,  there  result  the  well-known  Doppler  shift  of 
observed  frequencies  and  changes  in  the  directivity.  Both  of 
these  phenomena  are  amenable  to  straightforward  analysis,  but  are 
of  little  interest  here,  and  therefore  are  not  treated  further. 

Comparison  of  Theoretical  Predictions  with 

Measured  Data  < 

Comparison  of  theory  with  experiment  is  made  difficult  by 
the  complexity  of  the  calculations  (particularly  for  the  higher 
harmonics),  by  the  approximations  involved  in  the  derivations,  by 
uncertainties  in  the  experimental  conditions,  and  by  lack  of  suf¬ 
ficiently  detailed  data,  Since  the  concern  here  is  with  discrete 
frequency  noise,  one  needs  narrow-band  data  for  comparison,  or  at 
least  an  assessment  that  the  data  in  specific  bands  are  due  pre¬ 
dominantly  to  pure-tonr  impeller  contributions.  Because  of  the 
significant  dlrectivlt,,  effects,  one  also  needs  to  know  the 
source-receiver  orientation  rather  precisely.  Nearly  all  of  the 
useful  data  presently  available  is  from  stationary  operation 
(ground  run-up),  with  measurements  taken  at  several  azimuthal 
positions  around  the  airplane.  (Comparison  between  theory  and 
measurements  usually  is  most  convenient  for  positions  in  the 
plane  of  the  propeller;  however,  the  thrust-related  dipole  term 
vanishes  here,  and  its  contribution  'can  then  not  be  determined.) 

Figures  A4  and  A5  compare  the  results  of  calculations  with 
experimental  measurements  for  the  three-bladed  OV-1  propeller. 


Both  figures  pertain  to  the  same  observation  position  givai  from  the 
propeller  axis,  in  the  plane  of  the  propeller,  but  the  former 
figure  corresponds  to  a  lower  rotational  speed  and  shaft  horse¬ 
power  than  the  latter.  Dipole  sound  prediction  calculations  were 
carried  out  on  the  basis  of  Eq.  (13) >  taking  2as  =  1,  which  corre¬ 
sponds  to  a  uniform  distribution  of  aerodynamic  forces  over  the 
blade  width.  Results  are  shown  for  two  values  of  a^:  70 %  and 

80%  of  the  propeller  radius  a;  one  notes  that  these  two  results 
agree  well  for  the  fundamental,  but  diverge  increasingly  for  the 
higher  harmonics  —  indicating  that  the  choice  of  a<j  affects  the 
higher  harmonics  more  strongly. 

Also  Indicated  in  Figs.  A1*  and  A5  are  the  sound  pressure  levels 
predicted  for  this  same  propeller  due  to  monopole  radiation.  For 
the  prediction  calculations  corresponding  to  Fig.  A**  the  disk 
volume  Vp  was  approximated  from  a  graph  of  blade  thickness  versus 
radius.  uFor  the  lowest  two  harmonics,  the  calculation  was  car¬ 
ried  out  on  the  basis  of  Eq .  (18),  and  by  accomplishing  the  inte¬ 
gration  by  means  of  a  series  expansion,  with  only  the  largest 
terms  retained.  For  these  harmonics,  the  effective  radius  a™  was 
found  to  be  0.57a  and  0.6^,  respectively.  The  third  harmonic 
contribution  was  estimated  by  use  of  Eq .  (20),  with  an  effective 
radius  am  =  0.75a.  The  values  plotted  in  the  figure  are  based 
on  the  assumption  of  a  rectangular  blade  cross-section  profile 
(which  assumption  affects  the  Fourier  coefficients  8S).*  Since 
the  effective  radius  and  the  Fourier  coefficients  to  be  used  in 
the  approximate  Eq .  (20)  are  related  to  the  propeller  geometry, 
and  not  the  operating  conditions,  it  seems  valid  to  use  values  of 
these  parameters  which  yield  acceptable  results  at  one  set  of 
operating  conditions  to  predict  noise  levels  at  other  operating 
conditions  of  the  same  propeller  by  use  of  the  simpler  Eq .  (20). 
Figure  5  shows  the  monopole  sound  pressure  levels  calculated  in 
that  manner,  assuming  a  rectangular  blade  profile  and  a  =  0.55a 
for  s  =  1,2,  and  0.75»  for  s  =  3- 


*Levels  for  the  first  two  harmonics  were  also  calculated  on  the 
basis  of  a  triangular  variation  of  blade  thickness  with  angle  <f>  J  . 
At  the  fundamental,  the  triangular-  thickness  gave  nearly  the 
same  result  as  the  rectangular  thickness,  but  at  twice  that  fre¬ 
quency,  the  triangular  distribution  gave  a  level  which  was  lower 
by  nearly  10  dB.  By  the  nature  of  the  Fourier  series  represen¬ 
tations,  one  would  expect  the  difference  between  levels  to  in¬ 
crease  rapidly  with  higher  harmonics.  As  another  example  of  the 
sensitivity  of  the  theoretical  model  to  small  variations  in  the 
assumptions  necessary  for  numerical  evaluation,  the  level  of  the 
third  harmonic  was  calculated  using  a  rectangular  thickness  pro¬ 
file,  but  at  an  effective  radius  of  0.60a.  This  result  was 
found  to  be  about  13  dB  lower  than  the  value  obtained  with 


Combination  of  sound  levels  due  to  monopole  and  diple  sources 
is  accomplished  by  the  usual  method  of  adding  uncorrelated  sound 
levels,  since  the  pressures  due  to  the  two  sources  are  90°  out  of 
phase.  Results  of  combination  of  the  monopole  terms  with  the  di¬ 
pole  terms  for  a^  -  0.8a  are  shown  in  Figs.  A4  and  A5  as  well  as 

corresponding  measured  data.*  [It  should  be  noted  that  the 
measured  levels  of  Fig.  A5  are  obtained  from  octave-band,  not 
narrow-band  data;  however,  narrow-band  data  suggest  that  the 
octave-band  levels  are  predominantly  due  to  the  discrete  frequency 
contributions . 3  The  agreement  between  predicted  and  measured 
levels  rs  quite  good.  One  notes  that  although  the  dipole  term 
alone  predicts  the  first  harmonic  quite  accurately  and  approximates 
the  second  harmonic  reasonably,  inclusion  of  the  monopole  con¬ 
tribution  greatly  improves  the  prediction  for  the  third  harmonic. 


Figure  A6  presents  a  comparison  between  predictions  and  data 
for  another  propeller,  and  also  indicates  predictions  computed  by 
methods  from  Hubbard  (Ref.  9)  and  Dodd  (Ref.  10)  for  the  dipole 
and  dipole-plus  monopole  sources.  Insufficient  information  re¬ 
garding  the  propeller  geometry  is  available  to  allow  prediction 
of  the  blade  thickness  pressure  term  by  the  method  of  Eq.  (20) . 

The  measured  data  shown  in  Fig.  A6  are  seen  to  agree  reasonably 
well  with  the  levels  predicted  here  by  use  of  only  the  dipole 
term  based  on  ad  -  0.8a.  Estimates  of  the  monopole  source 

contribution  shown  in  Fig.  A6  have  been  prepared  on  the  basis  of 
the  assumption  that  the  propeller  of  Fig.  A6  is  geometrically 
similar  to  that  of  Figs.  A4  and  A5 . 

In  order  to  assess  how  well  one  may  predict  the  thrust  di¬ 
pole  contribution,  calculations  have  been  carried  out  for  observer 
positions  out  of  the  planes  of  the  two  previously  discussed  pro¬ 
pellers.  Calculated  levels  and  measured  data  are  compared  in 
Figs.  A7  and  A8.  For  the  U-10  propeller  (Fig  8)  agreement  between 
estimated  and  measured  values  is  seen  to  be  quite  acceptable  for 
positions  30°  ahead  and  behind  the  propeller  plane,  but  very  poor 
for  60°  behind  the  plane  of  the  propeller,  particularly  in  the 
second  and  third  harmonics.  The  agreement  for  the  OV-1  propeller 
(Fig,  A7  is  generally  poor.  Brief  further  study  shows  that  the 
high-frequency  predictions  are  even  worse  for  angles  closer  to 
the  axis  of  rotation;  in  particular,  the  data  suggests  the 
presence  of  strong  tonal  contributions  in  the  axial  direction, 
which  the  theory  predicts  to  vanish. 


*In  order  to  make  the  "f ree-f ield"  predictions  correspond  to  the 
measurements  made  near  the  ground,  6  dB  has  been  added  to  all 
predicted  values,  in  order  to  account  for  ground  reflection. 


What  physical  arguments  can  one  advance  regarding  the  dif¬ 
ference  between  the  predicted  and  actual  behaviors?  (1)  Non- 
uniform  flow:  It  seems  unlikely  that  parts  of  the  aircraft  up¬ 
stream  from  the  propeller  could  disturb  the  inflow  seriously.  It 
is  of  interest,  however,  that  the  aircraft  for  which  best  agree¬ 
ment  was  obtained  at  30°  out  of  the  propeller  plane  has  a  single 
engine,  whereas  the  other  craft  has  two  wing-mounted  engines. 
Conceivably,  the  forebody  of  the  airplane  might  disturb  the  inflow 
in  the  latter  case,  but  it  seems  doubtful  that  the  effect  would 
be  as  strong  as  observed.  (2)  Reflection  and  radiation  from  the 
ground,  and  from  wings,  fuselage,  etc. :  Again,  these  effects 
appear  incapable  of  producing  the  large  observed  difference  be¬ 
tween  calculated  and  measured  levels;*  radiation  of  sound  from 
airplane  surfaces  (e.g.,  due  to  structureborne  vibration)  cannot 
be  evaluated  from  the  available  data.  (3)  Approximations  in 
analytical  evaluation:  The  thrust  and  drag  lorces  on  the  propel¬ 
ler  blades  have  been  approximated  by  uniform  loading  distribu¬ 
tions,  and. in  the  evaluation  of  the  resulting  Fourier  coefficients 
the  propeller  solidity  was  taken  to  be  sufficiently  small  so  that 
essentially  point  loadings  result.  These  approximations  assure 
relatively  large  high-frequency  content  in  the  excitation  due  to 
dipole-related  forces.  Similarly,  expansion  of  the  blade  thick¬ 
ness  distribution  on  the  basis  of  a  uniform  cross-section  nrob- 
ably  overemphasizes  the  actual  high-frequency  spectral  content. 
Since  the  disagreement  between  theory  and  measurements  becomes 
greater  with  Increasing  frequency  at  positions  out  of  the  pro¬ 
peller  plane,  it  appears  that  the  above  approximations  are  not 
the  source  of  any  discrepancies. 

None  of  the  above  possibilities  presents  themselves  as  ob¬ 
vious  sources  of  discrepancies  between  theory  and  data;  deeper 
investigation  of  these  and  other  conceivable  mechanisms  is  re¬ 
quired,  but  is  beyond  the  scope  of  the  present  work. 

ESTIMATION  OF  TOTAL  ACOUSTIC  POWER 

One  may  determine  the  acoustic  power  radiated  by  the  dipole 
and  monopole  sources  by  integrating  acoustic  intensities  obtained 
from  Kqs .  (13)  and  (20)  over  a  sphere  of  radius  r  (i.e.,  over  all 
rinrie;'.  v)  and  adding  the  results.  However,  this  computation  may 
he  somewhat  complicated,  and  one  may  do  just  as  well  for  general 


* However,  ground  reflection  may  cause  the  directivities  measured 
with  the  aircraft  on  the  ground  to  differ  from  those  produced  by 
the  aircraft  in  flight. 


estimation  purposes  by  evaluating  the  pressures  at  an  "average” 
angle  v  =  45°  and  multiplying  the  resulting  average  intensity 
by  the  area  of  the  sphere. 


Thus,  one  may  estimate  the  acoustic  power  Wa  of  the  s th  har¬ 
monic  due  to  both  dipole  and  monopole  sources  from 
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where  pS2  represents  the  value  of  the  dipole-related  pressure,  as 
found  from  Eq.  (13)  for  v  =  45°  and  psi  represents  the  monopole- 
related  pressure,  as  obtained  from  Eq.  (20)  with  v  =  45°. 

ESTIMATION  SCHEME 


The  acoustic  power  W0  for  each  harmonic  component  must  be 
calculated  separately.  Generally  only  the  lowest  few  components 
(s  =  1,2,3)  are  of  interest. 

1.  Calculate  rp  from  the  modified*  Eq.  (13), 
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(13a) 


The  propeller  thrust  Ft  and  drag  torque  T  usually  are  avail¬ 
able  from  propulsion  system  characteristics.  Otherwise,  for 
a  propeller  with  radius  a,  constant  section  lift  and  drag 
coefficients  and  CD,  and  with  a  small  pitch  angle,  one 
may  use 


*Obtained  by  taking  v  =  4S°  and  2as  «  1.  This  approximation  for 
a3  apolies  for  uniform  distribution  of  the  aerodynamic  forces 
over  the  blade  width. 

Note  that  values  of  Jm(x)  are  available  in  standard  tables  of 
Bessel  functions  and  that  (see  Ref,  2) 
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for  x  >>  1 


Ft  .  0.17CLpU^ipAb 


T  *  0.13CDpU£ipAba 


where  Upip  is  the  blade  tip  speed  and  A^  the  total  blade 
area  (for  B  blades,  projected  onto  plane  of  rotation). 


2. 


Find  rpgl  from  the  modified  Eq.  (20), 
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where  one  may  find  the  propeller  disc  volume  Vp  by  integrat¬ 
ing  the  radial  thickness  distribution,  or  approximately  from 
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where  d  is  the  average  blade  thickness.  The  Fourier  coef¬ 
ficient  of  the  circumferential  blade  thickness  distribu¬ 
tion  may  be  calculated  from  a  known  thickness  function  S ( 4> ) 
for  the  propeller,  from 
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or  may  be  approximated  for  a  low  solidity-ratio  propeller  by 
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where  (b/a)ave  represents  the  average  propeller  chord  to 
radius/ratio.  Unless  better  information  is  available,  one 
may  take  a.,,  =  0,55a  for  s  =  1,  am  =  0.65a  for  s  =  2, 
am  =  0.75a‘  for  s  *  3. 

3.  Calculate  the  total  acoustic  power  W5  produced  at  frequency 
w  =  sB£2  from 
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(21a) 


pH 


i 


and  then  find  the  corresponding  power  level  from 


Lw  =  10  log(Ws/Wref> 

4.  Determine  the  frequency  u>s  at  which  the  s th  component  occurs 
from 


uj  (radians/sec)  =  sRft 
s 

or 

f  (Hz)  =  u)s/2tt  =  sBn/2Ti 

Illustrative  Calculation  ' 

Consider  a  three-bladed  propeller  of  4  ft  radius,  rotating 
at  2100  rpm  and  producing  a  thrust  of  340  lb  at  a  torque  of  350 
ft/lb,  at  1000  ft  altitude  (where  the  air  density  is  0.070  lb/ft3 
and  the  speed  of  sound  is  1100  ft/sec). 

To  use  Eq.  (13a),  find 

Jl  *  2100  rev/min  =  2100( 2n/60 )rad/sec  =  220  rad/sec 
ad  =  0.8a  =  0.8(4  ft)  =  3.2  ft 

=  adn/,c  =  (3.2  ft)  (220  rad/sec )/( 1100  ft/sec)  =  0.64 
BMd//?  =  3(0.64)//2  =  1.35 


BO 

T  Ft 

3(220  rad/sec) 

350  ft-lb  340  lb 

4"nc 

adMd  "  /2 

4ir(1100  ft/sec) 

13.2  ft)  (0.6*)  ft 

=  0.0525  |  170  -  240  |  =  3.7  lb/ft. 

Then, 

for*  s  =  i,  rnS2  =  3.7  J3(1.35)  =  .3 . 7  (4 . 5*10-2 )  =  0.17  lb/ft 

for  s  =  2,  rps'  =  2(3.7)  J6(2.70)  =  7 . 4  (6 . 5*10“  3  )  =  0.048  lb/ft 

for  s  =  3,  rp  2  =  3(3.7)  J9(4.05)  =11 . 1 ( 1.05*10-3 ) =  0.011  lb/ft 

S  2 


•Values  of  Bessel  functions  J3,J6  ,J9  found  from  Handbook  of 
Mathematical  Functions,  Ed.  by  M,  Abramowitz  and  I. A.  Stegun, 
National  Bureau  of  Standards,  Washington,  D.  C.  1964,  p.  398. 


To  use  Eq.  (20a),  one  needs  to  know  also  the  average  blade 
thickness  d  and  the  average  ratio  of  the  blade  chord  to  the  blade 
radius.  These  values,  which  may  be  obtained  from  a  drawing  of 
the  propeller,  here  are  taken  as  d  *  0.2  ft,  (b/a)  _  =  0.25. 

ctvg 

Then  find 


Tia2d  =  7r  ( 14  ft  )2  (0 . 2  ft)  =  10.1  ft3 


B 

7T 


(0.25)  * 


0.24 


p(RSi)2 

IHr 


V, 


— -  -  U3Uig07ra^^  (10tl  ft3)(0.2l|) 
32.2  lbm  ft/sec2lb 


=  180  lb/ft 


Bfia  =  3(220  rad/sec) (4  ft) 
c/2  /2('1100  ft/sec) 


and 


.  X 


for 

s  = 

1, 

rpSl 

for 

s  = 

2, 

rpSl 

for 

s  = 

3, 

rnsi 

1&0  J3[l. 70(0. 55)3  =  180  J j ( 0 ,9 3) 

l8o(.015)  =  2.7  lb/ft 

4(180)  J6  [2(1.70)  (0 .65)  3  -  720  J4(2.2) 

720(2. lxlO"3 )  =  1.5  lb/ft 

9(180)  J, C 3(1 . 70 ) (0 . 75 ) ]  =  1620  J9(3.8) 

1620(6.  lxio-:>)  =  1.0  lb/ft 


Since  pc  - 1  \  f  (1100)  =  2.4  lb  sec/ ft3 ,  so  that 


2tt  /p  c 

=  2.8  ft3 /lb  sec,  one  finds  from 

Eq.  (21a) 

that 

for 

s 

=  1,  Vs  -  2 . 8[ (0 . 17) 2  +  (2 . 7 ) 2 ]  = 

■  20.4 

ft 

lb /sec 

for 

s 

=  2,  W s  ■  2.8[(0.048)2  +  (1.5)2] 

=  6.3 

ft 

lb /sec 

for 

s 

=  3,  wc  =  2 . 8 [ (0 . Oil )2  +  (1.0)2] 

5 

=  2.8 

ft 

Ib/se  c 

Then  since  ^ref  a  0.7^  *  10“12  ft  lh/sec, 

for  s  -  l,  Lw  -  10  log  [20.V(0.74xlO-la)3  -  135  dB»re  10“12 

for  s  =  2,  Lw  -  10  log  T  6.3/(0.74*l0-12)]  -  129  dB.re  10“12 

for  s  =  3,  Lw  *  10  log  [  2. 8/(0. 74*10-12)]  -  126  dB,re  10"12 

Since  =  3(220  rad/sec)  =  660  rad./sec,  the  frequencie 

associated  with  the  foregoing  "pure  tone"  sound  levels  are 
a),  =  660,  u,  =  2(660),  (jj,  =  3(660)  rad/sec  or  fx  *  105,  f*  = 
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FIGURE  A1  FORCE  VARIATION  AT  TYPICAL  LOCATION  IN  PROPELLER 
PLANE,  b  =  PROPELLER  WIDTH  IN  PLANE,  AT  RADIUS  r 


FIGURE  A2  COORDINATES 
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FIGURE  A3  THICKNESS  FUNCTION  S(r1t$,)  IN  NON-ROTATING 
COORDINATES 


OV-1  PROPELLER 


FIGURE  A5  NARROW-BAND  SOUND  PRESSURE  IN  PLANE  OF  OV-1  PROPELLER 


figure  A6  narrow-band  sound 


(a)  MEASURED 

Cb)  DIPOLE  CALCULATION 

(c)  MONOPOLE  CALCULATION 

Id)  DIPOLE  PLUS  MONOPOLE  CALCULATION 


FIGURE  A 7  NARROW-BAND  SOUND  PRESSURE  OUT  OF  PLANE  OF  OV-1  PROPELLER 
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APPENDIX  B 

VORTEX  NOISE  OF  PROPELLERS 


VORTEX  NOISE  LEVEL 

For  low-speed  quiet  propellers,  the  broadband  noise  associ¬ 
ated  with  vortex  shedding  becomes  an  important  contributor  to  the 
acoustic  power  generated  (Ref.  1).  E.Y.  Yudin  (Ref.  2)  originally 
developed  a  theory  for  the  noise  generated  by  rotating  rods. 

H.H.  Hubbard  (Ref.  1)  extended  Yudin's  work  to  the  case  <  r  rotat¬ 
ing  airfoils,  postulating  that  propeller-generated  vortex  noise 
is  proportional  to  the  first  power  of  the  blade  area  and  the  sixth 
power  of  the  flow  speed,  Schlegel  et  al .  (Ref.  3)  have  developed 
an  empirical  prediction  scheme  for  the  vortex  noise  obtained  at 
fifteen  degrees  behind  the  propeller  plane. 

Schlegel 's  prediction  for  the  overall  vortex-noise  sound 
pressure  level  Lp  may  be  expressed  as 

Lp(9=15o)  -  "62  -  10  log  Ab  +  60  log  U0>7 

-  20  log  (R/3)  +  20  log  (Cl/0.4)  (1) 

where  Ab  is  the  area  of  the  propeller  blade  (in  ft2)  projected 
onto  the  propeller  plane;  U0>7  is  the  air  velocity  relative  to 
the  airfoil  at  the  seven  tenths  radius  (in  ft/sec),  R  is  the  dis¬ 
tance  to  the  propeller  hub  (in  ft),  and  Cr  is  the  blade  lift  co¬ 
efficient.  This  prediction  was  based  on  data  obtained  at 
Ct  =  0 .  1) ,  and  implies  that  there  is  net  vortex  noise  for  Cr  -  0; 
this  is  clearly  not  the  case  -  rotating  rods  provide  zero  lift 
but  do  make  vortex  noise. 

D.  Ross  (Ref.  4)  developed  a  scheme  for  dealing  with  vortex 
noise  involving  an  "acoustic  conversion  efficiency"  n ,  which  re¬ 
lates  the  acoustic  power  w  to  the  mechanical  drag  power  (i,e,,  to 
the  blade  section  drag  power)  Wm.  For  cylindrical  bodies  of  all 
shapes,  Ross  gives  the  conversion  efficiency  as 

r>  =  W/Wm  =  9.0  x  ,  (2) 

whore  M  is  the  Mach  number.  For  a  two-dimensional  airfoil,  the 
mechanical  power  ner  unit  length  is  given  by 


i\( i 


v  1 — «rr~ T*" -T 


(3) 


"m  ‘  I  pbU’°D 

where  b  is  the  foil's  chord  length,  U  the  local  flow  speed  and 
Cp  the  drag  coefficient.  The  acoustic  power  radiated  by  a  pro¬ 
peller  blade  of  radius  a  thus  may  be  written  as 

„  0, 

w  =  £JLI0_1  |  J  U6(r)b  (r)CD(r)dr  (4) 

°3  0 

where  c  =  sound  speed 

U(r)  =  flow  speed  at  radius  r 
b(r)  =  chord  length  at  radius  r 
CD(r)  *  drag  coefficient  at  radius  r. 

For  a  propeller  with  B  blades  and  a  constant  section  drag 
coefficient  CD  =  0.01,  operating  in  air  at  sea  level,  one  finds 
the  acoustic  power  to  be  given  by 


W(ft  lb/sec)  =  8  x  io~17B 


fo 


U 6 ( r ) £ ( r ) dr 


(5) 


for  U  in  ft/sec  and  b ,  r,  a  in  f t .  If  the  chord  length  is  con¬ 
stant,  the  blade  area  A^  =  Bab  may  be  Introduced,  and  the  pre¬ 
vious  expression  may  be  rewritten  as 


! .  0 

W  =  8  x  io" 17  J  U6 (r/a)d(r/a)  .  (6) 

0 

For  a  propeller  with  a  small  pitch  angle,  U(r)  is  very  nearly 
proportional  to  r,  and 


r  ut(^) i°  (-r)1  dtr) a 

where  U  is  the  blade  speed  at  the  seven  tenths 

0.7 


(7) 
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radius . 


By  substituting  Eq.  (7)  into  (6)  and  taking  the  logarithms,  one 
obtains 

10  log  W  =  -l60  +  10  log  Afa  +  60  log  U07  ,  (8) 

which,  when  converted  to  watts  and  referred  to  10-12  watts,  yields 
the  acoustic  power  level  as 

Lw(dB,  re  10~12  watts)  =  -40  +  10  log  Ab  +  60  log  U0?  (9) 

where  Ab  is  in  units  of  ft2  and  U07  in  ft/sec. 

If  the  vortex  noise  were  generated  by  a  family  of  dipoles, 
all  aligned  in  the  plane  of  the  propeller,  then  the  sound  pres¬ 
sure  would  vary  as  sin20,  where  0  is  the  angle  which  a  line  from 
the  observation  point  to  the  propeller  hub  makes  with  the  propel¬ 
ler  plane.  One  would  then  write  the  overall  sound  pressure  level 
at  the  observation  point  as 


L  (dB,  re  0.0002  microbar)  =  Lw(re  10~U  watts)  -  8 

+  20  log  sin0  -  20  log  ^ 

=  -48  +  10  log  Ab  +  60  log  7 

+  20  log  sin0  -  20  log  ^  .  (10) 

The  value  Cp  =  0.01,  for  which  the  foregoing  expressions  have 
been  derived,  is  estimated  to  be  typical  for  a  standard  propeller 
section  at  Cg  =  0.4,  A  comparison  on  this  basis  between  Eq. 

(10)  for  0  =  15°  and  Schlegel's  prediction  at  Cl  =  0.4  [Eq.  (1)] 
show;-,  that  Eq.  (10)  predicts  a  level  that  is  2  dB  greater  than  that 
given  by  Eq .  (1);  thus ,  the  agreement  between  the  theoretically 

and  experimentally  derived  prediction  schemes  is  seen  to  be  good. 

SPECTRUM  SHAPE 

Bchlegel  (Ref.  3)  gives  an  empirical  spectrum  shape  for  the 
vortex  noise.  The  peak  frequency  fp^of  this  spectrum  depends  on 


the  Strouhal  number  at  the  seven  tenths  radios ,  based  on  the  pro¬ 
jected  blade  thickness 

t 

dp  »  d  cosa  +  b  sina  (11) 


and  is  given  by 


0.28 


U 


(12) 


where  d  is  the  actual  blade  thickness,  b  is  the  chord  and  a  the 
blade's  angle  of  attack.  The  shape  of  the  spectrum  about  this 
peak  frequency  is  given  in  the  Table  below. 


OCTAVE  BAND  SPECTRUM  OF  PROPELLER  VORTEX  NOISE 


fOB/fpk 

1 

2 

1 

2 

8 

16 

loa~lob 

8 

4 

8 

9 

13 

1A 

fQB  =  center  frequency  of  octave  band 

fpk  =  peak  frequency  of  vortex  noise  spectrum 

Lqa  -  overall  level 

L„n  =  level  in  octave  band 

DIRECTIVITY 

If  all  of  the  dipole  sources  on  the  propeller  blades  were 
aligned  parallel  to  the  shaft  axis,  the  vortex  noise  field  would 
be  that  of  a  dipole  along  the  propeller  axis.  However,  in  fact 
the  dipoles  on  the  propeller  blades  are  not  aligned  from  blade  to 
blade,  due  to  the  pitch;  this  causes  a  smearing  of  the  directivity 
pattern.  An  aligned  dipole  field  would  radiate  nothing  in  the 
plane  of  the  propeller,  but  a  highly  pitched  propeller  can  gen¬ 
erate  considerable  noise  in  its  plane  of  rotation. 

At  an  angle  6  out  of  the  plane  of  rotation,  the  angle  be¬ 
tween  the  line  from  the  observer  to  the  propeller  hub  and  the 
normal  to  the  axis  of  a  dipole  on  a  blade  with  pitch  angle  8  is 


a  =  8  cosip-6 


(13) 


wnere  ip  is  the  angle  the  propeller  blade  makes  with  a  plane 
passing  through  the  shaft  axis  and  the  observation  point  (see 
Pig.  Bl).  The  sound  pressure  level  at  the  observation  point  (In 
the  farfield)  varies  as  sin2a;  the  average  level  for  a  rotating 
pitched  propeller  blade  may  be  written  as 


Lp(0  ,8 )  =  L  (n/2,0) 


1_ 

2tt 


•2  it 


sin2  (8  cosip-0)dip 


(14) 


in  terms  of  the  on-axis  sound  pressure  level  Lp(7r/2,0)  of  an  un¬ 
pitched  propeller  (at  the  same  distance  from  the  hub). 

By  expanding  the  integrand,  one  may  rewrite  the  integral  of 
Eq.  (l'i)  as 


f2a 

J  sin2 (-9  +  8 
0 


cosip)dip  =(sin2 


cos2 [3 


cosip]  dip 


2  tt 

+  (cosz0)  J  sinz[8  cosip]dtp 
0 

-  i  (sin  20)  J  sin[28  cosip]  dip  .  (15) 
0 

The  first  integral  on  the  right-hand  side  then  can  be  expressed 
in  terms  of  the  zero  order  Bessel  function,  J#: 


cos  2 [6 


costp]  dip 


J  =■  [1  +  cos ( 28 
0  ^ 


cosip)  ] dtp 


-  IT  [1  +  J  0  (  20  )  ]  , 


(16) 


from  which  it  follows  that  the  second  integral  obeys 


-2  TV 

J  sin2[B  cosip] dip  =  tt  [  1-J  ( 28  )  ]  .  (17) 

0 


The  last  integral  vanishes  because  sin[2B  costp]  is  an  odd  function 
of  ip. 


0  0 


By  use  of  the  foregoing  results  one  may  now  rewrite  Eq.  (1*1) 
as 


L  (6,P)  =  Lp(TT/2,0)  |  [1-J0(26)  cos  26] 


(18) 


The  modification  of  the  directivity  pattern  due  to  pitch  is 
particularly  significant  for  propellers  designed  to  operate  at 
very  low  rotation  speeds,  because  propellers  with  low  tangential 
to  axial  velocity  ratios  must  have  high  pitch.  The  flattened  di¬ 
rectivity  pattern  associated  with  high  pitch  results  in  a  lessen¬ 
ing  of  the  rise  and  fall  of  the  vortex  noise  recorded  on  the 
ground  as  a  propeller  driven  craft  flies  over  at  constant  alti¬ 
tude.  For  the  zero  pitch  case,  the  vortex  noise  level  is  great¬ 
est  when  the  angle  with  the  propeller  plane  is  45°  (peaking  once 
as  the  aircraft  approaches  and  once  as  it  flies  away).  For  a 
slow  speed  propeller  with  a  45°  pitch  angle,  the  sound  level  at 
45°  out  of  the  propeller  plane  is  about  the  same  as  for  the  un¬ 
pitched  propeller,  but  continues  to  increase  slowly  as  the  aircraft 
approaches  until  a  peak  is  reached  (only  a  few  dB  higher  than  that, 
at  45°)  at  26°  out  of  the  propeller  plane;  the  sound  level  then 
remains  constant,  within  0.5  dB,  until  the  26°  peak  is  passed  as 
the  aircraft  leaves . 

COMPARISON  OF  PREDICTIONS  WITH  EXPERIMENTAL  DATA 

By  combining  Eq .  (10)  with  (18)  and  including  a  drag  coeffi¬ 
cient  proportionality  relation,  one  obtains 


Lp  =  -48  +  10  log  Ab  +  60  log  tJ#<7  +  10  log 

+  10  log  {  \  [  1-J 0 ( 2  8)  cos  26]  |  -  20  log  .  (19) 

In  order  to  test  this  nrediction  relation,  a  search  of  the 
literature  was  made  for  well-documented  propeller  noise  data. 

Enough  must  be  given  so  that  one  can  evaluate  all  of  the  param¬ 
eters  which  appear  in  the  foregoing  relation.  In  addition,  in 
order  to  check  a  vortex  noise  prediction  against  measured  data, 
one  needs  a  narrow-band  spectrum,  so  that  one  can  separate  the 
broad-band  vortex  noise  from  rotational  and  engine  harmonic  noise 
(unless  the  vortex  noise  center  frequency  is  in  a  region  where 
vortex  noise  clearly  dominates  the  one-third  or  octave  band  levels). 
Of  noise  measurements  made  on  subsonic  propellers  and  helicopter 
rotors,  few  have  been  made  with  an  eye  toward  determining  the 
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vortex  noise.  Consequently ,  most  reported  measurements  do  not 
give  sufficient  data;  only  three  have  been  found  which  give  de¬ 
tail  to  permit  checking  the  vortex  noise  prediction  scheme:  one 
set  of  measurements  is  for  a  helicopter  rotor  (Ref.  5),  one  for 
a  hovercraft  pusher  propeller  (Ref.  6),  and  one  for  the  low-speed' 
propeller  of  a  "quiet"  aircraft. 

In  an  investigation  of  helicopter  rotor  noise  by  Stuckey  and 
Goddard  (Ref.  5),  a  3-bladed  NACA  0012  section  helicopter  rotor 
was  spun  on  a  whirl  tower  at  various  speeds  and  pitch  settings. 

The  rotor  thrust  and  sound  pressure  levels  were  r-'-’orded.  Tests 
were  conducted  at  rotational  speeds  between  180  and  ^50  revolu¬ 
tions  per  minute.  The  vortex  noise  levels  were  obtained  by  tak¬ 
ing  overall  levels  of  the  noise  between  100  and  2500  Hz,  which 
had  been  corrected  for  ground  reflection.  This  procedure  is 
justified,  because  this  frequency  range  extends  over  the  region 
of  major  vortex  noise  and  is  above  the  range  of  frequencies  where 
blade  rate  harmonics  contribute  significantly  to  the  noise  level. 
Figure  B2  shows  the  sound  level  measured  at  7-5°  behind  the  rota¬ 
tion  plane  and  151  ft  away  from  the  rotor  hub  (and  corrected  by 
Stuckey  and  Goddard  for  ground  reflection)  as  a  function  of  the 

mean  lift  coefficient  Cr  , 

Lm’ 


I,m  ■  L/ 


*b  2 


U2 ( r )dr 


*  L/k 2 


U2 ( 0 . 58  a) 


(20) 


which  relations  are  based  on  the  assumption  of  a  constant  lift 
coefficient  over  the  constant  chord  blade.  The  rotor  pitch  angle 
is  also  noted  in  the  figure.  The  dotted  curve  of  the  figure 
represents  the  data,  modified  by  use  of  Eq .  (18)  to  correspond  to 
a  zero  pitch  propeller.  For  the  rotor  blade  area  Ap  =  114.8  ft2 
and  for  U0<7  =  451  ft/sec  at  220  rpm,  the  zero-pitch  sound  pres¬ 
sure  level  at  the  measurement  location  is 


Lp(7.5°,  151  ft)  =  80  +  10  log  (21) 

Figure B3  shows  the  reported  drag  characteristics  of  the  NACA 
0012  section  (Ref.  7),  together  with  the  section  drag  coefficient 
which  is  required  to  make  the  Lp  calculated  from  Eq.  (21)  match 
tht:  zero-pitch  level  of  Fig.B2.  Since  the  mean  Reynolds  number 
for  drag  calculations  (at  0.58  a)  is  3.6  *  10s,  and  since  the  rotor 
is  likely  to  have  a  standard  roughness  finish,  this  calculated 
drag  curve  seems  to  be  not  an  unreasonable  one. 


In  a  series  of  noise  measurements  on  a  low  performance  liqht 
aircraft,  vortex  noise  data  was  obtained  for  two  propeller  condi¬ 
tions  and  at  two  different  angles  to  the  propeller  rotation  plane. 
The  tests  included  noise  measurements  in  the  propeller  plane  and 
on  the  rotation  axis  during  ground  run-up  tests,  and  an  additional 
measurement  as  the  aircraft  flew  past  a  microphone.  (Microphones 
were  mounted  over  woodchip-filled  boxes,  which  should  have  elimi¬ 
nated  ground  reflections  when  the  aircraft  was  at  a  relatively 
high  angle,  as  during  the  flyby.  The  woodchip  boxes  were  too 
small  to  affect  the  ground  reflection  path  from  the  aircraft  when 
tied  down.) 

For  a  highly  pitched  propeller,  the  axial  component  of  U0  7 
is  appreciable  and  cannot  be  neglected.  The  axial  velocity  for 
the  ground  run-un  test  was  estimated  from  the  given  shaft  horse¬ 
power  and  the  reasonable  efficiency  of  55%',  the  induced  axial 
velocity  was  computed  from  energy  considerations.  In  ground  run- 
un,  the  blades  of  this  propeller  are  under-pitched  at  the  tips, 
leading  to  an  increase  In  the  mean  lift  coefficient  by  an  esti¬ 
mated  15%.  Figure  B*t  shows  a  comparison  of  the  measured  vortex 
noise  with  that  predicted  by  use  of  Eq .  (19)  and  Schlegel's  spec¬ 
trum  share.  (The  measured  noise  for  the  tie-down  test  has  been 
corrected  for  3  dB  ground  reflection.)  The  prediction  for  the 
ground  run-up  is  probably  too  low  because  the  prone 1 lor  blade 
section  drag  was  underestimated  for  this  heavily  loaded  condition. 


Harrow-hand  sound  pm-. -sure  level  data  have  been  reported  for 
toe  oushor  propeller  of  a  hovercraft  development  I.td.  IID1  hover¬ 
craft  (rfef.  6) .  'i'h.-*  measurements  were  conducted  with  the  craft 

ti.-d  down  oil  a  concrete  runway.  Measurements  were  made  at  an 
angle  of  15°  behind  propeller  plane,  at  a  distance  of  100  ft  from 
the  hub ,  with,  the  microp.oone  5-5  ft  above  the  ground  surface.  The 
propeller  is  a  two-bladed  one,  with  a  7  ft  diameter,  a  total  blade 
ar°a  of  2.9  ft2,  a  blade  pitch  angle  at  0.75a  of  8°  1 0 T ,  and  a 
tfikness  at  0.7a  of  0.445  in.  Huns  were  made  at  1500  rpm  (23 
3 ■  s .  and  2*10  lb  thrust)  and  at  2500  rpm  (112  SUP  and  670  lb  thrust). 

The  blade  section  lift  coefficient,  based  on  the  foregoing 
information  and  on  the  assumption  of  constant  lift  coeffici  ‘fit 
over  the  length  of  the  1-lade  is  C,  *  0.6l  at  both  rotatit  speeds. 

.  i nee  the  propeller  pitch  angle  3  =  8°  10'  is  small  compared  with 
tno  out-of-plane  angle  6  =  15°,  the  directivity  function  can  be 
evaluated  as 


(22) 


10{ log  |  [1-J0(23)  cos  29 ] }  =  20  log  (sin9)  =  -11.8  . 

Because  of  the  hard  surface  over  which  the  microphone  was 
located,  the  predicted  sound  pressure  level  should  be  increased 
by  3  dB,  to  account  for  ground  reflection.  With  this  correction 
and  an  estimate  of  Cp  =  0.016  for  the  drag  coefficient  at  Cp  ~ 
0.61,  the  predicted  sound  pressure  levels  are. 


L  (dB,  re  0.002  microbar)  *  ! 

|  95, 

at 

1500 

rpm 

P  1 

1  82, 

at 

2500 

rpm 

For  an  assumed  4°  blade  angle  of  attack  of  a  1%  cambered 
section  operating  at  CL  =  0.61  (Ref.  7),  the  projected  blade 
thickness  is  found  from  Eq .  (11)  to  be  0.073  ft.  With  this  value 
Eq .  (11)  gives  the  center  frequencies  for  the  vortex  noise  to  be 
1*120  and  2360  Hz  for  1500  and  2500  rpm,  respectively.  Figure  B5 
shows  a  comparison  of  the  prediction  with  the  measured  data;  the 
agreement  is  seen  to  be  acceptable. 

ESTIMATION  SCHEME 

To  determine  the  spectrum  of  the  (broad-band)  noise  power 
generated  by  a  given  proneller,  one  may  proceed  as  follows; 

1.  Find  the  overall  acoustic  power  level  from 


Lw(dB,re  10“12  watts)  =  -  +  10  log  Afe(ft2)  +  60  log  U v ( ft/sec ) .  (9a) 

Here  A.  is  the  total  blade  area  and  U0  -  *  0.7aH  . 

2.  Find  the  projected  blade  thickness  from 


d  =  d  cos  a  +  b  sina  (11) 

c 

where  d  is  the  blade  thickness  and  b  the  blade's  chord 
length  (at  the  0.7a  radial  location),  and  a  its  angle  of 
attack.  Then  calculate  the  spectrum  peak  freauency  from 

fpk  -  0-28  U„/dp  (12) 
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3.  Use  the  table  below  to  develop  the-  octave  band  spectrum  of 
power. 

OCTAVE  BAND  SPECTRUM  OF  PROPELLER  VORTEX  NOISE 


f 


4  1  2  4  8  16 


OB  ok  2 


LOA“LOB  8  4  8  9  13  14 


-  center  frequency  of  octave  band 
=  peak  frequency  of  vortex  noise  spectrum 
=  overall  level 
=  level  in  octave  band 


Illustrative  Calculation 

Consider  the  three-b laded  propeller  described  at  the  end  of 
Apnendix  A,  again  rotating  at  2100  rpm.  The  total  blade  area  is 
Ab  =  4  ft2,  the  blade  thickness  is  d  =  0.2  ft,  the  blade  chord 
length  is  0.8  ft,  and  the  propeller's  angle  of  attack  is  15°. 

Then 

UM  =  0 . 7aft  -  0 . 7(  4ft )  ( 220  rad/sec)  =  615  ft/aec 

Lw  =  -  40  +  10  log  Ab  +  60  log  U  0(7 

=  -  4o  + 10  log  (4.0)  +  60  log(6l5)  =  -  40  +  10  (  0.6)  +  60(2.8) 
=  134  dB,re  10"14watts; 

d  =  d  cosa  +  b  sina  =  0.2  cos  15°  +0.8  sin  15°  =  0.40  ft 
n 

fpk  =  0.28  Uc>7/dr?  =  0.28(615  )/0 . 40  =  430  H z; 
and  from  the  table, 

fQB(Hz)  215  430  860  1720  3440  6900 

Lw(OB)(dB»re  10~12v;atts>  126  130  126  125  1?1  120 
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FIGURE  B1  DEFINITIONS  OF  ANGLES  FOR  DIRECTIVITY  ANALYS 


OVERALL  SOUND  PRESSURE  LEVEL 
(dB  RE  0.0002  MICROBAR) 
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FIGURE  B2  VORTEX  NOISE  AS  A  FUNCTION  OF  LIFT  COEFFICIENT  AND 
ANGLE  OF  ATTACK,  FOR  HELICOPTER  ROTOR  OF  REFERENCE 
5,  {AT  7.5°  BEHIND  ROTATION  PLANE) 


FIGURE  83 
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APPENDIX  C 


NOISE  OF  PISTON  ENGINES 


NOISE  SOURCES 

The  noise  of  piston  engines ,  like  that  of  any  air-breathing 
internal  combustion  engine, may  be  considered  in  terms  of  three 
components:  (1)  intake,  (2)  exhaust,  and  (3)  casing  noise.  In 
un-muffled  engines,  exhaust  noise  typically  predominates.  Casing 
noise  -  defined  as  the  noise  of  an  engine  with  "completely" 
muffled*  exhaust  and  intake  -  is  important  only  if  the  exhaust 
(and  possibly  the  intake)  indeed  is  well  muffled. 

The  exhaust  systems  of  Internal  combustion  engines  essen¬ 
tially  release  puffs  of  gas  (possibly  at  considerable  pressure) 
to  the  atmosphere  as  the  exhaust  valves  open  -  and  these  pres¬ 
sure  pulses  constitute  the  primary  components  of  exhaust  noise. 
Clearly,  exhaust  ducting  (such  as  manifolds,  pipes,  by-pass 
valves)  can  have  very  significant  effects  on  the  intensity  and 
shape  of  these  pressure  pulses  -  and  thus  on  the  exhaust  noise  - 
in  addition  to  adding  noise  (e.g.,  of  a  rushing  or  whistling 
nature)  due  to  the  flow  in  this  ducting.  (Ref.  1) 

The  foregoing  remarks  concerning  exhaust  noise  also  apply 
to  intake  noise,  except  that  the  intake  ingests  puffs  of  air, 
rather  than  emitting  them.  Acoustically,  ingested  puffs  have 
the  same  effect  as  emitted  ones;  however,  the  exhaust  pressure 
pulses  have  greater  pressure  differences  associated  with  them 
than  do  intake  pressure  pulses.  For  this  reason,  and  also  be¬ 
cause  high-oressure  pulses  tend  to  change  into  shock-waves  as 
they  propagate  along  a  duct,  unmuffled  exhausts  usually  are  much 
more  noisy  than  unmuffled  intakes. 

Casing  noise  may  be  ascribed  to  (1)  the  nearly  explosion¬ 
like  pressure  pulses  obtained  from  ignition  of  the  fuel,  (2) 
internal  mechanisms,  such  as  "piston  slap"  (lateral  impacts  of 
pistons  against  the  cylinder  walls)  and  impacts  associated  with 
valves  and  valve-lifters,  (3)  accessories,  such  as  fuel  pumps 
and  fuel  injectors,  and  ( <0  power  transmission  components,  such 
as  gears  and  bearings.  (Ref.  2) 


*mhe  exhaust  and  intake  of  an  engine  nay  he  considered  as  com¬ 
pletely  muffled  if  the  addition  of  further  attenuation  for  these 
sources  does  not  reduce  the  observed  noise  of  the  engines. 


DATA 


Availability  of  Data  on  Spark-Ignition  and  Diesel  Engines 

The  available  data  on  the  noise  of  spark-ignition  (SI)  en¬ 
gines  is  extremely  limited  and  so  inadequately  documented  that 
it  is  inadequate  to  serve  as  the  basis  for  development  of  a  noise 
prediction  scheme.  On  the  other  hand,  there  does  exist  a  useful 
body  of  information  on  the  noise  of  diesel  engines  (Pefs.  3,  *0 . 
Therefore,  it  appears  sensible  to  use  this  diesel  engine  data  as 
a  basis  for  estimation,  although  aircraft  are  more  likely  to  use 
SI  engines  than  diesel  engines. 

Similarity  of  Spark- Igni ti on  and  Diesel  Engine  Noise 

One  may  Indeed  argue  that  the  noise  of  SI  engines  should  be 
very  much  like  that  of  a  similar  diesel  engine.  Clearly,  the 
intake  process  does  not  depend  on  whether  the  fuel-air  mixture 
is  ignited  by  a  spark  or  by  compression  later  in  the  cycle,  so 
that  one  would  expect  the  intake  noise  of  a  given  engine  to  be 
the  same,  whether  that  engine  uses  spark  or  diesel  ignition. 
However,  the  intakes  of  SI  engines  may  be  designed  somewhat  dif¬ 
ferently  than  those  of  diesels  (because  SI  engines  may  have  car- 
burators,  for  example),  and  thus  may  sound  somewhat  different. 

The  same  sort  of  remarks  as  were  made  above  concerning 
intake  noise  also  apply  to  exhaust  noise.  Again,  the  exhaust 
puff  does  not  depend  on  the  ignition  process  -  but  rather  on 
the  exhaust  valve  timing,  on  the  cylinder  pressure  at  the  in¬ 
stant  the  valve  opens,  and  on  the  cylinder  pressure-time  charac¬ 
teristics  during  the  time  the  valve  is  open.  Although  "indicator 
diagrams"  (plots  of  cylinder  pressure  versus  piston  displacement 
or  crank  angle)  for  diesels  differ  considerably  from  those  of 
SI  engines  for  the  compression  and  combustion  narts  of  the  cycle, 
the  portions  of  these  diagrams  corresponding  to  the  exhaust  pro¬ 
cess  tend  to  be  very  similar.  (Ref,  1) 

On  the  other  hand,  one  might  expect  the  casing  noise  of  SI 
and  diesel  engines  to  differ  more  significantly.  Diesel  engines 
typically  are  designed  for  higher  compression  ratios  than  SI  en¬ 
gines,  and  the  (internal)  combustion  pressure  pulses  in  diesels 
tend  to  be  more  sharply  explosive  than  those  in  SI  engines.  The 
cylinder  pressure  spectra  of  diesels  thus  generally  are  higher 
than  those  of  SI  engines,  and  tend  to  have  much  more  high-frequency 
content  (Refs.  5-8).  However,  diesels  usually  are  built  much  more 
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massively  (to  contain  the  higher  pressures),  so  that  less  combus¬ 
tion  noise  tends  to  get  through  the  structure  (Refs.  9,10).  The 
severity  of  piston-slap  impacts  depends  on  the  cylinder  pressure 
to  some  extent,  so  that  these  impacts  are  more  severe  in  diesels 
(Refs.  4,11).  Again,  the  more  massive  structures  of  diesels 
tend  to  reduce  the  noise  radiated  as  a  result  of  these  Impacts. 
Piston-slap  is  likely  to  make  no  significant  contribution  to  the 
casing  noise  of  any  relatively  small  high-speed  engines  with 
tight-cylinder/piston  clearances,  which  one  might  use  in  aircraft 
(Refs.  4,11). 

Comparatively  little  is  known  about  the  accessory  noise  of 
the  various  types  of  engines  (Ref.  12).  There  appears  to  be  no 
reason  that  would  lead  one  to  believe  that  the  noise  of  a  power 
transmission  component  should  be  much  different  from  that  of  a 
similar  component  in  a  different  engine.  (Also,  these  compo¬ 
nents  usually  do  not  contribute  very  significantly  to  the  total 
casing  noise.)  Because  of  these  considerations,  and  in  the  light 
of  the  foregoing  discussion,  it  appears  logical  to  base  a  noise 
prediction  scheme  for  SI  engines  on  diesel  engine  data.  Such  a 
scheme  will  at  least  provide  a  first  approximation,  which  should 
prove  useful  until  adequate  SI  engine  data  become  available. 

Available  Data 

An  extensive  field  measurement  and  literature  review  program 
(Ref.  3)  has  yielded  noise  data  on  a  large  variety  of  diesel  and 
natural-gas  engines  in  electrical  power  plants  and  similar  sta¬ 
tionary  installations.  More  than  50  engines  (about  half  with 
4-stroke  and  half  with  2-stroke  cycles),  with  power  ratings  be¬ 
tween  12.5  and  5150  kw  (9.5  to  3850  Hp),  were  studied  at  actual 
power  outputs  between  0  and  5150  kw  and  at  speeds  between  225 
and  1800  rpm. 

The  "casing  noise"  of  these  engines  was  deduced  from  mea¬ 
surements  in  engine  rooms,  with  the  engine  intakes  and  exhausts 
either  outside  the  rooms  or  well-muffled.  Measured  acoustic 
characteristics  of  the  rooms  were  used  to  interpret  the  engine 
noise  data  in  terms  of  acoustic  power  spectra  of  the  sources. 

For  some  of  these  engines,  the  air  intake  or  engine  exhaust 
noise  was  measured  out  of  doors,  at  some  convenient  distance 
from  the  intake  or  exhaust  opening.  Where  no  muffler  was  present, 
these  measurements  were  directly  interpreted  in  terms  of  acoustic 
power ;  where  exhaust  mufflers  were  present,  the  measured  data 
were  corrected  by  means  of  estimated  attenuation  characteristics 
of  the  muffler  and  exhaust  piping  to  yield  the  approximate  acoustic 
power  of  the  unmuffled  exhaust. 
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The  data  of  Ref.  3  were  re-analyzed  in  order  to  provide  an 
improved  method  for  predicting  the  engine  noise  spectra  of  inter¬ 
est  here.  This  re-analysis  involved  discarding  all  data  on  gas 
engines  (because  these  are  inherently  quieter  than  diesels),  all 
data  on  engines  with  nonstandard  structures  (e.g.,  with  thermal 
insulation),  and  all  data  on  engines  with  Roots  blowers  (because 
these  blowers  add  to  the  noise  and  are  not  likely  to  be  used  in 
aircraft);  in  addition,  it  involved  re-examining  the  data  in 
terms  of  a  reduced  frequency  based  on  the  "firing  frequency". 

The  results  of  this  analysis  are  shown  in  the  attached  figures. 

DATA  ANALYSIS;  ESTIMATION  SCHEME 

Exhaust  Noise 

Figure  Cl  summarizes  octave-band  exhaust  noise  data  obtained 
from  measurements  on  nine  diesel  engines  with  widely  differing 
operating  characteristics.  In  addition  to  points  corresponding 
to  individual  measurements,  this  figure(as  well  as  all  the  other 
figures  of  this  section)  also  shows  "average"  and  "average  plus 
standard  deviation"  spectra,  which  were  obtained  by  carrying 
out  the  appropriate  arithmetic  operations  on  the  data  correspond¬ 
ing  to  the  points  indicated  at  each  frequency.*  Because  the 
acoustic  power  of  engine  noise  sources  has  been  found  to  be  pro¬ 
portional  to  the  engine  power,  the  sound  power  levels  in  Fig.  1 
and  in  the  rest  of  the  figures  of  this  section  have  been  reduced 
with  respect  to  the  engine  power. 

Because  the  exhaust  valve  of  an  engine  cylinder  opens  once 
per  cycle  (i.e.,  once  per  power  stroke  or  once  per  firing),  the 
exhaust  of  an  engine  emits  one  puff  per  firing  per  cylinder. 

One  would  thus  expect  the  exhaust  noise  spectrum  to  have  a  ma.ior 
peak  at  the  firing  frequency  fp,  which  for  an  N-cylinder  engine 
may  be  found  from 

f _ ( Hz )  =  engine  speed  (rev/sec)  *  n  •  N, 

r 


*The  averaging,  etc.,  computations  were  carried  out  not  on  the 
decibel  values  which  correspond  to  the  data  points,  but  on  the 
mean-square  pressures  they  represent.  The  "energy"  averages  ob¬ 
tained  in  this  way  are  likely  to  be  most  meaningful  for  estima¬ 
tion  purposes.  (Ref.  13) 


where  n  is  the  number  of  power  strokes  per  revolution.  For  a 
four-stroke  cycle  engine,  n  =  1/2;  for  a  two-stroke  eyrie  engine, 
n  =  1. 

The  data  of  Fig.  Cl  are  shown  in  Fig.  C2,  replotted  against 
the  ratio  of  octave  band  center  frequency  to  firing  frequency. 
Better  collapse  of  the  data  is  evident.  Also  shown  in  Fig.  C2 
is  a  curve  which  was  suggested  by  Franken  and  Beranek  (Ref.  14) 
for  estimating  the  average  exhaust  noise  levels  one  may  expect 
from  unmuffled  piston  engines;  this  curve  may  be  seen  to  encom¬ 
pass  most  of  the  peaks  of  the  average  data  curve,  but  to  over¬ 
estimate  the  high-frequency  noise  considerably.  An  improved  es¬ 
timation  curve,  which  follows  all  of  the  data  peaks  of  Fig.  C2, 
is  presented  in  Fig.  C7. 

Intake  Noise 

Intake  noise  data  obtained  from  measurements  on  six  engines 
with  unmuffled  intakes  are  summarized  in  Fig.  C3.  The  same  data 
are  shown  in  Fig.  as  a  function  of  firing  frequency.  Unlike 
for  exhaust  noise,  this  replotting  does  not  improve  the  collapse 
of  the  data  materially.  However,  since  it  does  bring  out  the 
existance  of  a  physically  meaningful  spectrum  peak  near  the 
firing  frequency,  it  appears  advantageous  to  base  an  estimation 
scheme  on  Fig.  C4,  rather  than  on  Fig.  C3.  The  corresponding 
curve  of  Fig.  C7  represents  a  smoothed  envelope  of  the  peaks 
of  the  "average"  spectrum  of  Fig.  C4. 

Casinn  Noise 

Figure  C5  shows  octave-band  data  on  casing  noise  for  27 
diesel  engines  with  widely  varying  engine  parameters,  as  indi¬ 
cated  in  the  figure  heading.  The  sound  power  levels  indicated 
again  are  normalized  with  respect  to  engine  output  horsepower 
(Hp)  and  include  A  and  C  correction  terms,  which  were  found 
useful  in  Ref.  3  for  improving  the  estimates  of  the  overall 
levels.  An  "average  curve",  which  Is  expected  to  be  useful  for 
general  prediction  purposes.,  (calculated  on  the  basis  of  energy 
averaging,  after  discarding  the  most  extreme  points)  is  also 
given  in  the  figure. 

figure  c6  differs  from  Fig.  C5  only  in  that  the  data  are 
plotted  against  the  ratio  of  octave  band  center  frequency  to 
engine  firing  frequency.  A  somewhat  better  collapse  of  the  data 
is  evident,  so  that  use  of  Fig.  C6  for  estimation  purposes  appears 
preferable.  The  casing  noise  estimation  curve  given  in  Fig.C7 
is  a  smoothed  envelope  of  the  peaks  of  the  "average"  curve  of 
Fig.  C6, 


ESTIMATION  SCHEME 


Pig.  C7  summarizes  an  estimation  scheme  for  all  three 
noise  components;  how  the  various  estimation  curves  were  de¬ 
veloped  was  indicated  above.  As  evident  from  Pig.  C7,  the 
noise  of  unmuffled  exhausts  always  predominates.  With  a  well- 
muffled  exhaust,  casing  noise  becomes  most  important.  Silencing 
of  the  intake  cannot  reduce  the  total  noise  of  an  engine  unless 
the  exhaust  and  casing  noise  components  have  been  reduced  dras¬ 
tically  . 

In  order  to  predict  the  noise  of  a  given  engine,  cue  may 
proceed  as  follows: 

1.  First  calculate  the  firing  frequency  from 

fp(.Iz)  =  engine  speed  (rev/sec)  •  n  •  N 

where  N  is  the  number  of  cylinders  and 
for  four-stroke  cycle  engines 

1  for  two-stroke  cycle  engines. 

Then  determine  what  actual  frequencies  f  correspond  to  the 
reduced  frequencies  f/'fp  of  Fig.  C7. 

2.  Find  10  log (Hp ) ,  A,  C,  using  the  relations  given  in  Fig.  C7. 
For  each  of  the  three  noise  components  represented  in  the 
figure,  use  the  values  of  Lw  -  10  log(Hp)  +  A  +  C  read  from 
the  figure  to  determine  the  corresponding  octave  band  spectra 

Note  that  one  may  obtain  these  spectra  simply  by  relabeling 
the  frequency  and  Ly,  scales,  and  by  shifting  the  casing 
noise  curve  downward  by  the  amount  A  +  C  with  respect  to 
the.  other  two  curves. 


Illustrative  Calculation 

Consider  an  unmuffled  "V-fi"  6-cylinder  four-stroke-cycle 
engine  producing  90  Ho  at  2100  rpm.  Here 

fp  =  rps  •  n  •  M  =  (p)  (6)  =  100  Hz. 

Thus  f/f-i  =  1  corresponds  to  f  =  105  Hz,  f/fp  -  \  to  f  =  52  Hz, 
f/fp  =  rto  1»20  Hz,  etc. 


Since  2100  rpm  >  1500  rprn,  one  finds  A  =  0,  and  for  a  "V" 
engine,  C  =  1;  and 

10  log (Hp )  =  10  log(90 )  =  19.5  dB 

Therefore  one  obtains  the  engine  noise  spectra  by  relabeling  the 
vertical  coordinate  scale  of  Pig.  C?  as  simply  I^(dB,r>-.  10" 1 2v;atts ) , 
increasing  all  numbers  indicated  along  this  scale  by  19.5  dB,  and 
shifting  the  casing  noise  curve'  downward  by  1  dB. 

One  thus  finds,  for  example,  that  the  peak  exhaust  noise 
power  level  is  131.5  dB  (at  52  and  105  Hz),  that  the  neak  intake 
noise  power  level  is  91.5  dB(re  10~12watts)  at  105  Hz,  and  that 
the  highest  casing  noi.se  level  of  108.5  dB  occurs  between  8^0 
and  1680  Hz. 
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APPENDIX  D 


NOISE  OF  ROTATING  COMBUSTION  ENGINES 


DATA 

Data  on  the  noise  of  Curtiss-Wright  Corp.  rotating  combustion 
(Wankel)  engines  are  presented  in  Ref.  1.  These  data,  adjusted 
for  engine  speed  and  horsepower  output  (so  as  to  make  the  data 
collapse  into  a  form  best  suited  for  prediction  purposes),  are 
summarized  in  the  figures  which  follow  this  text. 

As  in  virtually  all  engines,  intake  and  exhaust  noise  tends 
to  dominate  over  noise  from  all  other  sources.  This  noise  from 
"all  other  sources",  which  one  observes  if  the  intake  and  exhaust 
are  effectively  muffled,  is  called  casing  noise  and  is  labeled  as 
such  in  the  figures.  (Also  see  corresponding  discussion  under 
the  heading  of  "Noise  of  Piston  Engines".) 

AI R-COOLED  AND  WATER-COOLED  ENGINES 

The  data  points  indicated  in  the  figures  correspond  to 
measurements  carried  out  on  two  RC-2-60  water-cooled  engines*, 
onerating  at  between  2500  and  5000  rpm  while  producing  between 
30  and  103  Hp.  The  two  engines  to  which  these  data  points  apply 
differ  only  in  their  intake  port  configurations;  one  has  side- 
ports,  the  other,  peripheral  ports.  The  peripheral  port  con¬ 
figuration  tends  to  result  in  noise  levels  which  are  slightly 
higher  (by  at  most  0.5  db  for  casing  noise,  1.5  db  for  exhaust 
noise,  and  2.5  db  for  intake  noise)  than  those  for  side-ports. 

Although  detailed  data  on  the  noise  of  air-cooled  engines 
arc  not  given  in  Ref.  1,  that  document  does  present  comparisons 
between  the  noise  levels  of  an  air-cooled  RC-2-90  engine  and 
those  of  water-cooled  PC-2-60  engines.  The  curves  for  air¬ 
cooled  engines  in  the  attached  figures  were  obtained  from  the 
differences  betvieen  typical  air-cooled  and  water-cooled  engine 
data . 


*furt ' ss-Wright  uses  the  first  number  after  the  RC  to  designate 
the  number  of  rotors,  and  the  second  to  give  the  displacement 
(in  cubic  inches)  per  rotor. 


DISCUSSION 


The  data  available  were  measured  on  only  three  different 
types  of  engines,  all  having  two  rotors,  with  about  the  same 
displacement.  The  effects  of  adding  more  rotors  are  not  known, 
but  one  might  expect  the  noise  spectrum  to  shift  with  the  firing 
frequency . 

No  data  on  the  effects  of  engine  size  (rotor  displacement) 
are  available  at  all.  The  number  of  Wankel  engine  types  in  opera¬ 
tion  is  very  limited;  some  small  engines  of  lawh-mower  size  are 
in  production,  as  are  some  engines  for  small  automobiles,  but  no 
engines  rated  significantly  above  100  Hp  appear  to  be  available 
as  yet. 

ESTIMATION  SCHEME 

Clearly,  an  average  of  the  measured  spectra  should  serve 
for  prediction  purposes.  Probably  the  most  meaningful  average 
is  that  which  corresponds  to  the  mean-square  value  computed  for 
all  of  the  data  points  at  each  frequency.  Correspondingly,  the 
"average"  and  "average  plus  standard  deviation"  curves  indicated 
in  the  figure  have  been  obtained, not  by  performing  the  appropriate 
arithmetic  operations  on  the  decibel  numbers  represented  by  the 
various  data  points,  but  by  performing  those  operations  on  the 
corresponding  mean-square  values.  (See  Ref.  2.)  Because  of  the 
relatively  large  standard  deviations  calculated  for  most  of  the 
data,  "average  minus  standard  deviation"  values  generally  turn 
out  to  be  too  low  to  fit  on  the  plots  and  therefore  do  not  appear 
in  the  figures. 

Average  curves  for  all  three  types  of  noise  from  rotating 
combustion  engines  are  summarized  in  Fig.  DA. 

In  order  to  estimate  the  acoustic  power  produced  by  a  given 
rotary  combustion  engine,  calculate  15  log (rpm/1000 )  +  10  log(Hp) , 
add  this  value  to  the  numbers  given  on  the  vertical  scale  of  Fig. 

DA,  and  relabel  that  scale  simply  L^dBjre  10“ 1 2watts ) . 

Illustrative  Calculation 

For  an  air-cooled  rotating  combustion  engine  producing  90  Hp 
at  2100  rpm, 

15  log(rnm/1000 )  +  10  log(Hp)  =  15  log(2.1)  +  10  log(90)  *  2A.5  dB. 

Thus,  for  example,  the  110  dB  value  of  Fig.  DA,  at  which  the  exhaust- 
noise  oeak  occurs  (at  1000  Hz),  corresponds  to  13A.5  dB. 
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APPENDIX  E 


NOISE  FROM  TURBOSHAFT  ENGINES 


AVAILABLE  DATA 

Although  no  broadly  useful  collection  of  data  pertaining  to 
the  noise  of  aircraft  turboshaft  (turbo-prop)  engines  appears  to 
be  available,  such  a  collection  of  data  has  been  assembled  for 
stationary  gas  turbine  installations  (Ref.  1).  Accordingly,  it 
appears  useful  to  develop  a  noise  prediction  method  on  the  basis 
of  this  collected  stationary  engine  data,  until  more  directly 
relevant  data  become  available.  The  data  reported  in  Ref.  1 
were  obtained  largely  from  measurements  carried  out  on  operating 
installations  rated  between  2^0  and  1500  Hp  and  rotating  at  be¬ 
tween  1,200  and  ^3,030  rpm, 

NOISE  SOURCES 

The  noise  from  turbine  inlets  and  exhaust  tends  to  be  most 
prominent.  If  these  sources  are  well-muffled,  one  is  left  with 
the  "casing  noise"  of  the  turbine. 

As  air  flows  through  the  turbine  inlet  into  the  compressor, 
the  compressor  blades  "chop"  the  flow;  the  compressor  acts  like 
a  siren  and  produces  a  rather  pure-tone  noise  component.  This 
component  typically  results  in  a  peak  in  the  noise  at  the  "blade- 
passage  frequency"  of  the  rotating  blades  in  the  first  compressor 
stage.  The  blade-passage  frequency  fR(Hz)  is  given  by  the  pro¬ 
duct  of  the  number  of  blades  in  that  stage  and  the  compressor 
speed  in  revclutions/sec .  Another  pure-tone  component,  which 
often  results  in  a  minor  spectrum  peak,  usually  occurs  at  the 
shaft  rotation  frequency.  This  frequency  fs(Hz)  is  equal  to  the 
shaft  rotation  speed  in  revolutions/sec. 

The  flow  v/ithin  a  turbine  also  produces  noise  by  siren-action 
to  some  extent.  But  there  the  corresponding  pure-tone  components 
tend  to  be  masked  by  broader-band  noise  components  associated  with 
turbulence  and  turbulence-related  fluctuations  in  lift  and  drag 
of  the  turbine  blades.  Accordingly,  casing  noise  and  exhaust  noise 
tends  to  be  more  nearly  of  a  broad-band  character,  although  spec¬ 
trum  peaks  at  the  blade-passage  and  shaft-rotation  frequencies 
tend  to  be  present. 


COLLAPSE  OF  DATA 


Figure  El  summarizes  all  of  the  intake  noise  data  available 
from  Ref.  1,  in  terms  of  octave-hand  power  level,  reduced  with 
respect  to  rated  power,  plotted  against  actual  frequency.  (An 
attempt  at  referring  the  noise  to  the  actual  power  produced  by 
the  turbines  resulted  in  much  greater  scattering  of  the  data; 
turbines  make  considerable  noise  even  when  operating  at  zero 
power  output,  and  the  noise  of  a  turbine  often  increases  rela¬ 
tively  little  as  its  power  output  increases.  Although  other 
parameters,  such  as  intake  flow  rate,  possibly  may  affect  the 
noise,  insufficient  data  is  available  for  the  assessment  of 
their  importance.) 

In  addition  to  all  of  the  data  points.  Fig.  El  also  shows 
three  sample  spectra  which  illustrate  the  typical  double-peaked 
character  of  inlet  noise,  with  peaks  at  the  blade-passage  and 
shaft-rotation  frequencies. 

Figures  E2  and  E3  show  the  same  data  as  Fig.  El,  but  plotted 
against  the  ratio  of  frequency  to  blade-passage  and  to  shaft- 
rotation  frequency/,  respectively.  In  the  construction  of  Fig.E2, 
only  data  near  the  high-frequency  peak  of  spectra  like  those 
shown  in  Pig.  El  were  retained;  low-frequency  data  were  omitted. 

In  cases  where  not  enough  information  was  available  to  permit 
determination  of  the  blade-passage  frequency,  it  was  assumed  for 
purposes  of  constructing  Fig. E2  that  this  frequency  corresponds 
to  the  aforementioned  high-frequency  peak. 

Figure E3  similarly  shows  only  data  near  the  low-frequ'  ic 
peak  of  the  original  un-reduced  spectra.  Here,  however,  the 
shaft  rotation  frequency  was  known  in  all  cases. 

pigureEU  presents  exhaust  noise  data  In  terms  of  actual 
frequency,  and  Figs.E5  andE6  present  the  same  data,  but  in 
terms  of  ratios  of  frequency  to  blade-passage  and  shaft-rotation 
frequencies.  Figures E7 ,E8 ,  andE9are  similar  plots  for  casing 
noise . 

The  spread  of  the  inlet  noise  data  is  seen  to  be  consider¬ 
ably  less  in  the  reduced  frequency  plots  of  Figs .  E2  andE3  than 
in  the  actual  frequency  plot  of  Fig.  El.  Similarly,  the  exhaust 
noise  data  cluster  somewhat  better  in  Figs,E5  andE6  than  in  Fig. 
E^,and  the  casing  roise  data  collapse  better  in  Figs.  E8  and  E9 
than  in  Fig.  ET . 
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ESTIMATION  CURVES 


The  "average"  and  "average  plus  standard  deviation"  spectra 
shown  in  the  figures  have  been  calculated  by  performing  the  ap¬ 
propriate  arithmetic  operations  on  the  mean-square  values  which 
correspond  to  the  data  points  at  each  frequency.  Such  energy 
averages  (Ref.  2)  are  expected  to  be  most  meaningful  for  estima¬ 
tion  purposes.  (Note  that  no  "average  minus  standard  deviation" 
spectra  are  given  in  the  figures.  Because  of  the  relatively 
large  standard  deviations  found  here,  such  spectra  would  fall 
off  the  bottom  of  the  graphs.) 

The  various  "average"  spectra  are  most  applicable  for  gener¬ 
al  estimation  purposes,  whereas  the  "average  plus  standard  devia¬ 
tion"  spectra  are  useful  for  estimating  the  highest  noise  levels 
one  is  likely  to  obtain. 

Because  the  average  inlet  noise  spectrum  has  pronounced 
peaks  at  the  blade-passage  and  shaft-rotation  frequencies,  one 
may  base  an  estimation  scheme  on  Figures  E2  and  E3.  Such  a  scheme 
is  presented  in  Fig.  E10}  which  suggests  estimating  the  inlet 
noise  spectrum  of  any  turbine  engine  in  two  parts  -  a  low-frequency 
part  centered  at  the  shaft-rotation  frequency,  and  a  high-frequency 
part  centered  at  the  blade-passage  frequency*  -  and  interpolating 
between  these  two  parts  as  needed.  Comparison  of  Fig.  E10  ;« 'th 
Fig.  El  shov/s  that  the  scheme  of  Fig.  Eli)  preserves  the  essential 
character  of  the  unreduced-frequency  spectra. 

If  one  attempts  to  construct  a  prediction  scheme  for  exhaust 
noise  like  that  which  was  developed  for  inlet  noise,  than  one 
finds  that  either  the  two  partial  spectra  must  coalesce,  or  one 
cannot  reproduce  the  essential  character  of  the  average  spectrum 
of  Fig.  E4  .  Because  the  shaft-rotation  frequency  is  physically 
more  meaningful  for  exhaust  noise  than  is  the  blade-passage  fre¬ 
quency,  a  prediction  scheme  based  on  the  shaft-rotation  frequency 
is  suggested.  This  scheme  is  presented  in  Fig.  Ell  ,  which  is 
based  on  averaging  and  smoothing  of  the  spectra  of  Figs.  E5  and  E6 
(aligned  so  that  the  peak- frequencies  coincide).  Again,  Fig. 

Ell  may  be  seen  to  be  very  similar  to  the  average  spectrum  of  Fig. 
E4, which  is  plotted  in  terms  of  actual  ( non-reduced)  frequency. 


*The  partial  spectra  of  Fig.  ElO  were  obtained  by  smoothing  the 
average  spectra  of  Figs.  E2andE3and  by  omitting  those  portions 
of  these  spectra  which  are  so  far  from  the  peak-frequencies 
that  the  estimation  is  likely  to  be  unreliable. 
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Finally,  one  finds  that  the  two  reduced  casing  noise  spectra 
of  Figs .  E8  and  E9  cannot  be  made  to  coalesce  in  any  meaningful 
way  to  produce  an  estimation  which  resembles  the  average  spectrum 
of  Fig.  E7.  Also,  it  is  likely  that  casing  noise  may  be  dom..nated 
by  the  effects  of  auxiliary  components  and  casing  structural  para¬ 
meters,  rather  than  by  the  siren-like  blade-passage  effects.  Thus, 
the  average  spectrum  of  Fig.  E7,  which  spectrum  is  reproduced  in 
somewhat  smoothed  form  in  Fig.  E12,  is  likely  to  be  best  for  es¬ 
timation  purposes. 

ESTIMATION  SCHEME 

In  order  to  nredict  the  noise  produced  by  a  given  turboshaft 
engine,  one  should  deal  with  each  of  the  three  major  noise  compo¬ 
nents  (inlet,  exhaust,  casing)  separately. 

1.  Use  Fig.  E10  to  estimate  the  inlet  noise  power.  Calculate 
the  shaft  rotation  frequency  fg  from 

fs(Hz)  =  compressor  rotation  speed  (rev/sec). 

and  the  blade-passage  frequency  fg  from 

f  (Hz)  =  f o  •  number  of  blades  in  first 
compressor  rotor  stage 

Determine  the  actual  frequencies  that  correspond  to  the 
frequency  ratios  indicated  in  Fig.  E10  and  the  power  levels 
L-w  that  correspond  to  the  reduced  levels  shown. 

Plot  the  levels  corresponding  to  each  octave-band  center 
frequency . 

Re-plot  the  two  partial  snectra  obtained  by  this  process, 
and  connect  the  noint  corresponding  to  f/fs  “  ^  with  that 
for  f/fg  =  1/^  bv  a  straight  line  (corresponding  to  the 
dotted  line  of  Fig.  F,1D).  If  the  two  partial  spectra 
overlap,  use  a  smoothed  (upner  bound}  envelope  for  the  total 
spectrum. 

2.  Use  pig.  Ell  to  estimate  the  exhaust  noise,  by  replacing 
the  frequency  ratios  in  that  figure  with  the  corresponding 
frequencies  and  the  reduced  power  level  by  the  actual  power 
level , 

3.  Use  Fig.  E12  to  estimate  the  casing  noise,  by  replacing  the 
reduced  power  levels  indicated  by  the  corresponding  actual 
power  level. 


8? 


Illustrative  Calculation 

Consider  a  turboshaft  engine  rated  at  90  Hp,  with  a  com¬ 
pressor  which  has  22  blades  in  its  first  (rotating)  stage,  and 
which  rotates  at  27,000  rpm. 

Since  here  10  log(Hp)  =  19.5  dB  and  fs  =  27,000/60  =  450  Hz, 
fg  *  (450)  (22)  =  9,900  Hz,  one  finds  the  following  power  levels 
corresponding  to  Pig.  E10 : 


f/fs 

1/4 

1/2 

1 

2 

4 

f/fB 

1/4 

1/2 

12  4 

f  (Hz) 

112 

225 

450 

900 

1,800 

2,500 

5,000 

10,000  20,000  40,000 

Lw(dB) 

95 

103 

105 

102.5 

101 

103 

111 

118  116  113 

By  shifting  the  computed  values  to  the  nearest  standard  octave- 
band  frequency  and  connecting  those  two  partial  spectra,  one 
obtains  the  curve  shown  in  Fig.  E13. 

Using  the  same  values  of  fg  and  10  log(Hp)  in  conjunction 
with  Figs.  Ell  and  E12,  one  similarly  obtains  the  exhaust  and 
casing  noise  curves  also  shown  in  Fig,  E13. 
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APPENDIX  F 

DISCRETE-TONE  NOISE  PRODUCED  BY  TURBOJET  FANS 


INTRODUCTION 

Fans  and  compressors  of  turbojet  engines  typically  produce 
high-pitched  whining  sounds,  which  consist  primarily  of  one  or 
more  pure-tone  components.  Axial-flow  fans  resemble  multi-bladed 
propellers,  and  the  basic  mechanisms  responsible  for  fan  noise 
are  in  part  the  same  as  those  which  cause  propeller  noise.  Also, 
in  both  cases,  the  pure-tone  components  are  associated  with  pres¬ 
sure  fields  that  spin  about  the  axis  of  rotation. 

However,  propeller  noise  theory  is  inapplicable  to  fans  (and 
its  use  can  lead  to  gross  underestimates  of  fan  noise).  The  duct 
in  which  fans  are  housed  affect  the  aerodynamics,  as  well  as  the 
propagation  of  sound,  and  the  dominant  sources  of  fan  noise  are 
associated  with  the  interactions  of  successive  rows  of  (rotor  and 
stator)  blades. 

SOURCE  MECHANISMS  (Ref.  1-3) 

Two  sources  of  pure-tone  noise  are  of  primary  importance  for 
subsonic  rotors  -  i.e.,  for  rotors  which  move  subsonically  rela¬ 
tive  to  the  air  flow:  (1)  The  fluctuating  lift  generated  on  a 
downstream  blade  row  as  it  interacts  with  the  wakes  formed  by  tne 
upstream  blade  row  produces  an  effective  acoustic  force  field  in 
the  plane  of  the  blades.  (2)  A  siren-like  volume  source  effect 
is  produced  as  successive  blade  rows  present  a  periodic  obstruc¬ 
tion  to  the  flow;  the  obstruction  is  minimum  when  the  wakes  are 
in  line,  and  maximum  when  they  are  alternately  spaced. 

The  mechanisms  present  in  subsonic  fans  also  are  present  in 
supersonic  fans,  but  in  supersonic  fans  there  occur  also  other 
mechanism,  which  may  be-  of  greater  importance.  These  mechanism 
are  associated  with  the  rotor  shock-wave  field,  which  radiates 
directly,  and  which  also  produces  lift  fluctuations  on  the  stator 
blades.  Because  of  the  greater  noise  from  supersonic  rotors, 
such  rotors  should  be  avoided  in  "quiet  aircraft"  applications, 
and  will  not  be  discussed  further  here. 
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PREDICTION 


Much  theoretical  work  has  been  done  on  the  various  noise 
mechanisms  (e.g.,  Refs.  2-H) .  However,  most  of  the  relations 
that  have  been  developed  are  not  very  useful  for  prediction  of 
noise  because  they  require  information  on  various  aerodynamic 
details  which  is  generally  not  available  until  a  prototype  has 
been  built.  And  by  then,  of  course,  one  may  obtain  much  better 
data  on  the  noise  simply  by  measuring  it. 

Perhaps  the  best  means  for  predicting  the  level  of  the  fun¬ 
damental  pure-tone  component  generated  by  a  fan  consists  of  a 
relation  derived  in  Ref.  3,  which  gives  the  sound  pressure  level 
Ln ( 200 » )  Produced  at  200  ft  to  the  side  of  fan  engines,  at  the 

blade-passage  frequency,  as 

Lp(2001  )(dB>re  0.0002ubar)  =50  log  utl  (ft/sec)  +  20  log  D(in)-75 

(1) 

in  terms  of  the  tip  speed  U^ip  of  the  fan  blade  and  the  fan  dia¬ 
meter  D. 

If  one  assumes  the  directivity  to  be  uniform  -  an  assumption 
that  is  not  strictly  valid,  but  that  is  justified  for  the  simpli¬ 
fied  scheme  for  estimation  of  aural  detection  suggested  in  this 
report  -  then  one  may  easily  determine  the  sound  power  level  that 
corresponds  to  En.  (1).  Noting  that  lLin  =  SID  =  Dirfir/30,  where 
Q  represents  the  fan  rotational  speed  in  radianr./sec  and  SI  that 
same  speed  in  rpm,  and  by  changing  the  dimensions  of  D  to  rft , 
one  may  obtain 

Lw(dR,re  10“12  watts)  =  50  log  Or(rpm)  +  70  log  D(ft)  -  56  (2) 

In  Ref.  3  there  is  also  presented  a  comparison  of  Eq,  (1) 
with  experimental  data  obtained  on  6  different  fans.  For  tip 
sneeds  up  to  about  700  ft /sec,  Eq.  (1)  is  very  nearly  an  upner 
bound  to  the  data;  the  data  in  this  region  fall  largely  within 
+  1  dB  and  -  7  dB  of  the  values  given  by  Eq.  (1).  For  transonic 
tip  sneeds,  i.e.  between  700  and  1500  ft/sec,  the  data  are  more 
widely  scattered  and  fall  essentially  between  +  5  dB  and  -  10  dB 
of  Eq.  (1). 

Although  it  is  well  known  that  fan  noise  in  general  also 
contains  several  significant  harmonics  of  the  fundamental  blade- 
nassage  frequency,  no  adequate  means  are  available  for  predicting 
the  sound  levels  associated  with  these  harmonics.  Experimental 


data  (Ref,  2)  also  Indicate  wide  fluctuations  in  the  relative 
levels  of  the  first  few  harmonics,  even  in  successive  runs  of  the 
same  engine.  On  the  basis  of  the  very  limited  amount  of  data 
available,  it  may  suffice  for  aural  detectability  estimation  pur¬ 
poses  to  consider  only  the  first  three  harmonics  of  the  fundamen¬ 
tal,  and  to  take  the  levels  of  successive  harmonics  to  decrease 
by  3  dB. 

ESTIMATION  SCHEME 

To  estimate  the  discrete-tone  sound  levels  produced  by  a 
turboject  fan,  proceed  as  follows: 

1.  Calculate  the  blade-passage  frequency 

fB  =  Bfl 

where  B  represents  the  number  of  rotor  blades  in  the  fan's 
last  stage,  and  fi  the  fan's  rotational  speed  (in  radians/sec). 
Then  use  Eq.  (2)  to  determine  the  acoustic  power  Lw  produced 
at  that  frequency.  A 

2.  Then  take  -  3  dB  at  2  fg 

1^  -  6  dB  at  3  fB 

and  Lw  -  9  dB  at  4  fg 

as  the  remaining  significant  pure-tone  components  for  de¬ 
tectability  estimation  purposes. 

3.  Repeat  steps  1  and  2  for  the  fan’s  next-to-last  stage,  but 
reduce  all  levels  by  3  dB  (Ref.  1), 

Illustrative  Calculation 

Consider  a  1,65  ft  diameter  fan  with  22  blades  in  its  first 
and  29  blades  in  its  second  (and  last)  stage  rotor,  operating  at 
5,800  rpm  (  =  610  rad/sec). 

For  the  last  stage,  then,  fR  =  Bfl  =  29(610)  ~  17,500  Hz  and 
L  =  50  log  nr  +  70  log  D  -  56 

=  50  log(23, 300 )  +  70  log (1 . 65 )  -  56  =  171  dB.re  10"12  watts. 

Thus  one  obtains  171  dB  at  17.5  kHz,  168  dB  at  35  kHz,  165  dB  at 
52.5  kHz,  and  162  dB  at  70  kHz. 


Then  for  the  first  stage,  one  finds  that  at  fg  »  22(610)  * 
13.5  kHz,  L  =  171  ~  3  ■  168  dB,  and  at  27  kHz  one  obtains  165  dB 
etc.  w 

Note,  however,  that  frequencies  above  30  kHz  are  of  no  con¬ 
cern  in  relation  to  aural  detection.  Also  see  illustrative  calcU' 
lation  at  end  of  Appendix  H. 
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APPENDIX  G 


NOISE  OF  GAS  JETS 


The  primary  source  of  noise  from  any  well-muffled  and  acous¬ 
tically  enclosed  aircraft  propulsion  system  (e.g.,  from  a  turbo¬ 
fan  engine  "buried"  in  the  fuselage  of  an  aircraft  and  provided 
with  appropriate  inlet  and  exhaust  ducts)  is  due  to  the  thrust- 
producing  Jet  efflux  from  the  system. 

Much  work  has  been  done  on  the  noise  of  jets,  and  well- 
validated  prediction  methods  are  available,  particularly  for 
jets  with  relatively  small  velocities  (Refs,  1-3),  which  are 
likely  to  be  of  greatest  interest  in  relation  to  small  quiet 
aircraft . 

ACOUSTIC  EFFICIENCY  AND  MECHANICAL  POWER 

Perhaps  the  best  and  simplest  estimation  method  may  be  ob¬ 
tained  on  the  basis  of  the  summary  study  undertaken  by  Heckl 
(Ref.  4).  It  has  often  been  observed  that  similar  rockets  or 
jets  have  very  nearly  the  same  "acoustic  efficiency"  o  which 
is  defined  as 


O  -  W/V/  ,  (1) 

ac  m  ' 

where  W  represents  the  total  (overall)  acoustic  power  radiated 
and  Wn  the  mechanical  power  of  the  flow.  This  mechanical  power 
is  given  by 


Wm=PtUj  (2) 

whereF^.  denotes  the  thrust  and  Uj  the  jet  velocity.  The  thrust 
of  air-breathing  propulsion  systems  obeys*  (Ref.  5) 

?t  ■  -  V  +  ,nfD|  +  Vpj  -  V  (3> 


*Note  that  this  relation  applies  for  any  type  of  system  that 
produces  thrust  by  accelerating  air;  it  thus  applies  equally 
well  for  propeller  and  jetTengine  systems. 


where  ma  represents  the  mass-flow  rate  of  air  through  the  propul¬ 
sion  system,  which  flow  enters  the  system  with  a  relative  velocity 
U0  (typically  the  vehicle  air-speed)  and  leaves  with  velocity  U. 
Similarly,  mf  represents  the  mass-flow  rate  of  fuel,  which  enters 
the  propulsion  system  with  essentially  zero  velocity  and  exits 
with  velocity  Uj .  In  addition,  Aj  represents  the  area  of  the  jet, 

Pj  the  pressure  in  the  jet,  and  P0  the  ambient  (intake  air)  pres¬ 
sure.  In  most  propulsion  systems,  ma/mf  >  50,  and  P0  *  Pi,  so 
that 

\  ’  ">a(uj  -  V  "  "auj  "  »°jA  *  (1|> 

where  p  represents  the  density  of  air  in  the  Jet  and  the  usual 
condition  U.  >>  U0  has  been  assumed.  Thus,  the  mechanical  power 
may  be  approximated  by 

W  -  pA.U.3  .  (5) 

m  J  J 

Figure  G1  summarizes  the  observed  dependences  of  the  acoustic 
efficiency  on  Jet  velocity,  taken  from  Ref.  b.  The  band  Indi¬ 
cated  for  low  velocities  represents  the  scatter  of  data  in  this 
region;  higher-turbulence  jet  flows  here  have  significantly 
higher  acoustic  efficiencies.  Also  indicated  In  the  figure  are 
expressions  which  give  the  dependences  of  oac  on  Mach  number 
corresponding  to  the  various  straight  lines  shown  in  the  figure. 
The  dotted  line  within  the  low-frequency  band  represents  an  ap¬ 
proximate  energy  average  between  the  upper  and  lower  bounds  of 
that  band,  and  may  be  used  for  general  estimation  purposes. 

SPECTRUM  SHAPE  AND  DIRECTIVITY 

The  Jet  noise  spectrum  may  be  estimated  by  use  of  Fig.  G2 , 
which  shows  how  the  octave-band  levels  differ  from  the  overall 
level.  The  spectrum  has  the  typical  haystack  shape,  with  its 
peak  at  the  dimensionless  frequency  (Strouhal  number) 
fU4/D=»0.2,  where  f  represents  the  actual  frequency  and  D  the 
jet  (nozzle)  diameter. 

A  curve  for  the  estimation  of  directivity  effects  is  given 
In  Fig.  03.  It  shows  that  the  highest  sound  levels  occur  at  about 
30°  from  the  Jet  efflux  axis,  and  that  the  sound  levels  ahead  of 
the  Jet  (at  180°  from  the  efflux  axis)  are  nearly  10  dB  lower 
than  the  spatial  average  value. 


ESTIMATION  SCHEME 


In  "Tier  to  estimate  the  noise  produced  by  a  given  Jet  at  a 
given  location,  one  may  proceed  as  follows: 

1,  Use  the  mechanical  power  Wm  anfi  the  diameter  D  of  the  Jet 
to  calculate  the  Jet  velocity  Uj®  of  an  equivalent  Jet  at 
room  temperature  and  zero  altitude  from 

-  tVV*  ’  <5a> 

taking  p  *  0.07  lb/ft3  . 

2.  Use  Fig.  G1  to  find  oac,  and  from  oac  find  the  overall  acous¬ 
tic  power  of  the  equivalent  Jet  from 


w  =  °ac^m 


( la) 


Then  calculate  the  corresponding  power  level  from 
Lw(dB,re  10“12watts)  =  10  log(W/10_1 2watts ) 


3.  Add  to  this  overall  power  level  the  correction  AL^  for 
atmospheric  pressure  and  temperature 

SLw  =  10  log  (Pat/Patj0)  ♦  "0  log  <T,/T0)  -  35(Ta/Tc) 
where  P  .  =  actual  atmospheric  pressure 

3  v 

p  .  =  standard  atmospheric  pressure 

at  ,o 

T.  =  absolute  temperature  of  Jet 

T  =  absolute  temperature  of  ambient  air 
a 

T  =  absolute  reference  (room)  temperature 
o 


4. 

5. 


Apply  the  directivity  correction  from  Fig.  G3,  if  directivity 
information  is  of  interest.)* 


Calculate  the  Jet  diameter  from 


or 


D2  =  i4Ft/TrpU2 


D2  =  4W  /ttpU? 

m  J 


(4a) 

(5b) 


and  use  Fig.  G2  to  obtain  the  octave  band  noise  spectrum. 


•Directivity  effects  are  not  used  in  the  estimation  method  for  aural 
detection  suggested  in  this  renort.  Directivity  estimation  is  in¬ 
cluded  here  only  for  the  sake  of  completeness. 
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Illustrative  Calculation 


Consider  a  propulsion  system  which  generates  a  thrust  of 
155  lb  by  producing  an  air  jet  that  e  anates  from  a  1.65  ft  dia¬ 
meter  duct  at  300  ft/sec  and  250°F.  The  mechanical  power  of  this 
Jet  is 


Wm  =  FtUJ  =  ^155  lb)(30°  ft/sec)  =  46,500  ft  lb/sec 

and  its  area  is  Aj  =»  ttD2/4  *  2.44  ft2.  Hence 

Wm  (4,65x104ft  lb -/see)  /32.21b  ft\ 

U»  -  -P- - - - - £ -  - U  8.80  x  106 

Je  pAJ  (0. 071b m/ ft3) (2. 44ft2)  '  lbfsec2  / 


U,  «*  206  ft/sec. 

J  ° 

From  Fig.  Gl,  o  a  2  *  10“ 7 ; 

3.C 

W  *  o  W  *  ( 2*10“ 7 )  ( 4 . 65x10 “ft  lb/sec)(1.36  watt  sec/ft  lb) 
sc  m 


=  1,26  x  10”2  watts 


Lw  =  10  log(  1, 26x10“ 2 /10“  1  2  )  =  10  log(l. 26x10*  °)  ■  101dB,rel0“12watts 


If  the  ambient  air  is  at  T  =  4o°F  =  500°R  and  Pat,o  *  0.9atm, 
whereas  the  reference  air  temperature  is  * 

T0  =  70°F  =  530°R,  the  correction  required  is 


ALW  =  10  log 


/  p  \  /T  \  /  T  \ 

'__at__j  +  40  log^j-  35  log^J 


at  ,o 


-0.4  +5.0  +  0.8 


*  5.5  dB 

Therefore  the  overall  sound  power  level  is  Lw  +  ALW  * 
106.5  dB,re  10“l2watts, 


105 


i-  4«.  .  ;;■  »-  «  -  ,^w  •*?  ■*',  *•  '  .'* f^®3? 


Since  U./D  =  (300  ft/aec)/(1.65  ft)  -  l80/sect  fD/Uj  =  1 
corresponds  J  to  f  *  180  Hz.  One  may  thus  assign  frequencies  to 
the  Strouhal  numbers  shown  in  Fig.  G2.  Some  octave-band  power 
levels  found  from  the  previously  determined  overall  level  of 
106.5  dB,  as  obtained  from  Fig.  G2,then  are  as  follows: 


ft-hei?';.*  isjSEf3jfe&. 


' 


fOB(Hz) 


18 

36 

72 

180 

360 

720 

96 

100 

90 

9^ 

91 

87 

It 
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APPENDIX  H 

SOUND  ATTENUATION  IN  CIRCULAR  DUCTS 


INTRODUCTION 

The  attenuation  of  sound  propagating  in  ducts  is  a  function 
of  the  sound  source  and  of  the  geometrical  and  acoustical  para¬ 
meters  of  the  duct. 

The  source  may  be  described  in  terms  of  the  spatial  distri¬ 
bution  of  the  sound-producing  elements,  and  is  usually  expressed 
in  terms  of  wave  number  parameters  or  modes.  This  source  descrip¬ 
tion  also  involves  frequency-dependent  impedance  terms,  which  can 
be  related  to  the  various  individual  modes.  Although  flow  can 
affect  the  source  characteristics,  it  is  convenient  to  consider 
this  effect  separately. 

The  geometry  of  a  duct  can  be  described  in  terms  of  the  ra¬ 
tios  of  length  to  diameter  of  the  lined  duct,  the  ratio  of  lining 
thickness  to  duct  width  (or  diameter),  and  the  ratio  of  the  length 
of  the  lined  section  to  that  of  the  unlined  section  of  the  duct. 
One  then  needs  only  one  absolute  length  to  describe  the  geometry, 
and  one  may  compare  this  length  to  the  acoustic  wavelength. 

The  acoustical  properties  of  a  duct  are  commonly  described 
in  terms  of  the  propagation  constants  Tn  and  an  of  modes,  which 
constants  are  defined  by  the  axial  pressure  function 

-ikz(Tn-icn)  ikz(Tn-ion) 

P(z)  =  F  Pn+e  +  F  Pn-e  (1) 

n  n 

where  k  denotes  the  acoustic  wavenumber  (in  free  space).  The 
propagation  constants  are  related  to  characteristic  modal  impe¬ 
dances  Zn  by 


Z„  =  pc/(,n  -  lon>.  (2) 

In  designing  a  duct,  the  one  variable  over  which  one  has 
control,  in  addition  to  the  geometry,  is  the  wall  impedance  Z. 

For  a  rectangular  duct,  without  flow,  the  impedance  of  a  homo¬ 
geneous,  locally  reacting  wall  is  related  to  the  propagation  con¬ 
stants  as  (Refs.  1-4) 


(3) 


f  »,V, . 


_  (T^.iCT^jz  tan  (kH/l  -  (Tn-ion)*) 


where  H  denotes  the  width  of  the  rectangular  duct  (assumed  lined 
on  one  side).  For  a  circular  duct  of  diameter  D,  one  may  deter¬ 
mine  that 


if£  -  TTr-roTT1 


‘V1  V!] 


(4) 


where  J0,  denote  the  Bessel  functions  of  order  0  and  1.  (This 
result  is  derived  in  the  last  part  of  this  section.) 

Air  flow  through  the  duct  and/or  high-intensity  sound,  leads 
to  changes  in  the  effective  wall  Impedance  and  in  the  relation 
between  the  wall  impedances  and  the  axial  sound  propagation  (Ref. 
5).  The  changes  are  small,  however,  for  centerline  flow  speeds 
that  correspond  to  Mach  numbers  below  0.1  and  for  sound  pressure 
levels  below  1A0  dB  (re  0.0002  microbar).  Because  high  flow 
speeds  and  sound  pressure  levels  are  of  limited  interest  here, 
their  effects  are  neglected  in  the  subsequent  discussion. 

Sound  absorption  at  the  duct  walls  and  turbulent  air  flow 
cause  the  modes  to  interact.  Therefore,  the  concept  of  modes, 
v/hich  is  a  mathematical  device  introduced  to  describe  sound 
fields  in  rigid  enclosures  in  terms  of  orthogonal  functions,  is 
not  a  useful  one  here.  In  particular,  it  makes  no  sense  to  con¬ 
sider  modal  distributions  in  a  duct  v/hich  are  not  compatible  with 
any  source  configuration.  Because  of  these  problems  with  modes 
concept,  and  also  because  of  the  mathematical  difficulties  aris¬ 
ing  from  the  complex  transcendental  equations  which  describe  the 
relation  between  wall  impedances  and  modal  parameters  [Eqs.  (3) 
and  (A)],  the  analysis  of  sound  attenuation  presented  here  is 
based  on  approximations,  which  approach  the  idea  of  modes  only 
insofar  as  a  stabilized  cross-distribution  of  sound  pressure 
(with  minimum  attenuation  along  the  duct)  can  be  identified  with 
a  fundamental  mode, 

HIGH  FREQUENCY  ATTENUATION 

It  is  convenient  to  consider  two  limiting  frequency  regions; 
a  high-frequency  region  in  which  the  acoustic  wavelength  X  is 
small  compared  to  the  duct  diameter  D,  and  a  low-frequency  region 
in  which  X  >  D. 


For  high  frequencies,  the  sound  propagation  in  ducts  may  be 
analyzed  in  terms  of  rays.  Figure  HI  shows  a  cylindrical  duct  of 
length  L  and  diameter  D  with  incoherent  point  sources  assumed  lo¬ 
cated  at  one  end,  in  the  z  -  0  plane.  The  source  density  is  taken 
to  be  a  function  q(r)  of  the  radius  r,  but  not  of  the  angular  po¬ 
sition  4>. 


Directly  Emitted  Rays 

If  the  duct  is  lined  with  an  absorbing  material,  most  of  the 
sound  radiated  out  of  the  duct  is  due  to  rays  which  propagate 
through  the  duct  without  being  reflected  from  the  walls.  Such 
unreflected  rays  which  emanate  from  an  element  rdrd<j>  of  the  source 
region  are  confined  to  a  spherical  angle  ft0 .  If  no  duct  were  pre¬ 
sent,  the  same  source  element  would  radiate  into  the  spherical 
angle  2rr.  Therefore,  the  acoustic  intensity  due  to  rays  emitted 
from  the  element  rdrdtf  and  propagating  without  reflections  is 


dl  -  dl0 


(5) 


where  dl  denotes  the  intensity  obtained  with  a  very  short  or 
completely  rigid  duct.  The  intensity  emitted  from  tne  entire 
source  region  is  a  function  of  the  source  distribution  q(r)  and 
obeys 


D/2 

f  q(r) 


dr 


1  =  1, 


D/2  i 


U/  C  r  O.-. 

f  ( r )  /'sin-^-J 

J  o 

The  duct  geometry  affects  the  angle  eo ; 
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dr 


sin 
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-  u  + 


ft 


2  )-V2 


(6) 


(7) 


Considerable  reduction  of  the  sound  intensity  can  be  obtained 
only  if  0  is  small;  then 

e0  „  e„  „  D 


T  *  T  ~  2L  1 


(8) 


11 1 


sin 


for  any  source  distribution. 

Reflected  Rays;  Attenuation  of  Short  Ducts 

/ 

In  order  to  take  account  of  the  imperfect  absorptivity  of 
the  duct  lining,  one  must  consider  the  rays  that  are  reflected 
from  the  walls  before  they  leave  the  duct.  The  rays  which  suffer 
one  reflection  are  confined  to  the  spherical  angle  -  Q0 »  where 
=  2n  sin6./2  corresponds  to  rays  that  are  emitted  from  a  duct 
with  a  diameter  of  3D,  as  shown  in  Pig.  H2.  If  one  again  limits 
oneself  to  a  relatively  long  duct  for  which  Eq.  (8)  holds,  then 
one  finds  chat  the  contribution  from  rays  which  are  reflected  be¬ 
fore  they  leave  the  duct  is  given  by 

I  =  (10) 

where  a  represents  the  absorption  coefficient  of  the  duct  wall. 
From  Eqs,  (9)  and  (10),  one  may  then  determine  the  attenuation 
(in  decibels)  to  obey 

AL  =  -  10  log  £  (|-  -ot);  (11) 


this  relation  provides  one  with  an  estimate  of  the  performance  of 
short  ducts  lined  with  an  absorbing  material.  Because  the  absorp¬ 
tion  coefficient  a  usually  increases  with  increasing  angle  of  in¬ 
cidence  or,  equivalently,  because  rays  which  suffer  more  reflec¬ 
tions  are  attenuated  more  highly,  one  usually  may  neglect  the 
contributions  from  multiply  reflected  rays. 

Attenuation  of  Long  Ducts 

One  may  analyze  the  multiple  reflections  of  rays  in  long 
ducts  by  considering  the  mean  free  path  l  of  a  ray  between  suc¬ 
cessive"  reflections  at  the  walls.  The  intensity  in  the  duct  then 
decays  according  to 


=  I 
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(12) 


For  weakly  attenuated  rays,  which  propagate  almost  axially, 
the  mean  free  path  length  £  is  found  to  vary  directly  as  the  duct 
diameter  D  and  inversely  as  the  (small)  angle  e  between  the  wall 
and  the  initial  direction  of  propagation 


£  =  E 


(13) 


The  set  of  rays  which  results  in  a  pressure  minimum  at  the  wall 
suffers  the  least  attenuation.  For  this  set  o-f  rays,  the  angle 
e  is  determined  by  the  ratio  of  the  radial  to  the  axial  wavenum¬ 
bers  , 


k  k 

e  «  tan  e  =  r~  s  X 
kz  2  it 


(14) 


For  circular  ducts,  the  pressure  minimum  at  the  wall  is  de¬ 
termined  by  the  first  zero  of1  the  Bessel  function  J0(krr);  this 
zero  occurs  for  krr  *  2.4,  corresponding  to  which 


kr  =  4.8/D 


(15) 


Accounting  only  for  rays  which  emanate  from  the  source  nearly 
axially  (e  <<  ,L),  one  obtains  the  attenuation  AL  (in  decibels), 
as 


AL  =  4.35 


oL  _X_ 
D  2tt 


4.8  _  _  ,  LX 
D  3'3aZ  D  D 


(16) 


The  absorption  coefficient  a,  for  such  rays  usually  is  small,  and 
can  be  approximated  by 

<*. =  }  •  <17> 

where  Re{Z/pc}  is  the  real  part  of  the  wall  impedance  for  nor¬ 
mally  incident  sound,  divided  by  the  characteristic  impedance  pc 
of  the  gas, in  the  duct. 

Since  the  attenuation  in  a  longer  duct  turns  out  to  be  pro¬ 
portional  to  the  square  of  X/D  -  which  is  small  in  the  region  of 
high  frequencies  -  no  great  improvement  can  be  obtained  by  ex¬ 
tending  a  short  duct.  The  Initial  section,  for  which  the  attenu¬ 
ation  is  approximately  given  by  Eq.  (11),  contributes  almost  all 


of  the  attenuation.  If  one  desires  to  increase  the  attenuation 
substantially,  one  needs  to  curve  the  duct  or  add  obstructions 
to  interact  with  the  axial  rays. 


LOW-FREQUENCY  ATTENUATION 

Analysis  of  the  performance  of  lined  ducts  at  frequencies  at 
which  the  ducts  are  not  very  long  compared  to  the  wavelength  in¬ 
volves  three  problems.  They  are  associated  with  the  source  re¬ 
gion,  with  transitions  from  a  rigid  duct  to  a  lined  duct  (and 
vice  versa),  and  with  the  duct  termination. 

Source  Distribution 

As  previously  mentioned,  the  spatial  distribution  of  sources 
can  be  described  only  inadequately  by  modes.  At  low  frequencies, 
at  which  the  acoustic  wavelength  is  greater  than  the  duct  diam¬ 
eter,  the  spatial  distribution  becomes  stable  in  a  very  short 
distance  along  the  duct.  The  stabilized  radial  distribution  suf¬ 
fers  the  least  attenuation  with  distance;  one  may  therefore  ob¬ 
tain  a  conservativ.  estimate  of  the  attenuation  performance  of 
lined  ducts  by  cons  '.dering  only  the  stabilized  distribution. 

This  approach,  whi  <  .1  does  not  account  for  the  actual  source 
characteristics,  ii  employed  below. 

Duct  Terminations  and  Lining  Discontinuities 

The  exact  calculation  of  reflection  and  transmission  coeffi¬ 
cients  at  the  transition  from  a  rigid  duct  to  a  lined  duct  is  ex¬ 
tremely  difficult.  Heins  and  Feshbach  (Ref.  6)  applied  a  Wiener- 
Hopf  technique  to  develop  series  approximations.  The  first  term 
of  their  series  agrees  witn  results  obtained  from  the  simplest 
considerations,  using  the  basic  modal  impedances  Zp  =  pc  for  rigid 
ducts  and  Zc  =  pc/(x-lo)  for  lined  ducts.  Such  approximations 
are  appropriate  in  the  case  of  moderate  attenuation,  and  will  be 
employed  in  the  present  analysis.  Effects  of  duct  terminations 
have  been  studied  by  Levine  and  Schwinger  (Ref.  7),  who  have  cal¬ 
culated  the  radiation  from  an  unflanged  pipe.  One  finds  that  one 
may  obtain  an  adequate  approximation  by  considering  only  the  mass 
reactance  in  parallel  with  the  plane  wave  impedance  at  the  open 
end  of  a  rigid  duct  (Ref.  8).  Results  for  the  radiation  from 
lined  ducts  generally  show  somewhat  higher  values  than  obtained 
from  this  approximation;  this  difference  can  be  ascribed  in  part 
to  the  curvature  of  the  wavefront  in  the  duct. 

The  geometry  of  discontinuities  in  the  duct  lining  and  of 
unflanged  open  tailpipes  to  be  studied  is  shown  in  Fig.  H3,  which 
indicates  a  lined  duct  of  length  Ls ,  connected  to  a  tailpipe  of 
length  Lp  with  rigid  walls;  the  duct  and  tailpipe  are  taken  to 
have  the  same  diameter  D. 
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The  transmission  of  sound  through  individual  homogeneous 
parts  of  the  duct  may  conveniently  be  described  by  means  of 
transfer  matrices,  which  involve  the  propagation  constant  and 
the  characteristic  impedance  of  the  stabilized  radial  distribu¬ 
tion  (or  fundamental  mode).  The  effects  of  series-arrangements 
of  different  parts  may  then  be  determined  by  multiplication  of 
these  transfer  matrices.  Thus,  the  sound  pressure  and  velocity 
at  the  inlet  of  the  absorbing  duct  are  related  to  those  quanti¬ 
ties  at  the  open  end  of  the  rigid  tailpipe,  to  the  first  order 
of  approximation,  by 


sin[kLs(T-io)) 


Pi  \  /  cos[kLs(T-io)]  i  - -  \ 

^pc  /  \  i  (t-ic )  sin[kLs  (t-io)  ]  cos[kLs  <  t-i  o  )  ]  ' 


cos  kL 


R 


i  sin  kL 
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pcv, 


(18) 


At  the  open  end  of  an  unflanged  pipe  of  diameter  D,  the 
pressure  and  velocity  are  related  by 


v2pc 

~p7" 


1  + 


4 

fkD 


(19) 


One  may  then  calculate  the  attenuation  of  the  duct  in  terms  of 
the  ratio  of  the  pressure  at  the  inlet  of  the  absorbing  duct  to 
that  at  the  open  end  of  the  tailpipe,  as 


I 

AL  =  20  log 


dB 


(20) 


The  inlet  pressure  depends  on  the  source  conditions,  as  dis¬ 
cussed  below.  ' 

Source  Connected  to  Absorbing  Duct  via  Unlined  Duct 

Here  the  pressure  pt  results  from  the  complex  amplitudes  p1+ 
and  p  of  an  Incident  wave  and  of  the  wave  reflected  at  the 
wave-impedance  discontinuity  at  the  inlet  of  the  absoi’blng  duct; 
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(21) 


Pi  =  Pi+  +  P 


1  - 


For  the  plane  wave  in  the  rigid  inlet  duct,  there  applies  the 
relation 


vlPc  =  p1+  -  Pl_ 


(22) 


If  the  source  Is  not  specified  otherwise,  one  may  assume  that  its 
internal  impedance  is  equal  to  the  characteristic  impedance  pc. 
Then  the  pressure  p(i)  is  equal  to  the  pressure  p1+  of  the  Inci¬ 
dent  wave,  and  from  Eqs.  (18)  to  (22)  one  obtains  the  attenuation 
(in  decibels)  as 


AL  =  8.7koL  +  10  log  (l  +  — - — 
s  \  (kD)2 


(23) 


+  20  log 


1  + 


-2ikLs( t-1o ) 

- (  l-T  +  i0) 

4(f-io) 


e-21kLR 

1-T+iO  +■ — - (1  +  T-io) 

l+ikD/2 


The  first  term  represents  the  attenuation  of  a  homogeneous  lined 
duct.  The  second  term  gives  the  reflection  loss  at  the  end  of 
an  unflanged  pipe.  The  part  of  the  third  term  which  involves 

1-e  2ikLs(T  io)  accounts  for  waves  propagating  back  and  forth  in 
the  lined  duct,  which  result  from  reflections  at  the  discontinui¬ 
ties  in  the  duct  lining;  the  part  of  the  third  term  with  the  co¬ 
efficient  e  accounts  for  waves  reflected  betv/een  the  end  of 

the  lined  duct  and  the  open  end.  (Note  that  the  latter  reflected 
waves  disappear1  for  kD  ■*“,  which  corresponds  to  an  anechoic  termi¬ 
nation  of  the  rigid  tailpipe), 

For  strong  reflections  from  the  open  end  of  the  tailpipe 
(kD-*-0),  and  weak  interaction  between  the  two  ends  of  the  absorb¬ 
ing  duct,  the  third  term  of  Eq .  (23)  reduces  to 


20  log 


i  +  istua  fi .  T+l0  ♦  e"21kLR 

(T-iO ) L 


(  1  +  T-iC ) 


which  for  t  >  1  >>  a  becomes 
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Thus,  quarter  wavelength  resonances  in  the  tailpipe  (kLR  =  n/2) 
result  in  positive  values  and  half  wavelength  resonances  (kLR  =  u) 
in  negative  values  of  the  third  term  of  Eq.  (23)»  as  long  as  t  >1, 
which  is  the  case  below  the  first  resonance  in  the  duct  lining. 

For  very  short  absorbing  ducts  (kLs  ■*  0),  the  third  term  of 
Eq.  (23)  becomes 


20  log  1  + 


-(  1-T  +  iO  ) 
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-2ikLR 

x  +  io  +  ^ - (1+  -lo) 
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Sound  Source  Located  Directly  at  the  Inlet  of  the 
Absorbing  Duct 

For  a  low-impedance  source  (pressure  source),  the  pressure 
i(i ^  which  characterizes  the  attenuation  performance  of  the  duct 
s  equal  to  the  pressure  p! ,  Equations  (18)  and  (19)  here  yield 
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Within  the  approximation  valid  for  x  *  1 ,  the  third  term  here 
accounts  for  reflections  from  the  discontinuity  of  the  duct  lin¬ 
ing  and  from  the  end  of  the  tailpipe,  in  agreement  with  Eq .  (23). 

For  a  high-impedance  source  (velocity  source)  at  the  inlet 
of  the  absorbing  duct,  the  reference  pressure  P(j)  results  from 
the  velocity  Vj  multiplied  by  the  characteristic  Impedance 
pc/(r-io),  ana  one  obtains 
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This  result  differs  from  Eq.  (24)  only  in  the  algebraic  sign  re¬ 
lating  to  the  phases  of  reflected  waves  in  the  lined  duct. 

The  source  impedance  does  not  affect  the  attenuation  perfor¬ 
mance  of  the  duct,  if  the  reflected  wave  from  the  discontinuity 

of  the  duct  lining  can  be  neglected  (e“2kLsCr<<  1  ). 

PROPAGATION  CONSTANTS  FOR  CIRCULAR  DUCTS 

The  linearized  wave  equation  in  cylindrical  coordinates  r,z 
is 
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With  the  axial  dependence 
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P<3>  '  e"ikz(T-ia) t 


(27) 


one  obtains  Bessel’s  differential  equation  for  the  radial  depen¬ 
dence  , 


+  i  IS  +  k2 [i  -  ( T-io ) 2 ]p  =  0. 


9r2  r  9r 

The  solution  compatible  with  a  finite  axial  pressure  is 
p(r)  ~  J# [kr/1  -  ( x -1  a ) 2 ] , 
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where  J#  is  the  zero  order  Bessel  function.  Boundary  conditions 
at  a  homogeneous,  locally  reacting  wall  at  r  *  D/2 
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can  be  described  by  the  wall  impedance 
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Thus,  the  propagation  constants  t ,c  may  be  calculated  from 
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The  various  solutions  of  this  transcendental  equation  are  common¬ 
ly  attributed  to  various  modes. 

A  reasonable  approximation  for  the  fundamental  mode  (which 
has  the  smallest  attenuation  constant  e)  can  be  obtained  by  a 
finite-difference  approximation: 

If  (rn}  3  S[p(rn+^i)  '  p(rn^)J 


(rn)  *  (¥)p(rn+?i)  "  2p(rn)  +  p(rn“^)]  * 


where  m  denotes  the  number  of  elements  considered  in  0  <  r  <  D/2. 
The  subscript  n  refers  to  a  particular  location. 

Figure  H7  shows  the  system  considered  for  m  *  2.  The  finite- 
difference  approximations  of  Eq.  (28)  at  the  points  1  and  2  may 
be  written 


(33) 

(34) 
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~  (-P!+P2)  +  —  (Pj-Pj)  +  k2[l  -  (t-1o)2]  p,  *  0 


D 


(35) 


D 


~  (P3-2p2+P!)  +  ( P 3 ~P j )  +  k2[l  -  (T-ia)2]  p2  =  0  (36) 
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The  boundary  conditions  Eqs.  (30)  and  (31)  can  be  approxi¬ 
mated  by 
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The  coefficients  of  the  homogeneous  system  of  Eqs.  (35)  to  (37) 
vanish  for 
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For  |8Z/kDpc|  >>  1,  the  result 


t  —  ic  ~  ^1  "  1  5^  TT 


(39) 


agrees  with  the  general  approximation  for  propagation  constants 
in  ducts  of  arbitrary  shape,  with  weak  effects  of  the  absorbing 
wall  (Ref.  3 ) 


x  -  io  =  */i  -  i  5a  ££~ 

T  V  kS  Z  » 


(40) 


'■/here  Ua  denotes  the  absorptive  part  of  the  duct  perimeter,  S  the 
area  of  the  duct  cross  section. 


For  | 8Z/kDpc |  <<  1,  Eq.  (38)  yields 
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There  is  no  wave  propagation  in  a  duct  with  soft  walls  below 
the  cut-off  frequency  at  which  .Vl'-  \; 

i 1 

\\ 

kD  =  ) 

"T  ” 


This  value  closely  approximates  the  first  zero  of  the  zero-order 
Bessel  function. 

As  evident  for  the  extreme  cases  of  rigid  and  soft  duct  lin¬ 
ings,  Eq.  (38)  gives  a  good  approximation  for  the  propagation 
constant  of  the  fundamental  mode.  In  the  vicinity  of 


a  vanishing  square-root,  in  Eq.  (38)  indicates  the  highest  attenu¬ 
ation. 

ESTIMATION  SCHEME 

In  order  to  estimate  the  attenuation  provided  by  a  given 
duct,  one  may  proceed  as  described  below. 

1.  Determine  the  frequency 

fD  =  c (T)/D 

above  which  the  duct  may  be  considered  as  large  compared  to 
a  wavelength  (that  is,  D>1).  Here  c(T)  denotes  the  speed  of 
sound  in  the  gas  in  the  duct  (which  speed  is  proportional  to 
the  square-root  of  the  absolute  temperature). 

2.  Calculate  the  duct  diameter/length  ratio  D/L. 
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For  frequencies  f  >  f  p  (where  X<D), 

(a)  If  D/L<1,  find  the  attenuation  AL  from 

al  =  10  log  (yt jrrir)  •  (lla) 

(b)  If  D/L>1,  find  the  attenuation  from 

kL  ■  3.3  55“  10  fe  (fe)2  Re{^}  (l6a) 

For  frequencies  f  <  f ^  (where  X  >  D),  the  lengths  Lg  and  Lr 
of  the  lined  and  unlined  duct  portions  (if  any),  must  be 
considered,  in  addition  to  the  type  of  sound  source. 

Consider  a  duct  geometry  as  shown  in  Fig.  H3,  with  a  duct 
lining  configuration  as  shown  in  Fig.  H4.  Radial  partitions 
are  used  in  order  to  obstruct  axial  sound  propagation  inside 
the  lining.  The  absorptive  material  is  assumed  to  have  a 
rigid,  open-pore  structure;  the  acoustical  properties  re¬ 
quired  for  its  specification  are  the  flow  resistance  5b/pc, 
the  structure  factor  X  and  the  porosity  q p.  For  a  layer  of 
absorptive  material  that  is  thin  compared  to  the  wavelength, 
X  *  t'  «  1. 

(a)  Calculate  the  wall  Impedance  Z  of  the  lined  duct  from 
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and  k  denotes  the  wavenumber  in  free  space. 

(b)  Find  the  approximate  propagation  constants  c,t  of  the 
duct  by  use  of  Eq .  (38), (39),  or  (41),  whichever  is 
applicable,  depending  on  the  value  of  8Z/kDpc. 


(c)  Calculate  the  attenuation  that  results  over  the  length 
.  Lg  of  a  lined  duct  from 

ALL(dB)  =  8.7  koLs 


(d)  Calculate  the  attenuation  that  results  due  to  reflec¬ 
tions  at  the  end  of  the  duct  from 

AL£(dB)  =  10  log  [1  +  (2/kD)2] 

(e)  Determine  the  attenuation  due  to  interaction  between 
a  lined  section  and  an  unlined  tailpipe  from 

10,  if  entire  duct  is  lined 

third  term*  of  Eq.  (23),  otherwise 

(f)  Determine  the  attenuation  due  to  source  effects  from 

(third  term  of  Eq.(24),  for  low-impedance  source 
third  term  of  Eq.(25),  for  high-impedance  source 

(g)  Find  the  total  attenuation  from 
AL  =  ALl  +  ALe  +  ALj  +  ALg 


Sample  Results 

The  results  of  some  sample  calculations  are  plotted  in  Figs. 
H5  and H6. Figure  H5  pertains  to  a  duct  which  is  lined  over  a  length 
L  =  3D  and  which  has  an  urilined  tailpipe  of  the  same  length; 
the  lining  thickness  is  (d  =  D/2)  and  the  absorbing  material  is 
a  thin  layer  of  thickness  b  *  d/10  with  a  flow  resistance  that  is 
twice  the  characteristic  impedance  pc  of  the  gas.  The  attenuation 
8.7  koLs  over  a  length  Ls  of  a  homogeneous  lined  duct  and  the 
reflection  loss  10  log  [1  +  (kD/2)2]  of  the  tailpipe  are  plotted 
in  Fig.  H5,  Figure  H6  shows  the  excess  attenuation  AL  -  8.7  koLs, 
which  accounts  for  interactions  in  the  case  of  a  pc-source  imped¬ 
ance.  The  excess  attenuation  can  be  approximated  by  the  reflec¬ 
tion  loss  of  the  tailpipe  in  the  frequency  region  where  the  total 
length  Ls  +  Lr  exceeds  a  quarter  wavelength. 


*or  approximations ,  see  text 


Illustrative  Calculations 


High-Frequency  Effect 

Evaluate  the  effect  of  placing  a  13  ft  long  duct  of  1.65  ft 
exit  diameter  at  the  outlet  of  a  small  turbofan  engine,  if  the 
flow  leaving  the  duct  is  at  250°.  (This  duct  obviously  will 
affect  only  the  noise  emanating  from  the  engine  outlet;  a  duct 
at  the  inlet  will  be  required  to  reduce  the  noise  emanating  from 
that  end  of  the  engine.) 

The  speed  of  sound  at  70°F  is  1100  ft/sec;  thus  at  250°F, 
c  =  1100  =  1270  ft/sec.  Then 

fD  =  c (t )/D  =  (1270  ft/sec )/( 1. 6  ft)  =  770  Hz 

At  all  frequencies  above  770  Hz,  the  duct  is  large  compared  to 
the  acoustic  wavelength. 

Since  D/L  =  1.65/13  s  0.127,  Eq.  (11a)  applies.  If  the  duct 
is  unlined,  a  =  0  and  the  attenuation  provided  by  the  duct  is 

4L  -  10  ■  10  loe  ^5#  '  ^  dB 

at  all  frequencies  above  about  800  Hz.  (If  the  duct  contains  a 
lining  with  an  absorptive  coefficient  a  =  0.4,  one  obtains 
AL  =  8.5  dB.  The  value  of  a  generally  depends  on  frequency,  how¬ 
ever,  so  that  one  would  need  to  carry  .out  the  foregoing  calcula¬ 
tion  at  several  frequencies  in  the  range  of  interest.) 

Ducts  Muffling  only  One  of  Two  Like  Sources 

If  an  unmuffled  (and  unducted)  fan  produces  171  dB,re  10~12 
watts  at  17.5  kHz,  one  may  assume  168  dB  from  the  outlet  and  168 
dB  from  the  inlet.*  Adding  a  duct  with  AL  =  7  dB  to  the  outlet 
would  reduce  the  noise  from  that  source  to  161  dB,  but  would  not 
affect  the  inlet  noise  -  and  thus  would  result  in  a  level  of  169 
dB,re  in-12  watts  (which  corresponds  to  a  decrease  of  only  2  dB). 
However,  adding  a  duct  with  AL  =  7  dB  to  both  inlet  and  outlet 
would  reduce  the  171  dR  by  7  dB ,  i.e.  to  164  dB. 


*Bee  Fig.  1  for  combined  acoustic  levels. 
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Low-Frequency  Effect  (  n 

Consider  now  what  effect  the  previously  described  duct  would 
have  on  sound  at  150  Hz,  if  half  of  the  duct  is  lined  with  an 
acoustically  absorptive  system  having  a  wall  impedance  Z  =  2.5pc 
(at  150  Hz),  and  half  is  unlined  (see  Fig.  H5  for  general  con¬ 
figuration)  . 

Since  fp  =  770  Hz,  one  must  use  item  4  of  the  preceding  es¬ 
timation  scheme.  If  Z  is  given  (as  obtained  from  experimental 
data,  for  example)  one  need  not  calculate  it.  Since 

k  -  a)/c  =  2irf/c  =  2tt  (15  0/sec)/ (12  70  ft/sec)  =  0.743/ft, 
one  finds 


8Z/pc  _  8(2.5)  _  20.0 

kD  (0. 743/ft  V(  I.65  ft)  1.23 


16.3 


This  value  is  much  greater  than  unity;  therefore  Eq .  (39)  applies. 
With  (VkD)(pc/Z)  =  (4/1.23)  (1/2.5)  *  1.30  that  equation  becomes 

,  -  10=  [1-1(1. 30)]*  -  [Vi+u'.To'jV1  arc  tan  1-30]1'4 

=  [l.6»e-152-1,”],/S  -  1.28e-126-2° 

=  1.28(cos  26 . 2°-i  sin  26.2°)  =  1.15  -  10.565 
Then,  since  the  lined  length  is  Ls  =  6.5  ft. 


/ 

/ 


ALl  =  8.7koL$  =  8.7(0.743/ft)(0.565)  (6.5  ft.)  =  23.5  dB 
ALe  =  10  log[ 1+ (2/kD) 2 ]  =  10  logC 1+ ( 2/1.23) 2 ]  =  5.5  dB 


and  from  Eq.  (23),  letting  K  =  t  -  i c. 
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-2ikLR 

where  ^  =  1  -  C  +  5 -  (1  +  5) 

l+ikD/2 


one  finds  after  much  arithmetic  that  s  -1,08  +  1(2.18)  and 
ALj-  »  n  dB. 

For  a  high-impedance  source  (which  produces  essentially  the 
same  acoustic  pressure  with  and  without  the  added  duct  —  probably 
a  reasonable  approximation  for  engine  noise  sources)  one  finds 
from  Eq.  (25)  that 

e-2ikLR  1  _  e-2ikLs5 
l+ikD/2  25 

The  total  attenuation  at  150  Hz  therefore  is  23.5  +  5.5  -  0  +  4 
■  33  dB. 


ALg  =  20  log  1  - 
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APPENDIX  0 
GEAR  NOISE 


MECHANISM  OF  NOISE  GENERATION 

Tooth  Contact  Forces 

Gear  sets  transmit  torque  and  power  by  virtue  of  forces 
transmitted  via  mating  gear  teeth.  When  a  given  tooth  does  not 
make  contact  with  a  mating  tooth,  it  is  subject  to  zero  force; 
when  contact  is  made,  the  inter-tooth  force  increases,  reaches 
a  maximum,  decreases,  and  again  returns  to  zero  when  tne  tooth 
disengages.  These  dynamic  tooth  forces,  which  are  inherent  in 
the  operation  of  gears  (and  which  may  be  modified  or  aggravated 
by  inaccuracies  in  gear  tooth  profiles  and  gear  alignment), 
cause  the  gear  to  vibrate.  These  gear  vibrations  may  radiate 
sound  directly,  or  they  may  cause  vibrational  energy  to  travel 
via  the  shaft  and  bearings  to  the  gearbox,  which  may  then  radiate 
sound  into  the  surrounding  air. 

Characteristics  of  Noise 

One  may  deduce  some  qualitative  properties  of  the  gear  noise 
spectrum  by  considering  the  dynamic  tooth  interaction  forces. 

The  force  acting  on  a  gear,  and  thus  the  gear  noise,  may  be  ex¬ 
pected  to  have  a  periodic  component  associated  with  the  tooth 
contact  frequency  ft.  In  addition,  one  would  expect  periodic 
components  at  the  shaft-rotation  frequencies  of  the  various  inter 
acting  gears,  due  to  any  asymmetries  and  other  geometric  inaccura 
cies  of  the  gears. 

For  a  geometrically  perfect  gear  set,  the  force  spectrum 
(and  therefore  also  the  sound  spectrum)  would  have  components  at 
only  the  tooth  contact  frequency  and  Its  harmonics,  with  the  am¬ 
plitude  of  these  spectral  components  depending  on  the  shape  of 
the  force-time  curve  corresponding  to  a  single  tooth  contact. 

The  amplitudes  of  these  components  may  be  expected  to  be  approxi¬ 
mately  constant  up  to  the  frequency  which  corresponds  to  the  re¬ 
ciprocal  of  the  effective  duration  Teff  of  a  contact  force  pulse 
(see  Fig.Jl);  beyond  this  frequency,  the  amplitudes  are  expected 
to  decrease  rapidly. 


For  an  imperfect  gear,  or  with  unsteady  driving  moment,  the 
shapes  -  as  well  as  the  amplitudes  -  of  the  tooth-force  pulses 
may  vary  with  a  period  that  corresponds  to  the  shaft-rotation 
frequency  fs,  as  sketched  in  Fig.  Jl.  Because  of  this  "amplitude 
modulation"  of  the  tooth  force  pulses,  the  spectra  here  may  be 
expected  also  to  contain  components  at  the  shaft  rotation  fre¬ 
quency  and  at  harmonics  of  that  frequency.  Measured  spectra  of 
gear  noise  (e.g.,  see  Fig.  10  of  Ref.  1)  are  in  agreement  with 
the  foregoing  qualitative  considerations. 

DATA 

Available  Data 

Although  much  work  has  been  done  on  gear  noise,  only  quali¬ 
tative  understanding  of  it  rests  on  a  relatively  firm  basis. 

The  present  state  of  the  technology  does  not  permit  one  to  predict 
the  noise  of  a  gear  set  from  the  basic  physical  parameters  of  that 
set  (except  from  correlations  of  empirical  data),  and  generally 
also  does  not  permit  one  to  predict  with  much  confidence  the  ef¬ 
fects  of  design  changes  on  noise.  A  number  of  relevant  recent 
publications  are  indicated  in  the  appended  reference  list. 

The  only  available  collection  of  data  on  the  noise  of  many 
different  gears,  measured  under  well-known  conditions,  appears  to 
tie  that  of  Refs.  1  and  2.  These  references  summarize  the  results 
of’  noise  measurements  on  76  sets  of  power  gears  (including  spur 
and  bevel  gears  and  planetary  sets)  transmitting  between  10  and 
25,000  kw  (7.5  and  18,500  Hp)  at  peripheral  speeds  between  3.6 
and  160  m/sec  (12  and  525  ft/sec)  with  gearing  ratios  up  to  256. 

Dependence  of  Overall  Noise  on  Transmitted  Power  and 

Tooth  Force 

A  careful  statistical  analysis  of  the  measured  gear  noise 
data  revealed  a  strong  correlation  (0.87  correlation  coefficient) 
between  the  total  acoustic  power  produced  by  a  gear  set  and  the 
mechanical  power  transmitted  by  that  set.  A  similar  correlation 
was  found  between  noise  and  the  average  tooth  force. 

Figure  J2  (taken  from  Ref.  2)  shows  how  the  measured  overall 
sound  power  levels  of  gear  sets  vary  with  the  transmitted  power 
(which  for  these  measurements  was  between  75%  and  100%  of  the 
design  power  for  the  sets  tested).  Also  indicated  in  that  figure 
i n  line  which  corresponds  to  the  mean  of  the  measured  points, 
a.-,  well  as  dashed  lines  which  correspond  to  one  standard  devia¬ 
tion  from  the  mean.  An  equation  for  the  mean  line  is  also  given 
In  the  figure. 


Noise  Spectra 


Octave-band  noise  spectra,  derived  from  Refs.  1  and  2, 
are  presented  in  F1rs.J3.J4  andJ5  for  the  three  classes  of 
gears  and  speeds  treated  in  the  original  references.  Fipures  J3 
and  J4.  which  may  be  seen  to  be  very  similar,  pertain  to  plane¬ 
tary  and  spur  gear  sets  operating  at  relatively  high  speeds; 
Fig.J5  pertains  to  gears  of  various  types  operating  at  lower 
speeds. 

The  spectra  in  Figs,J3-5  are  presented  in  terms  of  differ¬ 
ences  between  overall  levels  and  octave  band  levels.  In  addi¬ 
tion  to  spectra  corresponding  to  the  means  of  the  data,  curves 
are  given  which  correspond  to  the  mean  +  standard  deviation. 

PREDICTION  CURVES 

By  combining  the  average  dependence  of  overall  acoustic 
power  on  transmitted  mechanical  power,  as  given  in  Fig.J2,  with 
the  spectra  of  Figs.  J3-5  one  may  arrive  at  octave-band  spectra 
of  sound  power  referred  to  mechanical  power. 

Figure  J6  shows  mean  spectra  obtained  in  this  manner,  and 
is  readily  applicable  for  general  noise  estimation  purposes. 
Because  of  the  similarity  of  the  mean-curves  of  Figs,J3  and  J4, 
these  two  have  been  combined  (averaged  and  smoothed)  into  a 
single  "high-speed"  gear  noise  curve;  the  "low-speed"  curve  is 
obtained  directly  from  Fig,  J5. 

ESTIMATION  SCHEME 

In  order  to  estimate  the  noise  produced  by  a  gear  set  (two 
mating  gears),  one  may  proceed  as  follows: 

1.  Select  the  curve  of  Fig.  J6  that  corresponds  to  the  rota¬ 
tional  sneed  (rpm)  of  the  smaller  (:.  L.,er-sneed)  gear. 

2.  Calculate  10  log(Hp),  where  Hp  is  the  mechanical  horse¬ 
power  transmitted  by  the  gear  set,  and  add  this  value  to 
the  numbers  shown  along  the  vertical  scale  of  Fig.  J6. 

The  resulting  numbers  then  represent  the  octave-band 
sound  power  levels  for  the  gear  set. 

3.  Calculate  the  tooth  contact  frequency  from 

ft(Hz)  =  Nt(rpm)/60 

where  Nt  is  the  number  of  teeth  on  either  gear  and  rpm  is 
the  rotational  speed  (in  revolutions/min )  of  that  same  gear. 


Then,  for  aural  detection  estimation,  take  the  levels  (ob¬ 
tained  from  the  foregoing  steps)  at  frequencies  ft,  2ft, 
3ft  to  correspond  to  pure  tones,  and  the  levels  for  octave 
bands  above  3ft  to  correspond  to  broad-band  noise.  Delete 
all  parts  of  the  spectrum  of  Fig.  J 6  below  ft.  (No  appre¬ 
ciable  noise  is  produced  at  these  low  frequencies.) 

Illustrative  Calculation 

Consider  a  gear  set  transmitting  100  Hp ,  with  the  smaller 
gear  (which  has  18  teeth)  rotating  at  3000  rpm.  Here 

10  log(Hp)  =  20 

f.  =  N.  (rpm)/60  =  l8(3000)/60  =  900  Hz 
X>  t/ 

This  is  a  "high-speed"  gear  set  (rpm  >  1500);  the  dashed  curve 
of  Fig.  J6  applies.  For  aural  detection  estimation,  one  thus 
obtains  the  following  pure  tone  noise  levels: 


f  (Hz ) 

200 

1800 

2700 

Lw(dR,re  10“  12  watts) 

91 

89 

86 

olus  the  following  octave-band  levels: 


f  (Hz) 

4000 

8000 

Lw(dR,re  10” 12  watts) 

85 

78 

with  no  significant  noise  below  900  Hz. 
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APPENDIX  K 


SOUND  RADIATED  FROM  FLOW  OVER  RIGID  SURFACES 


SOUND  FROM  TURBULENT  BOUNDARY  LAYER 

Infinite  Surface 

The  sound  power  radiated  to  the  atmosphere  (per  unit  area) 
from  an  infinite  turbulent  boundary  layer  supported  by  a  rigid 
flat  surface  provides  a  lower  bound  for  the  sound  power  radiated 
by  a  finite  system.  Because  of  this  fact,  and  because  analysis 
of  the  infinite  case  is  relatively  simple  and  exhibits  the  im¬ 
portant  parameters,  presentation  of  a  summary  of  the  analysis  is 
justified  here. 

Mathematical  Model  of  Boundary  Layer 

Extensive  experimental  and  theoretical  work  has  been  done  to 
obtain  descriptions  of  the  pressures  on  surfaces  under  turbulent 
boundary  layers  (Refs.  1-6).  These  pressures  usually  are  des¬ 
cribed  in  terms  of  space-time  correlation  functions 

4>(xj,Xj,t)  =  <P(x10  ,x30  ,t0  )p(x,  0  +  Xj  ,x30+x3  ,t0+t)>  (1) 


where  x10  and  x30  represent  coordinates  on  the  surface,  x]0  in 
the  direction  of  the  flow  and  x30  perpendicular  to  the  flow  di- 
r  ,.ion;  x3  and  x3  represent  the  separation  between  two  observa¬ 
tion  noints,  the  brackets  <...>  indicate  averaging  over  x10>x30> 
and  time  t0 .  One  may  calculate  the  wavenumber- frequency  spectrur 
of  the  cressure  as: 


=  — ~  fff  *  (x)t)e"i(-o‘-“u)t;  dx  dt  (2) 

p  ^  ( 2  7T )  3  JJJ  P 

where  kp  is  the  wavenumber  vector,  with  the  components  k3  and  k3, 
and  x  is  the  separation  vector,  with  components  x3  and  x3,  and  w 
represents  the  radian  frequency. 

It  has  been  shown  (Ref.  1)  that  the  observed  behavior  of  the 
wavenumber-frequency  spectrum  may  be  represented  by 


$  (k  ,u>)  =  4>f(w)  4>d  k j  ( uj )  - 


where  4>f(w)  represents  the  frequency  spectrum  one  measures  with  a 
transducer  that  is  fixed  at  one  location  on  the  surface,  ‘t’dCkj) 
describes  the  decay  of  eddies  in  the  flow  direction,  and  4>3(k3) 
is  the  wavenumber  spectrum  in  the  cross-stream  direction. 

The  two  wavenumber  spectra  typically  obey 


W  =  <Ld/Tr)(1  +  k?Ldrl 

(4) 

*3(k,)  -  (L3/n)(l  +  kjL ^ )-1 

(5) 

where  Ld  represents  the  eddy  decay  length,  which  may  be  approxi¬ 
mated  by 

Ld  *  9Uc/w  (6) 

and  L3  is  the  transverse  correlation  length,  which  may  be  approxi¬ 
mated  as 


La  *  U  /0.7w  (7) 

1  c 

in  terms  of  the  convection  velocity  Uc,  which  is  approximately 
0.7  times  the  free-stream  velocity  U. 

The  "fixed  transducer"  spectrum  may  be  determined  from 

the  corresponding  normalized  spectrum, xwhich  is  given  in  Fig.  K1 
and  further  discussed  below. 

Sound  Radiation 

Only  those  fluctuating  pressure  components  for  which  the 
wavenumbers  kD  are  smaller  than  the  acoustic  wavenumber  k  =  w/c 
(where  c  represents  the  speed  of  sound)  contribute  to  the  sound 
radiation.  Accordingly,  the  spectrum  of  the  radiated  acoustic 
pressure  at  the  wall  is  given  by 

♦ac<«>  -  VU)  /  f  *d(kl  -  iJ-K(k3)dkldk3  *  (6) 

(k j  +k3  <k)  c/ 


1  'id 


If  one  substitutes  Eqs.  (4),  (5)»  (6)  and  (7)  into  Eq.  (8),  and 
limits  oneself  to  the  usual  case  of  interest  here  where  the  con¬ 
vection  velocity  Uc  is  much  smaller  than  the  speed  of  sound 
(Uc  <<  c),  one  obtains  •  ' 


*ac(u))  *  *f(u))  J  J  lorfe.  dkidk3 


(k  +k  <k) 

13 


~  $f(w) 


6 . 3tt  2  u» 2 


irk2  . 


’With  Uc  *  0.7U,  this  result  reduces  to 


<*ac(w)  *  2 .  b  ■<  2  0-2M2  $ ^(l})  (101 

where  M  =  U/c  represents  the  Mach  number  of  the  free-stream  flow. 

Equation  (10)  gives  the  acoustic  pressure  spectrum  in  terms 
of  the  boundary  layer  pressure  spectrum.  From  this  result  one 
may  readily  obtain  the  relation 


p2/pbl  *  2"5  x  10_2m2 


which  holds  between  the  mean-square  acoustic  pressure  p2  in  a 
given  frequency  band  Aw  and  the  mean-square  pressure  p£]_  in  the 
boundary  layer  in  the  same  band;  these  pressures  are  related  to 
the  spectra  as 


*  /  4ac(")d“  • 

Aw 


=  /  *fu 


)  du>  . 


The  acoustic  power  W  radiated  from  a  surface  of  area  A  may 
be  estimated  from 

V/  =  Ap2/4pc  ,  (1 

where  p  represents  the  density  of  the  air,  as  usual. 


Frequency  Distribution 

The  fixed-transducer  spectrum  of  boundary  layer  pressure 
fluctuations,  which  also  determines  the  spectrum  of  radiated 
sound  in  view  of  Eq.  (10),  may  be  determined  from  the  normalized 
fixed-transducer  spectrum  <t^(u))U/<5 *  shown  in  Pig.  K1  as  a  function 
of  the  dimensionless  frequency  (Strouhal  number)  S  =  u)6*/U,  where 
6*  represents  the  displacement  thickness  of  the  boundary  layer. 

The  spectrum  values  are  adjusted  so  that  <t>J.(co)  satisfies  the 
relation  1 


$l  (<j))du 


1 


(1*0 


then  the  actual  spectrum  4>f(w)  is  related  to  the  normalized  spec¬ 
trum  <t^.(ui)  as 


Vw)  =  p0A$f(w)  (15) 

in  terms  of  the  overall  mean-square  fluctuating  pressure  p*A  in 

the  boundary  layer.  From  Fig.  K1  and  Eq .  (15)  one  may  deduce  how 
the  fluctuating  pressure  p^  in  third-octave  bands  varies  with 

frequency;  the  result  is  indicated  in  Fig.  K2 . 

The  boundary  layer  displacement  thickness  6*  is  a  function 
of  the  Reynolds  number  and  Mach  number.  For  low  Mach  numbers, 
and  for  Reynolds  numbers  in  the  range  between  106  and  107,  which 
is  of  primary  interest  for  small  aircraft,  one  may  estimate  6* 
from 


6#  z  0.0025  X  (16) 

where  X  is  the  distance  from  the  front  (leading  edge)  of  the  body 
of  concern  (Ref.  7). 

In  order  to  estimate  the  actual  levels,  one  still  requires 
the  mean-square  overall  boundary  layer  pressure  p^.  There 

exists  much  evidence  (Ref.  8)  that  the  corresponding  root-mean- 
square  pressure  is  proportional  to  the  free-stream  dynamic  pres¬ 
sure  q.  The  constant  of  proportionality  depends  on  Reynolds  num¬ 
ber,  Mach  number,  and  surface  roughness,  but  for  most  smooth  air¬ 
craft  surfaces  one  may  take 

t  0.006  q  .  (17) 


*  '  '  1 


Effects  of  Surface  Edges 


Whereas  pressures  acting  on  a  rigid  surface  can  give  rise  to 
no  appreciable  acoustic  particle  velocities,  pressures  near  an 
edge  do  not  encounter  the  high  impedance  of  the  surface  and  thus 
cause  greater  velocities,  and  therefore  more  noise  radiation. 

Thus,  the  noise  associated  with  edges  is  likely  to  dominate  over 
the  noise  from  boundary  layers  on  relatively  rigid  surfaces.  The 
noise  associated  with  edges  is  discussed  in  the  following  section. 

SOUND  FROM  FLOW  PAST  FINITE  BODIES 

Airfoil  in  Large-Scale  Turbulent  Flow 

A  body  moving  through  a  turbulent  fluid  experiences  unsteady 
lift  and  drag  forces;  the  related  reaction  forces  on  the  fluid 
are  responsible  for  sound  radiation.  In  most  practical  cases  of 
interest,  e.g.,  of  an  aircraft  flying  through  atmospheric  turbu¬ 
lence,  the  spatial  scale  of  the  pressure  fluctuations  is  consid¬ 
erably  larger  than  the  airfoil  dimensions.  There  then  occur  lift 
and  drag  fluctuations,  which  act  essentially  like  acoustic  point 
dipole  sources.  (The  point  source  approximation  is  valid  because 
the  associated  acoustic  wavelengths  are  much  greater  than  the 
characteristic  surface  dimensions.) 

The  acoustic  power  W  radiated  by  a  point  dipole  at  frequency 
f  is  given  by 


W  =  irF2f2/3pc3 


(18) 


where  F2  represents  the  mean-square  fluctuating  force,  p  the  den¬ 
sity  of  the  ambient  air,  and  c  the  speed  of  sound  in  air. 


One  may  express  the  root-mean-square  fluctuating  lift  or 
drag  force  in  terms  of  the  corresponding  lift  or  drag  coefficient 

CI,(or  D)  a‘; 


FL(or  D)  ~  2pu2/1‘CL(or  D)(urms/U) 


(19) 


where  U  represents  the  free-stream  (mean)  velocity  ,  A  a  reference 
area,  and  urms  the  root-mean-square  fluctuating  velocity.  If  one 
knows  the  turbulence  spectrum  urms(f)/U,  one  may  then  calculate 
the  corresponding  force  spectra,  and  from  these,  the  associated 
acoustic  power. 


Measurements  by  Clark  and  Ribner  (Ref.  9)  have  verified  the 
proportionality  of  sound  radiated  from  an  airfoil  in  large-scale 
turbulent  flow  to  the  mean-square  fluctuating  force  on  the  air¬ 
foil,  and  Sharland  (Ref.  10)  has  observed  that  the  directivity 
pattern  of  the  sound  radiation  from  such  an  airfoil  is  like  that 
of  a  classical  dipole.. 

Small-Scale  Turbulent  Flow  over  Airfoil 

If  the  disturbances  on  the  airfoil  and  the  associated  acous¬ 
tic  wavelengths  A  are  small  compared  to  the  airfoil  dimensions 
(e.g.  the  chord  b),  one  no  longer  has  the  effect  of  point  dipoles 
Instead,  for  A  <  b,  one  has  essentially  arrays  of  unsteady  forces 
acting  at  the  various  edges  (Refs.  11,  12). 

Chanaud  and  Hayden  (Refs.  11-13)  have  considered  two  ceparat- 
trailing-edge  noise  sources*:  (1)  Interaction  of  turbulent  bound¬ 
ary  layer  with  the  edge,  and  (2)  forces  due  to  wake  vortices  act¬ 
ing  on  the  edge.  By  analysis  of  the  problem  in  terms  of  dimen¬ 
sionless  groups  of  parameters  and  application  of  empirical  data, 
they  have  derived  the  following  prediction  expression  for  the 
overall  acoustic  power  level: 


Lw(0A)  ~  X  +  10  log  ( <5wU6 )  -  5 


(20) 


where  Lw(0A) 
w 
u 
6 

X 


overall  acoustic  power  level  (dB,  re  10“12  watts) 

spanwise  dimension  of  edge  (ft) 
free-stream  velocity  (ft/sec) 

boundary  layer  thickness  or  wake  thickness  (ft) 

(-27  dB  for  boundary  layer/edge  interaction 
(-23  dB  for  wake  vortex  effect. 


The  corresponding  octave  band  spectra  (Ref.  13)  may  be  obtained 
from  Pigs.  K3  and  K4. 

The  turbulent  boundary  layer  thickness  6  at  the  trailin'? 
edge  of  an  airfoil  depends  on  the  angle  of  attack,  the  surface 
condition  of  the  airfoil,  and  the  smoothness  of  the  inflow. 


*No  corresponding  leading  edge  studies  have  been  undertaken,  but 
because  the  pressure  fluctuations  at  the  trailing  edge  generally 
exceed  those  at  the  leading  edge  when  the  inflow  is  undisturbed, 
trailing  edge  noise  tends  to  dominate. 


Unless  one  has  better  information,  one  may  estimate  6  from  the 
corresponding  value  for  a  smooth  flat  plate,  for  which 

6  «  0.^4  bV5  (v/U)1/4  (21) 

where  5  and  the  chord  b  are  in  units  of  ft,  v  represents  the 
kinematic  viscosity  of  the  air  (ft2/sec)  and  U  is  the  free-stream 
velocity  (ft/sec). 

The  wake  thickness  6^  also  depends  strongly  on  the  condition 

of  the  airfoil.  For  thin  tapered-edge  airfoils  at  small  angles 
of  attack  one  may  estimate  6^  as  the  sum  of  the  airfoil  edge 

thickness  and  26.  For  airfoils  at  considerable  angles  of  attack 
(say,  | oc  |  >5°),  one  may  estimate  6W  from 

6^  =  b  sina  (22) 

where  b  is  the  chord  length.  Note  that  narrow-band  wake  noise 
seldom  occurs  for  |a|  >5°,  and  only  boundary  layer/edge  inter¬ 
action  noise  is  significant  for  such  cases. 

ESTIMATION  SCHEME 

Flow  Over  Fuselage 

In  order  to  estimate  the  broad-band  noise  produced  by  turbu¬ 
lent  flow  over  a  fuselage  (or  similar  body),  one  may  proceed  as 
follows : 

.1 .  Divide  the  fuselage  into  a  number  of  convenient  regions  over 
which  the  boundary  layer  is  relatively  uniform. 

2.  For  each  region, 

(a)  Determine  the  boundary  layer  displacement  thickness  6* 
by  using  available  aerodynamic  data  or  estimating  on 
the  basis  of 

6*  =  0.0025  X 

where  X  represents  the  distance  from  the  aircraft  nose 
to  the  middle  of  the  region.  Then  calculate  the  actual 
frequencies  f  that  correspond  to  the  reduced  frequen¬ 
cies  f6*/U  shown  along  the  horizontal  axis  of  Fig.  K2. 


(b)  Determine  the  corresponding  octave  band  sound  power 
ievei  Lw(qb)  by  adding* 

51  +  20  log  q  (lb/ft2)  +  10  log  A  (ft2) 

to  the  numbers  shown  along  the  vertical  coordinate  of 
Fig.  K2 . 

4.  Combine  the  octave-band  levels  for  all  regions  (band  by  band), 
by  the  method  discussed  in  the  main  text  and  Fig.  1  of  the 
main  text . 

Flow  Over  Ai rfoi 1 s 

In  order  to  estimate  the  noise  produced  by  flow  over  airfoils, 
proceed  as  follows: 

1.  Find  the  (broad-band)  noise  due  to  turbulent  flow  by 

(a)  Estimating  the  boundary  layer  thickness  <$  at  the  trail¬ 
ing  edge,  either  from  aerodynamic  data  or  by  use  of 

Eq.  (21) 

(b)  Calculating  the  overall  acoustic  power  level  from  Eq . 
(20),  with  x  =  -27  dB 

(c)  Determining  the  octave-band  spectrum  by  use  of  Fig.  K3. 

2.  Find  the  (essentially  pure-tone)  noise  associated  with  vor¬ 
tex  shedding  by 

(a)  Estimating  the  wake  thickness,  either  from  aerodynamic 
data  or  from 

!6te  +  26  for  |a|  <  5° 
b  sina  for  |a|  >  5° 


*By  combining  Eqs .  (11),  (13)  and  (17),  taking  p  =  0.07  lb/ft3 
and  c  =  1100  ft/sec,  and  using  the  standard  definitions  and  re¬ 
ference  value  for  power  level,  one  finds  that 


Lw  =  10  log  [1.3  *  105  q2  ( Ib/ft 2 )  A  (ft2)]  +  20  log  [p03/r«0Al 


where  6W  represents  the  airfoil's  trailing  edge  thick¬ 
ness.  (Note  that  vortex  noise  occurs  only  very  rarely 
for  | a |  >  5°.) 

(b)  Calculating  the  overall  acoustic  power  level  from  Eq. 
(20),  with  x  =  -23  dB. 

(c)  Determining  the  octave-band  spectrum  by  use  of  Pig.  K4 . 
(For  aural  detection  estimation,  use  the  value  at 
f<5w/U  =  0.2  as  a  pure  tone,  and  consider  remainder  of 

spectrum  as  broad-band  noise.) 

3.  If  the  (predominantly  low-frequency)  noise  due  to  flight 

through  large-scale  turbulence  is  of  concern,  estimate  this 
noise  for  a  specified  turbulence  spectrum  u  (f)/U  by 

X  ITlo 

(a)  calculating  the  fluctuating  lift  and  drag  forces  from 
Eq.  (19) 

(b)  using  Eq.  (18)  and  Lw  =  10  log  (W/Wpef)  to  determine 
the  sound  power  level  at  each  frequency. 
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FIGURE  K 1  NORMALIZED  FIXED-TRANSDUCER  SPECTRUM  OF  PRESSURE 
FIELD  UNDER  TURBULENT  BOUNDARY  LAYER 
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APPENDIX  L 


SOUND  RADIATION  FROM  BEAM-REINFORCED 
PLATES  EXCITED  BY  POINT  FORCES 

INTRODUCTION 

An  unsteady,  force  applied  to  an  aircraft  structure  (e.g., 
to  a  stringer  or  directly  to  the  skin),  sets  this-  structure  into 
motion,  and  such  motions  lead  to  the  radiation  of  sound.  This 
section  presents  an  approach  toward  estimating  th.e  sound  radi¬ 
ated  from  typical  aircraft  structures  ex.cited  by  a  normal  point 

force . 

/ 

Aircraft  surface  structures'  typically  consist  of  a  thin 
skin,  reinforced  at  intervals  by  stringers,  rings,  frames,  etc. 
For  the  present  approximate  prediction  purposes  it  suffices  to 
consider  all  such  structures  as  thin  flat  plates  with  straight 
reinforcing  beams,  and  to  address  the  two  simplest  and  most  sig¬ 
nificant  corresponding  sound  radiation  problems.  Accordingly, 
the  following  paragraphs  deal  with  sound  radiation  from  beam- 
reinforced  plates  excited  by  oscillating  localized  (point) 
forces,  which  act  normal  to  the  piano  of  the  plate,  either  on 
the  skin  or  on  a  reinforcing  beam. 

POWER  INPUT  FROM  POINT  FORCE 

Admittance  of  Beam-Plate  System 

The  point  incut  admittance  Yc  of  a  system  which  consists  of 
an  infinitely  long  uniform  reinforcing  beam  connected  continu¬ 
ously  all  along  its  length  to  a  uniform  plate  of  infinite  extent 
is  given  (Ref.  1)  by 

(3-r2)co5~‘r  +  ( 3+r2 ) sink"1 r  ^  ^  ^ 

2r3 ( 1 -r2 ) 1/2  J 

Tills  expression  applies  for  the  simplest  case  In  which  the  ex¬ 
citing  force  is  so  located  and  the  beam  section  is  of  such  a 
shape  that  the  beam  does  not  twist  as  it  flexes.  The  first  term, 

Yb  =  C2mbcb(l+i) ]"‘  .  (2) 


i  . 
TTB^u)mb 


...  A# 


represents  the  point  input  admittance  of  the  beam  in  absence  of 
the  plate,  whereas  the  rest  represents  the  contribution  due  to 
the  plate.  The  symbols  employed  in  the  foregoing  relations  have 
the  following  definitions: 

w  -  2irf  =  radian  frequency 
rnp  =  mass  per  unit  length  of  beam 
Gj  =  B/D  »  12I/h3 
r  =  cp/cb  =  /h7^ 

B  =  El  =  bending  stiffness  of  beam 
^  D  *  Eh3/12  =  flexural  rigidity  of  plate 
E  =  Young's  modulus 

I  =  centroidal  moment  of  inertia  of  beam  cross-section 
h  =  plate  thickness 

Cp  =  [a)hcjj//l2’]  ^  speed  of  bending  waves  in  plate 
cb  =  v/wcxIcL  =  sPee<*  of  hendin-  waves  in  beam 
CL  =  ^/oa  ~  longitudinal  wave  velocity 
Pg  =  density  of  material  of  beam  and  plate 
ctj  =  radius  of  gyration  of  beam  cross-section 

Admittance  of  Plate 

The  point  input  admittance  of  an  infinite  plate  is  given 
(Ref.  2)  by 


where  mp  is  the  mass  per  unit  area  of  the  plate.  This  expression 
also  represents  the  frequency-average  of  the  admittance  of  a 
finite  plate,  provided  that  the  averaging  interval  encompasses 
several  resonances  (Ref.  3). 

Power  Input 

The  power  W^n  supplied  to  a  system  by  an  oscillating  force 
of  amplitude  P  depends  only  on  the  real  part  of  the  input  admit¬ 
tance  Y, 

i  Wln  =  i|F|2Re(Y)  j  ('0 


lb  ^ 


—  hzpcT 

V5  LJ 


hence  one  obtains 


.  |2  o  for  excitation  on  beams  of 

I  P I  //°mbcb  beam-plate  systems 

(5) 

|  P|  2/l6/D'm'p  for  excitation  plates 

EXCITATION  ACTING  ON  BEAM  OF  BEAM-PLATE  SYSTEM 
Components 

The  sound  radiation  here  may  be  considered  as  composed  of 
three  parts : 

(1)  radiation  from  the  vicinity  of  the  excitation  point  (near¬ 
field  radiation), 

(2)  radiation  from  waves  propagating  along  the  beam,  and 

(3)  radiation  from  reverberant  vibration  field  on  the  plate, 

Near-Field  Radiation 

The  acoustic  power  radiated  by  the  near-field  in  the  vicin¬ 
ity  of  the  excitation  point  may  be  approximated  (Ref,  7)  by 

W ,  *  0 . 3^ ( A  +w IpcX^v2  (6) 

where  «  vhc,//3  c  =  flexural  wavelength  of  plate  at  critical 

frequency 

=  cb/f  =  bending  wavelength  of  beam 
w  =  width  of  beam 
P  =  density  of  ambient  air 
c  =  speed  of  sound  in  ambient  air 
v  =  velocity  amplitude  at  excitation  point. 

V.  ii.ii  v  -  P  |  Y  c  |  and  Yc  as  given  by  Eq .  (l),one  may  then  readily 
determine  the  power  . 

Radiation  from  Waves  Propagating  along  Beam 

because  reinforcing  beams  typically  are  strongly  coupled  to 
the  plate,  energy  in  waves  travelling  along  the  beam  tends  to  be 
1.  *  *  .•  i  1 1 :  •  m  ittod  Into  the  plate  within  a  relatively  short  distance. 

■*’!,'  radiation  from  these  propagating  naves  tends  to  be  neg- 

!  i i  bl<- . 


Radiation  from  Plate;  Estimation  Scheme 

The  plate  flexural  wavenumber  in  the  direction  parallel  to 
the  beam  is  equal  to  the  wavenumber  k^  =  oj/c^  for  waves  on  the 
beam.  This  wavenumber  component  is  small  compared  to  the  (total) 
wavenumber  kp  =  w/cp  for  free  bending  waves  on  the  plate.  There¬ 
fore,  free  plate  flexural  waves  propagate  in  directions  nearly 
perpendicular  to  the  beam  (Ref.  4). 

If  the  plate  is  infinite  and  if  cp  <  c,  the  free  flexural 
waves  radiate  no  sound  at  all;  if  cp  >  c,  they  radiate  well. 
However,  the  plates  in  actual  aircraft  structures  cannot  be  con¬ 
sidered  infinite,  so  that  infinite-plate  results  here  are  rela¬ 
tively  meaningless. 

As  plate  waves  propagate  along  a  finite  structure,  they  soon 
encounter  beams,  where  they  are  partly  reflected  and  partly 
transmitted.  In  any  bay  between  reinforcing  beams  the  initial 
and  multiply-reflected  waves  build  up  a  reverberant  wave  field, 
in  which  amplifications  can  occur  as  the  result  of  resonance  ef¬ 
fects,  and  where  sound  is  radiated  because  of  interaction  of  the 
plate  waves  with  the  boundaries.  In  order  to  analyze  this  case 
most  simply,  one  may  make  the  very  reasonable  assumptions  that 
all  of  the  power  supplied  to  the  structure  is  transmitted  to  the 
plate,  and  that  in  the  steady  state  the  power  input  to  the  plate 
must  equal  the  power  lost  by  the  plate. 

The  loss  of  power  from  a  finite  panel  may  be  considered  as 
consisting  of  two  parts:  acoustic  radiation  (W)  and  mechanical 
dissipation  (W^) .  Hence  one  may  write  an  energy  balance 

Wj,n  =  2W  +  =  A<vz>(mpU)n  +  2pco)  (7) 

where  A  =  area  of  plate  (one  side) 
n  =  loss  factor  of  plate 
<v2>  =  mean-square  velocity  of  plate 
0  =  acoustic  radiation  efficiency. 


The  factor  2  appears  in  the  above  equation,  because  acoustic  en¬ 
ergy  generally  is  radiated  from  both  sides  of  a  plate.  The  acous¬ 
tic  power  radiated  from  one  side  may  then  be  written  as 


V.'  = 


A<vz>pco  * 


W 


in 


2  + 


mpu)n 

pea 


\  ^  > 


W^n  may  It 


its, 


where  the  second  form  was  obtained  by  use  of  Eq.  (7). 
calculated  from  Eq.  (5). 


EXCITATION  ACTING  ON  PLATE 


If  a  force  acts  on  a  finite  plate,  there  again  results  a 
reverberant  flexural  wave  field  on  the  plate,  as  previously  dis¬ 
cussed.  Whereas  the  plate  wave  field  due  to  excitation  on  a  beam 
tends  to  consist  primarily  of  wave  components  propagating  perpen¬ 
dicular  to  the  beam,  tnat  due  to  point  excitation  directly  on  the 
plate  tends  to  be  more  homogeneous  in  direction. 

However,  the  same  energy  balance  considerations  apply  here 
as,  before,  so  that  Eq.  (8)  applies  for  the  present  case  as  well; 
the  only  difference  is  that  W^n  and  ara(j  are  different  for  the 
two  excitation  conditions.  W^n  here  is  given  by  Eq.  (5)>  and 
°rad  here  is  computed  on  the  basis  of  the  entire  plate  perimeter 
(see  Pig.  LI). 

ESTIMATION  OF  LOSS  FACTOR  AND  RADIATION  EFFICIENCY 

Loss  Factor 

In  spite  of  considerable  effort  that  has  been  expended  on 
t,ho  problem,  there  still  exists  no  reliable  general  means  for 
predicting  the  loss  factor  n  of  realistic  aircraft  structural 
panels.  Pome  such  rediction  means  are  suggested  in  Ref.  5,  but 
structure:-  whir:,  are  as  complex  as  those  in  actual  aircraft, 
i  predict )  -in  menu.',  are  not  much  more  reliable  than  simply 

a. ■.•liming  n  50  0.11. 

Radiation  Efficiency 

'!'! i* ■  acoustic  radiation  efficiency  a  ,.afj  of  t  eam-reinforced 
;  I  ■  •  t  f  ■  ■  bus  been  analyzed  in  Ref:  .  3  and  '(>  and  is,  also  discussed 
'  i.  i.i  nooks  (Refs,  A ,  7).  Figure  LI  summarizes  the  dependence 
!’  o of  a  firi’te  pane]  on  fre  quency  and  the  various  panel 
:  'r:it':-  k,.-rs  and  may  ne  used  to  estimate  a  ra  . 

ESTIMATION  SCHEME 

in  order  to  estimate  the  sound  radiated  from  a  panel  excited 
iv  an  ,  >sc i L  I  a tor. v  point  force  of  amplitude  |F|  acting  either  di- 
:  ■ '  •  c  1 1 y  on  the  panel  or  on  a  reinforcing  beam  (or  rib)  at  the 
pan-  i's  edg,  ,  one  may  proceed  as,  follows: 

i.  b.u-  Pi  1,1  to  estimate  the  udiutiori  efficiency,  as  a  func- 
i  i  mi  o’'  frequency .  This  i.  est  done  bv  calculating  the 

vuiues  of  the  various  parameters  indicated  in  the  figure 
,i,d  i  i  ut.it if.  a  graph  like  that  in  ’-’iq.  LI,  based  on  those 
vs  I ues  . 


Find  the  power  input  W^n  from  the  appropriate  form  of  Eq.  (5). 

Calculate  the  radiated  power  W  (at  each  frequency)  from  Eq . 

(8)  (using  the  structural  loss  factor  value  n  *  0.01  unless 
better  information  is  available)  and  determine  the  sound 
power  level  from  L  »  10  log  (W/Wref). 
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APPENDIX  M 


ATMOSPHERIC  PROPAGATION  EFFECTS 


INTRODUCTION 

The  aural  detectability  of  an  aircraft  Is  determined  not 
only  by  the  acoustic  characteristics  of  the  source  and  by  the 
acoustic  environment  and  sensitivity  of  the  receiver,  but  also 
by  the  propagation  characteristics  of  the  intervening  atmosphere. 
This  section  summarizes  the  propagation  effects  which  may  affect 
the  detectability  of  a  small  "quiet"  aircraft  approaching  a 
ground  observer  at  an  altitude  of  a  few  hundred  feet.  The  slant 
ranges  of  interest  here  are  of  the  order  of  a  few  thousand  feet, 
and  the  angles  of  elevation  6e  typically  are  small  (usually  less 
than  20'1). 

Most  of  the  available  literature  on  sound  propagation  is 
concerned  with  either  overhead  flight  (6e  >  20°)  of  aircraft  or 
with  ground-to-ground  propagation  (6  =  0°).  For  the  case  of 

overhead  flight,  only  the  effects  of  spreading  and  atmospheric 
absorption  are  significant.  For  ground-to-ground  propagation, 
the  effects  of  ground  absorption  and  atmospheric  refraction  must 
also  be  considered  and  may,  in  fact,  make  the  principal  contri¬ 
bution  to  attentuatlon . 

one  may  consider  two  types  of  propagation  factors:  (1) 

"prec  I y table"  factors,  which  occur  at  all  times  and  at  all  loca¬ 
tions,  and  (2)  "variable"  factors,  which  vary  with  time  and  lo¬ 
cation.  The  predictable  factors  include  spreading  and  atmospher¬ 
ic  absorption;  spreading  is  independent  of  all  meteorol ogical 
variables,  and  absorption  can  be  predicted  from  given  temperature 
and  humidity  profiles.  The  variable  factors  include  terrain 
attenuation,  attenuation  due  to  turbulent  scattering,  attenuation 
btv  fog  anu  rain,  focusing  and  formation  of  shadow  zones  by  re¬ 
fract  ion;  all  of  these  depend  on  the  local  terrain  and/or  pre¬ 
vailing  meteorological  conditions. 

in  addition  to  attenuating  an  acoustic  signal,  propagation 
effect.:  may  alter  the  acoustic  "signature"  (i.e.  the  waveform, 
spectrum,  or  time-variation  of  the  spectrum)  of  a  source.  Al¬ 
though  this  consideration  is  of  greater  Importance  for  helicop¬ 
ter  audibility  than  for  quiet  aircraft,  it  will  also  be  con¬ 
sidered  here  briefly. 


"PREDICTABLE"  EFFECTS 
Spreading  Loss 

The  energy  from  an  acoustic  source  spreads  over  an  increas¬ 
ing  area  as  it  propagates  away  from  the  source.  This  spreading 
results  in  a  decrease  in  the  acoustic  intensity  with  increasing 
distance  from  the  source,  or  in  "spreading  loss"  AL  which  obeys 

ALg(dB)  =  20  log  (R/R#  )  (1) 

where  R  represents  the  range  from  the  source  to  the  observation 
point  and  R0  represents  the  range  to  a  reference  point  (usually 
at  unit  distance  from  the  source).  ALS  thus  indicates  how  much 
reduction  in  the  acoustic  pressure  one  measures  as  the  result 
of  moving  the  observation  point  from  R0  to  R. 

Atmospheric  Absorption 

Atmospheric  absorption  is  usually  described  in  terms  of  an 
absorption  coefficient  a ,  which  generally  is  measured  in  units 
like  dB/1000  ft  and  accounts  for  "excess  attenuation"  beyond 
that  due  to  spreading.  The  classical  absorption  coefficient 
acl  accounts  for  losses  due  to  viscosity,  heat  conduction,  heat 
radiation  and  diffusion;  it  is  negligible  at  frequencies  below 
several  kHz.  The  molecular  absorption  coefficient  amol  accounts 
for  absorption  of  energy  from  the  sound  field  by  internal  modes 
of  the  gas  molecules,  and  is  negligible  at  frequencies  up  to 
several  hundred  Hz. 

The  vast  literature  on  this  subject  has  been  discussed  in 
a  recent  review  (Ref.  1)  and  empirically  derived  corrections 
have  been  added  to  the  theoretical  results,  resulting  in  formu¬ 
lae  which  match  the  available  data  almost  perfectly.  The  re¬ 
sults  are  outlined  below  and  are  illustrated  in  Fig.  Ml, 

The  classical  absorption  coefficient  at  a  temperature  of 
59°F  obeys 

a  , (dB/1000  ft)  =  5.3  *  10“2 (f/1000)2  (2) 

c  1 

where  f  represents  the  frequency  in  Hz.  This  absorption  coeffi¬ 
cient  is  independent  of  humidity  and  increases  by  about  1%  for 
each  20CF  increase  in  temperature. 


The  molecular  absorption  coefficient  is  given  by 
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reference  atmospheric  pressure  =  14.7  psl 
absolute  temperature  (°R) 
reference  temperature  -  519nR 
frequency  (Hz) 


1 1  s  relation  is  nlotted  in  Fig.  Ml  for  the  particular  temperature:: 


;>e°  >  and  100°F.  Kxamlnation  of  the  curves  reveals  that  for  a 

tropical  climate  (temperature  approaching  100°F,  ,  humidity  greater 
than  40C()  the  coefficient  may  be  approximated  rather  precisely  by 


a  ,(dB/l'>:>0  ft)  =  2(f/1000) 
mol 


(4) 
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f  amain  reivresent.'  the  frequency  in  Hz. 
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"VARIABLE"  EFFECTS 
Focusing 

The  possibility  of  anomalously  high  intensities  over  dis¬ 
tant  regions  has  been  a  subject  of  great  interest  in  the  study 
of  ground-to-ground  propagation  of  rocket  launch  noise.  When 
wind  and  temperature  profiles  are  favorable,  propagation  gains 
of  around  15  dB  at  all  frequencies  may  be  observe.,  relative  to 
normally  observed  levels.  Such  gains  occur  when  the  effective 
sound  velocity  decreases  with  height  up  to  a  certain  altitude 
and  then  increases  with  height  above  that. 

In  general,  the  probability  that  focusing  will  decrease  the 
aor-ai  detection  ranp;e  of  small  low-flying  aircraft  is  small,  for 
tut-  following  reasons: 

1,  Even  for  ground-to-ground  propagation,  the  conditions  for 
focusing  occur  rarely.  (For  example,  in  static  firings 
of  a  rocket  at  NASA's  Marshall  Space  Flight  Center,  fo¬ 
cusing  occurred  only  in  about  l'i%  of  the  cases.  Ref.  2) 

As  the  source  is  raised,  the  probability  of  focus  forma¬ 
tion  is  decreased. 

2.  For  ground-to-ground  propagation,  the  distance  to  a  focal 
region  Is  rarely  less  than  two  or  three  miles  (Ref.  2). 

Any  ai r-to-ground  focal  regions  would  occur  at  even  greater 
distances.  If  an  aircraft  Is  detected  under  non-focusing 
conditions  at  a  distance  of  the  order  of  a  mile  or  less, 
then  some  signal  enhancement  at  a  distance  of  several  miles 
will  be  of  little  concern,  because  it  will  largely  be  can¬ 
celled  by  the  increased  spreading  and  absorption  losses. 


If  low-altitude  flights  are  envisioned  in  areas  where  the 
meteorological  conditions  may  present  an  appreciable  probability 
of  focusing,  then  it  may  be  worthwhile  to  obtain  meteorological 
data  in  order  to  minimize  this  probability.  Often  appropriate 
choices  of  aircraft  approach  azimuth  nay  be  available  to  minii  ' ze 
the  probability  of  focusing  in  the  target  area.  (See  Ref.  2.) 
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Terrain  Attenuation 


For  low  angles  of  elevation,  additional  attenuation  may  be 
provided  by  ground  cover.  Reference  3,  which  reports  measure¬ 
ments  of  ground-to-ground  propagation  through  a  jungle  in  Panama 
is  the  classical  paper  on  this  subject.  The  principal  results 
are  shown  in  Fig.  M2. 

The  important  question  for  aircraft  detectability  considera¬ 
tions  is  how  these  results  are  affected  as  the  source  elevation 
is  increased.  A  number  of  studies  Indicate  that  terrain  attenua¬ 
tion  decreases  monotonically  with  source  height  and  becomes  neg¬ 
ligible  at  an  elevation  angle  of  around  6°.  For  example,  Parkin 
and  Scholes  (Ref.  4)  found  that  the  attenuation  of  a  helicopter 
at  a  range  of  3000  ft  began  to  Increase  rapidly  as  it  descended 
below  300  ft  height.  Hubbard  and  Maglieri  (Ref.  5)  found  that 
the  detectability  of  a  single-propeller  aircraft  decreased  sub¬ 
stantially  with  decreasing  elevation  angle  and  attributed  this 
to  the  increased  terrain  attenuation;  they  also  conducted  a 
number  of  measurements  to  determine  the  functional  dependence 
of  the  terrain  loss  coefficient  on  elevation  angle  (Fig.  M3). 
Loewy  (Ref.  6)  has  attempted  to  develop  a  general  procedure  to 
account  for  the  effect  of  elevation  angle,  but  his  procedure 
appears  to  be  based  on  insufficient  Information  and  has  not  been 
validated. 

Although  a  great  deal  of  work  remains  to  be  done  before 
reliable  quantitative  predictions  can  be  made,  it  is  clear 
qualitatively  that  the  aural  detectability  of  an  aircraft  can 
be  reduced  by  flying  along  a  path  along  which  the  ground  atten¬ 
uation  is  high  (e.g.,  where  there  is  heavy  foliage),  and/or 
by  flying  at  as  low  an  altitude  as  possible.  If  a  choice  is 
available,  detectability  can  be  decreased  by  approaching  the 
receiver  or.  an  azimuth  along  which  the  foliage  is  heaviest.  For 
low-altitude  flight  over  a  forest  or  Jungle,  terrain  attenuation 
may  far  outweigh  the  combined  effects  of  spreading  and  atmospher¬ 
ic  absorption. 

Attenuation  Due  to  Turbulent  Scattering 

At  low  freouencies,  additional  attenuation  may  be  caused 
by  scattering  from  turbulent  fluctuations  in  the  wind  and  tem¬ 
perature  field.  For  Isotropic  turbulence,  the  scattering  at¬ 
tenuation  coefficient  is  proportional  to  frequency  and  to  the 
gustiness  (i.e.,  the  root-mean-square  wind  speed  fluctuation, 
divided  by  the  mean  wind  speed). 


For  an  rms  wind  fluctuation  of  one  m/sec,  Horiuchi  (Ref.  7) 
has  estimated  the  attenuation  coefficient  to  be  on  the  order  of 
.002  f  dB/1000  ft.  Based  on  experiments  conducted  at  the  Marshal 
Space  Flight  Center,  the  following  suggested  minimum  values  were 
obtained  for  the  average  meteorological  conditions  prevailing 
there  (Ref.  1): 

a  . (dB/1000  ft)  =  0.03  for  f  <  7  H2 

o  Cd  v 

=  0.0042f  for  7  <  f  <  60  Hz 
=  0.25  for  60  <  f  <  200  Hz 


At  these  low  frequencies,  the  attenuation  due  to  turbulent 
scattering  may  be  greater  than  that  due  to  absorption;  at  higher 
frequencies,  absorption  effects  tend  to  predominate. 

Shadow  Zones 

Refraction  by  wind  and  temperature  gradients  can  sometimes 
lead  to  "shadow  zone"  regions  in  which  the  sound  levels  are 
greatly  reduced.  The  distance  Xg  to  the  beginning  of  the  shadow 
zone  is  given  approximately  by 


where  hg  is  the  mean  height  of  source  and  receiver,  c  is  the 
sound  velocity  and  du/dz  is  the  gradient  of  effective  sound  ve¬ 
locity  (Ref.  8).  For  purposes  of  minimizing  detectability,  one 
wishes  to  decrease  X  as  much  as  possible.  This  can  he  accom¬ 
plished  by  flying  atslow  altitude  into  the  wind.  However,  even 
for  a  reasonably  large  gradient  of  0.01  ft/sec  Tier  ft,  one  ob¬ 
tains  a  shadow  zone  at  a  distance  as  large  as  about  two  miles, 
for  a  source  height  of  ^00  ft. 

An  ideal  shadow  zone  occurs  only  in  the  limiting  case  of 
geometric  acoustics  (hirh  rre.iuer:c  1  es ) .  In  practice,  appre¬ 
ciable  energy  is  diffracted  and  scattered  into  a  shadow  zone. 
Wave  acoustics  predicts  an  attenuation  coefficient  inside  t no 
shadow  which  is  proportional  to  f1/3  (Ref.  9)  so  that  high  fre¬ 
quencies  are  attenuated  most. 


Attenuation  Due  to  Fog  and  Rain 

The  presence  of  suspended  water  droplets  may  cause  addi¬ 
tional  attenuation  through  the  mechanisms  of  viscous  dissipation 
and  heat  conduction.  The  theoretical  and  experimental  facts 
have  been  reviewed  in  Ref.  10.  Attenuations  on  the  order  of 
O.h  dP/1000  ft  (increasing  with  frequency)  may  be  obtained  for 
propagation  through  heavy  natural  fogs.  Analysis  of  measurements 
reported  in  Ref.  10  indicates  that  the  effect  of  light  rain  on 
sound  propagation  is  very  small. 

EFFECTS  OF  PROPAGATION  ON  PRESSURE  SIGNATURES 

Aircraft  detectability  may  be  increased  by  the  existence 
of  easily  recognizable  regularities  in  the  acoustic  (pressure 
vs.  time)  "signature"  of  the  source.  These  regularities  are  of 
primary  importance  for  helicopter  detection  (Ref.  6),  but  may 
also  play  a  significant  role  for  other  types  of  aircraft.  There- 
fore/a  brief  description  of  the  effect  of  the  propagation  chan¬ 
nel  on  acoustic  signatures  is  also  included  here. 

Random  inhomogeneities  in  the  refractive  index  cause 
frequency-dependent  fluctuations  in  the  amplitude  and  phase 
acoustic  waver.  Expressions  for  these  fluctuations  have 
been  derived  by  Tatarski  (Ref.  11).  The  effect  on  the  signal 
<"orm  can  be  determined  by  introducing  these  expressions  into 
the  spectral  representation  of  the  signal  and  then  calculating 
a  statistical  average.  This  procedure  has  been  carried  out  by 
.'hirkova  (Ref.  12).  For  simple  wave  forms,  such  as  a  rectangular 
ruise,  he  obtained  analytical  formulas  for  the  distortion.  The 
distortion  of  an  arbitrary  waveform  can  be  calculated  by  numeri¬ 
cal  methods. 

A  major  criticism  of  Shirokova's  work  is  that  tie  assumes  a 
normal  distribution  for  the  correlation  coefficient  of  the  re¬ 
fractive  index  fluctuations.  however,  Tatarski  has  shown  that 
a  correlation  proportional  to  the  2/3  power  of  distance  provides 
■i  pore  accurate  description  or  the  fluctuations.  'therefore, 
ohirokova's  formulae  should  be  revised  to  account  for  this. 

Tn  general,  one  may  expect  larger  distortion  of  signal 
shapes  when  the  conditions  near  the  ground  are  unstable.  Insta¬ 
bility  is  usually  characterized  by  gusty  winds  and  superadiabat io 
lanse' rates ,  which  result  in  turbulence  near  the  ground.  Experi¬ 
mental  evidence  of  this  effect  lias  been  obtained  by  observing 
ft,,,  distortion  of  pressure  wave  forms  from  sonic  booms  (Ref.  13). 
Karlv  in  the  morning  when  conditions  near  the  ground  are  general¬ 
ly  stable,  the  wave  forms  closely  resembled  the  theoretical  "N" 
■-.naoe.  In  the  afternoon  when  superadiabati c  conditions  prevail, 
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wide  variations  ranging  from  spiked  to  rounded  signatures  were 
noted.  In  fact,  the  magnitude  of  the  correlation  coefficient 
to  be  used  in  the  formula  discussed  above  can  be  expressed  di¬ 
rectly  in  terms  of  the  gradients  of  wind  and  temperature  near 
the  ground  (Ref.  11). 

ATTENUATION  CALCULATIONS 

Sound  Pressure  at  a  Given  Range 

The  reference  quantities  for  sound  power  level  Lw  and  sound 
pressure  level  LD  are  such  that  the  sound  pressure  level 
Lp ( Ro  =  1  ft)  at  1  ft  from  a  point  source  is  related  to  the  power 
level  of  that  source  as  (e.g..  Ref.  1*0 

Lp(R0)  =  Lw  -  0.5  dB  (6) 

The  sound  pressure  level  at  any  range  R  from  the  source  may  then 
be  calculated  from 


VR)  -  yv  -  y  m 

where  the  total  attenuation  AL^  is  given  by 

AL  =  AL  +  a..  ♦  R.  (8) 

p  s  tot 

Here,  AL  represents  the  loss  due  to  spreading  and  atot  the  ab“ 
sorptionRcoef ficient  (in  dB/unit  distance)  resulting  from  all 
other  effects. 

By  combining  Eqs.  (1)  and  (6)-(8)  and  expressing  R  in  feet, 
one  obtains  the  sound  pressure  level  Lp(R)  at  a  range  of  R  ft 
from  a  source  with  acoustic  power  level  Lw  as 

Lp(R)  =  Lw  -  20  log  R  -  c*tot(c!B/ft)  •  R  -  0.5  (9) 

Range  for  Given  Sound  Pressure 

By  substitution  into  the  foregoing  equation  one  may  readily 
determine  the  sound  pressure  that  results  at  a  given  distance 
from  a  given  source.  However,  determining  the  range  at  which  a 
prescribed  sound  pressure  is  obtained  from  a  given  source  is 


more  difficult,  because  Eq.  (9)  cannot  be  solved  algebraically 
for  R.  Figure  M4 has  been  prepared  to  facilitate  the  numerical 
solution  of  this  inverse  problem;  it  permits  one  to  obtain  a 
"corrected"  value#R  which  corresponds  to  a  given  value  of  a^ot 
from  the  easily  calculated  "uncorrected"  value  Ru  that  one 
obtains  from  Eq.  (9)  for  *  0. 


•The  corrected  and  uncorrected  ranges  are  related  by 
20  log  Ru  =  20  log  Ec  +  atotRc 

and,  of'  course, 

20  log  Ru  =  Lw  -  Lp  -  0.5 

where  both  Lw  and  Lp  are  riven. 
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FIGURE  Mia 


VALUES  OF  AIR  ABSORPTION  COEFFICIENTS  AT  59°F  AND 
VARIOUS  RELATIVE  HUMIDITIES  (FROM  REFERENCE  1) 
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GURE  M2a  TERRAIN  LOSS  COEFFICIENTS  FOR  GRASS  AREA.  (A)  OVER  TERRAIN 
THINLY  COVERED  WITH  SHORT  (6-12IN.)  GRASS;  (B)  OVER  TERRAIN 
THICKLY  COVERED  WITH  SHORT  ( 1 81 N .)  GRASS ;  (C)  THROUGH  A 
THICK  STAND  OF  SHORT  GRASS;  (D)  THROUGH  SHRUBBERY  AND  OVER 
SHORT  GRASS;  (E)  THROUGH  THICK  TALL  { 6 FT .  )  GRASS  (FROM 
REFERENCE  3) 
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FIGURE  M2c  A  CHART  FROM  WHICH  TO  ESTIMATE  TERRAIN  LOSS  COEFFICIENTS 
FOR  TROPICAL  JUNGLES.  ZONE  1,  VERY  LEAFY,  ONE  SEES  A 
DISTANCE  OF  APPROXIMATELY  20  FT.,  PENETRATION  BY  CUTTING; 
ZONE  2,  VERY  LEAFY,  ONE  SEES  APPROXIMATELY  50  FT., 
PENETRATED  WITH  DIFFICULTY  BUT  WITHOUT  CUTTING;  ZONE  3, 
LEAFY,  ONE  SFES  A  DISTANCE  OF  APPROXIMATELY  10(1  FT.,  FREE 
WALKING  IF  CARE  IS  TAKEN;  ZONE  4,  LEAFY,  ONE  SEES  A  DIS¬ 
TANCE  OF  APPROXIMATELY  200  FT.,  PENETRATION  IS  RATHER 
EASY;  ZONE  5,  LITTLE  LEAFY  UNDERGROWTH,  LARGE  BRACKETED 
TRUNKS,  ONE  SEES  A  DISTANCE  OF  APPROXIMATELY  300  FT., 
PENETRATION  IS  EASY  (FROM  REFERENCE  3). 
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AURAL  DETECTABILITY 


INTRODUCTION 

\ 

In  the  absence  of  any  ambient  noise,  an  aircraft  may  be 
heard  if  the  level  of  its  noise  at  the  location  of  the  listener 
exceeds  the  threshold  of  hearing  of  the  listener  in  any  frequency 
range.  The  presence  of  the  ambient  noise  tends  to  raise  the 
threshold  of  detectability  to  a  higher  level,  which  may  be  called 
the  "detection  threshold".  Increased  intensity  of  the  ambient  noise 
in  a  given  frequency  band  causes  the  detection  threshold  to  exceed 
the  threshold  of  hearing  by  increasing  amounts. 

Generally  speaking,  an  aircraft  is  aurally  detectable  if  its 
acoustical  signal  at  the  location  of  the  listener  exceeds  both 
the  threshold  of  hearing  and  the  masking  threshold.  Consequently, 
a  necessary  condition  for  aural  detection  is  that  the  signal  ex¬ 
ceeds  the  threshold  of  hearing  of  the  observer  in  at  least  one 
frequency  band;  a  sufficient  condition  is  that  the  acoustical 
signal  exceeds  also  the  detection  threshold  corresponding  to  the 
ambient  noise  level  present  at  the  time  of  observation  at  the 
location  of  the  listener. 

Since  the  ambient  noise  varies  from  place  to  place  -  and 
also  varies  with  the  time  of  day  at  any  given  location  -  the 
aural  detectability  of  a  given  aircraft  differs  for  different 
situations . 

In  addition,  aural  detection  is  affected  by  such  psycholo¬ 
gical  factors  as  the  observer's  familiarity  with  the  aircraft's 
acoustic  signature  (i.e.,  whether  he  knows  what  he  is  supposed 
to  listen  for),  his  other  duties  (which  may  interfere  with  his 
listening  concentration),  and  the  penalties  associated  with 
failure  to  hear  an  approaching  aircraft  as  early  as  possible. 

Because  of  the  many  variables  involved  in  determination  of 
detectability  criteria  for  a  complex  signal  in  the  presence  of 
fluctuating  ambient  noise,  the  problem  is  a  very  complex  one. 

In  order  to  arrive  at  an  engineering  approximation  of  detecta¬ 
bility  described  in  this  report,  it  was  necessary  to  rely  on 
laboratory  data  of  the  threshold  of  hearing  for  pure  tones  and 
bands  of  noise,  and  on  the  concept  of  critical  bandwidth  to 
evaluate  the  masking  effect  of  the  ambient  noise.  The  approach 
suggested  here  is  likely  to  lead  to  significant  errors  in  pre¬ 
dicting  detectability  in  actual  field  situations,  but  is  expected 
to  be  adequate  for  the  purpose  of  comparing  the  detectabilities 
of  alternate  aircraft  designs. 


The  acoustical  signature  of  an  aircraft  generally  consists 
of  pure  tones  and  broadband  noise.  Since  either  type  of  compon¬ 
ent  can  be  responsjble  for  aural  detectability,  both  must  be 
considered.  They  are  discussed  separately  in  the  following 
paragraphs . 

DETECTABILITY  OF  PURE  TONES 

Hearing  Threshold 

As  has  been  mentioned,  the  level  at  which  a  signal  becomes 
detectable  depends  on  the  threshold  of  hearing  and  on  the  masking 
level  provided  by  the  ambient  noise.  The  threshold  of  hearing 
for  pure  tones  has  been  investigated  extensively  in  the  labora¬ 
tory.  Figure  N1  shows  the  threshold  of  hearing  for  pure  tones,  as 
a  function  of  frequency,  as  recommended  by  the  International 
Standard  Organization  (Ref.  1).  This  figure  pertains  to  an 
"average"  young  listener,  and  indicates  that  the  ear  is  most 
sensitive  to  frequencies  between  about  3000  Hz  and^  5000  Hz. 

Masking  By  Noise 

Hawkins  and  Stevens  (Ref.  2)  measured  the  masking  effect  of 
white  noise  on  a  pure  tone  signal  by  exposing  test  subjects  simul¬ 
taneously  to  the  pure-tone  signal  and  the  masking  noise  and  noting 
at  which  noise  levels  the  subjects  could  just  recognize  the  tone 
as  having  a  definite  pitch.  Figure  N2  shows  the  observed  differ¬ 
ence  between  the  pure  tone  level  and  the  masking  white  noise 
level  (in  1  Hz  bands). 

There  exists  some  evidence  that  the  human  auditory  system 
senses  separately  the  energies  in  certain  "critical  bands"  (whose 
bandwidths  are  such  that  the  energy  in  the  masking  noise  in  these 
bands  matches  the  energy  in  the  just  audible  pure  tone  at  the 
band's  center  frequency).  Accordingly,  one  may  treat  various 
pure-tone  components  separately  and  use  Fig.  N2  as  the  basis  for 
an  acceptable  detectability  prediction  scheme. 

One  may  calculate  the  detection  level  L^f),  which  is  the 
level  of  a  pure  tone  which  is  just  detectable  in  the  presence  of 
noise  with  a  spectrum  level  -^(f),  simply  from 

Ld(f>  -  LN>1(f)  +  ec(f)  (i) 

where  0  , ( f )  is  a  correction  factor  which  is  plotted  in  Fig.  N2 
and  which  gives  one  the  noise  in  critical  bands.  The  spectrum 
(1  Hz  bandwidth)  level  L^j  ]_  is  related  to  the  noise  level  Ln  ^f 
measured  in  a  frequency  band  Af  (with  center  frequency)  as  * 


-  10  log  Af 


LN,1  =  LN , Af 

~  LN  oct  "  10  log  f  +  1.5  for  octave  bands 

=  ln)i/3  oct  “  10  log  f  +  ^.5  f°r  1/3-octave  bands.  (2) 
Detectability 

A  pure-tone  signal  thus  may  be  considered  detectable  in  the 
presence  of  noise,  if  the  signal  exceeds  both  the  detection  level 
Ld(f)  given  by  Eq.  (1)  and  the  hearing  threshold  of  Fig.  N1 . 

DETECTABILITY  OF  BROADBAND  NOISE 

Hearing  Threshold 

The  threshold  of  audibility  of  octave  bands  of  white  noise 
(constant  spectral  level  within  the  pass  band)  was  determined  by 
Robinson  and  Whittle  (Ref.  3).  Their  result  is  shown  as  the 
dotted  curve  of  Fig.  Nl.  Comparing  the  threshold  curves  for  pure 
tones  with  that  for  octave  bands,  one  finds  that  the  two  curves 
are  almost  identical.  The  difference  between  the  two  thresholds 
is  less  than  the  errors  in  estimation  procedures,  and  also  is 
much  smaller  than  the  variation  Ln  the  threshold  of  hearing  for 
individual  listeners.*  (It  should  be  noted,  however,  that  the 
octave  band  threshold  curve  was  obtained  for  constant  spectral 
level  within  the  pass  band  of  the  octave  filter.  Accordingly, 
these  results  are  strictly  applicable  only  for  similar  conditions.) 

Masking  By  Noise 

Like  a  pure-tone  signal,  a  broad.-band  signal  also  may  be 
masked  by  noise.  Unfortunately,  there  appears  to  be  available 
no  experimental  data  directly  applicable  to  predicting  the  mask¬ 
ing  of  a  broad-band  signal  by  broad-band  noise.  The  addition  of 
a  broad-band  signal  in  a  given  band  to  a  similar  noise  in  that 
band  would  be  sensed  by  a  listener  merely  as  an  increase  in  loud¬ 
ness,  unless  the  signal  has  a  spectrum  which  differs  from  that  of 
the  noise  —  in  which  case  the  listener  would  eventually  notice  a 
change  in  the  frequency  content  of  what  he  hears. 

One  may  perhaps  obtain  the  best  estimate  of  masking  effects 
by  comparing  the  levels  of  the  signal  and  noise  in  critical  bands. 


*Zwicker  and  Heinz  (Ref.  )  give  threshold  variances  of  5*5  dB  at 
50  Hz,  3.8  dB  at  1  kHz,  and  10  dB  at  12  kHz. 
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In  other  words,  if  one  has  a  signal  spectrum  Lo  *f  (f)  given  in 

frequency  bands  Afs  with  center  frequencies  f,  and  a  noise  signal 
Ln s A f n ( f ^  in  bands  Afn,  one  may  compute  the  corresponding  sound 

pressure  levels  in  critical  bands  from 

Ls,e<f>  '  Ls,if  <«  -  10  1os  Sfs  +  eo<f> 

(3) 

LN,e(f)  "  LN,Afn<f)  -  10  lQS  *fn  +  9c(f) 

Detectability 

/i 

The  signal  Lg  Af  (f)  may  be  considered  as  audible  if  its 
octave  band  level  1  s 

-  LS,oct  ■  LS,M's<f>  +  10  l0?  <f/lfs) 

exceeds  the  threshold  of  hearing  for  octave  bands,  and  if  in  ad¬ 
dition  the  sound  pressure  level  Lg  9(f)  of  the  signal  in  any 
critical  band  exceeds  the  level  LN|e(f)  in  the  same  critical  band. 

Note  that  if  the  signal  and  noise  are  given  in  the  same 
bands,  so  that  Afg  =  Afn  =  Af,  then 

Ls,e(f^  “  LN,6(f^  =  LS,Af(f)  "  LN,Af(f)‘ 

Therefore,  if  signal  and  noise  are  given  in  the  same  bands,  the 
signal  may  be  considered  detectable  if  it  exceeds  in  any  band 
both  the  hearing  threshold  and  the  noise.  In  other  words,  here 
the  band  detection  level  is  equal  to  the  noise  level  in  the  same 
band; 


Complications 


Superposition  on  ambient  noise  of  a  broadband  signal  which 
is  not  of  completely  random  character,  but  which  exhibits  an  am¬ 
plitude  modulation,  results  not  only  in  an  increase  in  intensity 
but  also  in  a  change  of  the  character  of  the  composite  signal. 
This  change  in  character  may  enable  a  listener  to  detect  the 
signal  at  a  lower  level  than  for  a  completely  random  signal. 

Feldkeller  and  Zwicker  (Ref.  5)  have  measured  the  effect  of 
modulation  on  the  detection  of  amplitude-modulated  pure  tones  and 
bands  of  random  noise.  For  bands  of  random  noise  with  sound 
pressure  levels  above  20  dB  and  center  frequencies  above  1  kHz, 
one  can  notice  amplitude  modulations  with  modulation  indexes  ex¬ 
ceeding  0.06. 

Since  in  these  experiments  only  a  single  band  of  amplitude- 
modulated  noise  was  presented  to  the  listener,  the  results  are 
not  immediately  applicable  to  the  aural  detection  of  aircraft, 
where  ambient  noise  is  present  in  all  bands  and  may  compete  for 
the  attention  of  the  listener.  However,  the  extreme  sensitivity 
of  the  ear  to  amplitude-modulation  should  be  taken  into  account 
in  designing  aircraft  for  minimum  aural  detectability. 

BACKGROUND  NOISE  LEVELS 

Data  on  the  ambient  noise  levels  in  Jungles  and  forested 
areas  are  presented  in  the  papers  of  Eyring  (Ref.  6),  McLaughlin 
and  Hand  (Ref.  7)  and  Saby  and  Thorpe  (Ref.  8).  The  data  given  in 
these  papers  excludes  transitory  sounds,  such  as  wind,  rain  and 
animal  calls  and  may  thus  be  used  to  obtain  conservative  esti¬ 
mates  of  aural  detectability. 

Eyring  presents  detection  levels  (^ig.  N3)  at  various  times 
of  day  for  a  dense,  leafy  jungle  and  for  a  forested  area  in 
Panama.  The  frequency  range  covered  is  from  100  to  6000  Hz.  His 
detection  levels  were  calculated  essentially  according  to  Eas.  (1) 
and  (2),  and  thus  apply  to  the  detection  of  pure  tones.  In  view 
of  Eqs .  (1),  (2),  and  (6),  the  corresponding  detection  levels  for 
bands  of  width  f  are  higher  than  those  shown  in  Fig.  N3  by  an 
amount  Ld  Af  ”  Ld  =  10  lo&  ^f  ~  0c^f^'  [Values  of  the  correction 
(Lh  -’Lh)  from  pure-tone  to  octave-band  masking  levels  are 
shcwnCin  Fig.  N2.J 

McLaughlin  and  Hand  present  one-third  octave  band  background 
levels  at  various  locations  and  times  in  a  Thailand  jungle,  for  a 
frequency  range  from  100  to  1000  Hz.  Although  this  paper  is  not 
readily  available,  some  representative  data  for  average  day  and 
night  conditions  have  been  reproduced  in  Ref.  9,  Saby  and  Thorpe 
present  spectrum  levels  of  background  noise  in  a  Fanama  jungle, 
for  high  frequencies,  from  8,000  to  25,000  Hz. 
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In  Pig.  N4  are  shown  the  combined  results  derived  from  the 
aforementioned  sources,  all  in  terms  of  detection  levels  for  pure 
tones.  The  two  figures  apply  to  average  daytime  and  nighttime 
conditions*  respectively.  (For  purposes  of  comparison.  Pig.  N4 
also  includes  a  curve  for  the  minimum  nighttime  detection  levels 
for  a  residential  area,  as  derived  from  the  background  noise  data 
of  Ref.  10.) 

The  various  sets  of  data  are  quite  consistent  and  indicate 
that  (a)  detection  levels  increase  with  increasing  density  of 
vegatation,  at  all  frequencies;  (b)  nighttime  detection  levels  are 
lower  than  daytime  levels  at  low  frequencies,  but  are  much  higher 
at  high  frequencies.  The  latter  effect  is  attributable  to  the 
nearly  continuous  presence  of  high-pitched  inset  sounds  during 
the  night.  Over  most  of  the  frequency  range,  the  detection  levels 
are  highest  during  early  evening  hours  and  lowest  around  midday. 

The  variation  of  detection  level  with  time  of  year  was  not  investi¬ 
gated,  but  would  be  expected  to  be  small  relative  to  the  diurnal 
variations . 

An  interesting  question  concerns  the  effect  of  the  human 
listener  on  these  detection  levels,  since  the  reported  measure¬ 
ments  were  all  made  with  a  microphone  placed  a  large  distance  from 
human  observers,  whereas  the  aural  detection  problem  involves  the 
detection  level  in  the  immediate  vicinity  of  an  observer.  Ac¬ 
cording  to  Ref.  8,  however,  the  quieting  of  nearby  insects  in  the 
presence  of  an  observer  does  not  measurably  alter  the  detection 
levels . 
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REPRESENTATIVE  DETECTION  LEVELS  IN  JUNGLES  AND  FORESTS 


APPENDIX  P 

ALTERNATE  METHOD  FOR  DETERMINING  UNCORRECTED 
DETECTION  RANGE 


INTRODUCTION 

The  "Detection  Level  Spectrum"  discussion  in  the  main  body 
of  this  report,  as  well  as  the  corresponding  aural  detectability 
estimation  approach  presented  in  Appendix  N,  are  based  on  the  re¬ 
sults  of  psychoacoustic  studies  of  the  detectability  of  stationary 
(i.e.,  time-independent)  pure-tone  and  broadband  signals.  Since 
different  detectability  criteria  apply  to  pure-tone  and  to  broad¬ 
band  signals,  one  must  approximate  the  noise  from  any  given 
source  in  terms  of  stationary  pure-tone  and  broadband  components, 
in  order  to  apply  the  suggested  approach. 

However,  the  noise  that  reaches  a  listener  on  the  ground  due 
to  an  aircraft  flyover  is  nonstationary  and  typically  consists  of 
many  broadband  and  "pure-tone"  components  that  chance  both  in 
apparent  frequency  and  amplitude.  Thus,  one  faces  the  often  very 
difficult  problem  of  estimating  equivalent  stationary  pure-tone 
and  broadband  levels.  The  various  noise  prediction  schemes  pre¬ 
sented  in  the  previous  appendices  suggest  how  these  levels  may 
be  estimated  for  the  various  noise  sources  considered,  and  may  be 
adequate  for  many  purposes  -  particularly  in  the  light  of  the 
often  considerable  uncertainties  associated  with  the  noise  source 
predictions  however,  an  alternate  method,  which  does  not  re- 
ouire  the  user  to  differentiate  a  priori  between  nure -tone  and 
broadband  signal  components,  has  been  found  useful  for  determin¬ 
ing  detectability,  especially  from  recorded  data  (Ref.  1).  This 
alternate  method  is  described  in  the  present  appendix. 

RESULTS  OF  SAILPLANE  AURAL  DETECTION  STUDY 

Reference  1  reports  the  results  of  a  study  of  the  flvby 
noise  (as  perceived  on  the  ground)  produced  by  three  different 
sailplanes  at  various  altitudes  and  speeds.  This  study  served  to 
characterize  the  noise,  but  also  Included  some  subjective  deter¬ 
minations  of  when  the  sailplanes  could  just  be  heard. 

It  was  found  that  the  observed  detection  ranges  agreed  well 
with  predictions  obtained  by  comparing  the  spectrum  levels  of 
the  received  sound  to  an  aural  detection  spectrum  for  pure  tones, 
where  the  spectrum  levels  were  determined  (by  means  of  an  appro¬ 
priate  Courier  analyzer)  with  an  effective  averaging,  time  of 
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50  milliseconds,*  and  were  obtained  by  reducing  the  bandwidth  of 
the  signal  analysis  until  the  level  remained  constant.  Since  this 
analysis  procedure  resulted  in  signal  levels  that  were  9  to  12  dB 
greater  than  the  spectrum  j.  tvels  calculated  from  the  third-octave 
levels  using  the  usual  bandwidth  conversion  relation,**  which  is 
based  on  the  assumption  that  the  energy  is  uniformly  distributed 
in  each  band,  it  appears  that  for  audibility  estimation  purposes 
one  should  increase  the  "constant  energy"  spectrum  levels  L^,  ^ 

obtained  from  broadband  spectra  by  about  10  dB  in  order  to  account 
for  the  ability  of  a  listener  to  detect  rapid  changes  in  frequency 
and  amplitude. 

ESTIMATION  OF  UNCORRECTED  DETECTION  RANGE  FROM  RECORDED 
FLYBY  SOUND 

In  order  to  determine  the  uncorrected  detection  range  of  an 
aircraft  from  tape-recorded  flyover  noise  data,  one  may  proceed 
as  described  below.  (An  example,  taken  from  Ref.  1,  is  shown  in 
Fig.  PI.) 

1.  Determine  the  pure-tone  detection  level  spectrum  L^(f)  from 

the  nure-tone  hearing  threshold  and  the  background  noise 
spectrum,  as  indicated  in  Appendix  N  and  Fig.  2  of  the  body 
of  this  report. 

2.  Determine  the  maximum  sound  pressure  level  in  each  one-third 
octave  band  during  the  flyover,  using  an  averaging  time  no 
greater  than  0.3  sec.  Listen  to  the  signal  during  the  data 
reduction,  and/or  carry  out  preliminary  narrow-band  analyses, 
judge  whether  pure  tones  are  dominant. 


*This  averaging  time  is  within  the  20  to  250  millisecond  range 
of  integration  times  of  the  ear,  as  reported  in  Ref.  2. 

**i.e . , 

LS,1  -  Ls,if  -  10  106  af  (1> 

where  Lg  ^  represents  the  level  measured  in  a  band  of  width  Af 
(e.g.,  a’third-octave  band),  and  Lg  1  denotes  the  corresponding 
"constant-energy"  spectrum  level. 
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3.  Find  the  one-third  octave  band  in  which  the  difference  be¬ 
tween  the  detection  level  (Step  1)  and  the  measured  third- 
octave  band  signal  level  (Step  2)  is  greatest.  (In  the  ex¬ 
ample  of  Fig.  PI,  this  greatest  difference  occurs  in  the 
band  centered  at  315  Hz  and  amounts  to  30.5  dB.) 

4.  Conduct  a  narrow-band  analysis  of  the  signal  over  a  fre¬ 
quency  range  that  includes  the  band  selected  in  Step  3, 
using  an  averaging  time  for  the  analysis  that  does  not  ex¬ 
ceed  0.1  sec. 

Repeat  this  analysis  with  narrower  and  narrower  analysis 
bandwidths,  until  the  peak  level  remains  essentially  un¬ 
changed.  Note  this  peak  level  and  the  frequency  at  which 
it  occurs.  (In  the  example  of  Fig.  PI,  this  spectrum  level 
is  37  dB,  at  285  Hz.) 

5.  Determine  the  difference  between  this  peak  level  and  the 
detection  level  spectrum  L(j(f)  at  the  peak  frequency.  (In 
Fig.  PI,  this  difference  is  20.5  dB.) 

6.  If  the  received  signal  is  not  dominated  by  pure  tones,  and 
a  narrow  band  analysis  cannot  be  conducted. 

(a)  Determine  the  effective  signal  spectrum  levels 

from  the  one-third  octave  band  spectrum  levels 

oct  found  in  Step  2,  by  use  of 


LS,1  =  lS,1/3  oct  -  10  log  f  +  16.5  (2) 


where  f  represents  the  center  frequency  of  the  third- 
octave  bands  (in  Hz). 

(b)  Consider  the  levels  Lg  i  to  occur  at  the  center  fre¬ 
quencies  of  the  corresponding  third-octave  bands,  and 
determine  the  greatest  difference  between  these  signal 
levels  and  the  pure-tone  detection  level  .  (In  the 
example  of  Fig.  PI,  the  greatest  difference,  22  dB  Is 
obtained  at  315  Hz). 

7.  Find  the  uncorrected  detection  range  Ru  from 

20  log  (Ru/Rmin>  =  LS,1  -  ^d  (3) 


-r - -  1 
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where  Rmin  is  the  minimum  distance  between  the  aircraft  and 

the  microphone  during  the  flyby  (and  is  equal  to  the  alti¬ 
tude  In  the  case  where  the  aircraft  passes  directly  over  the 
microphone),  and  Lg  ^-Ld  is  the  greatest  difference  between 

the  signal  spectrum  and  the  detection  spectrum  levels,  as 
obtained  from  Steps  3  and  k.  (In  the  example  of  Pig.  PI, 
uncorrected  detection  ranges  of  13^0  and  1580  ft  were  ob¬ 
tained.  ) 

Note  that  this  calculation  of  Ru  does  not  take  account  of 

atmospheric,  and  terrain  attenuation  effects.  Corrections  for 
these  effects  may  be  made  by  means  of  the  data  given  in  Appendix 
M.  If  Rmlp  is  large,  one  must  also  consider  that  these  attenua¬ 
tion  effects  can  affect  the  recorded  signals. 

ESTIMATION  OF  UNCORRECTED  DETECTION  RANGE  FROM 
PREDICTED  NOISE  SOURCE  CHARACTERISTICS 

In  order  to  determine  the  uncorrected  detection  range  from 
the  predictions  of  aircraft  noise,  as  given  in  the  various  fore¬ 
going  appendices,  one  may  proceed  as  follows: 

1.  Determine  the  pure-tone  detection  level  spectrum  Ld(f)  from 

the  pure-tone  hearing  threshold  and  the  applicable  back¬ 
ground  noise  spectrum,  as  described  in  Appendix  N  and  Fig.  2 
of  the  body  of  this  report. 

2.  Determine  the  power  level  and  frequency  for  each  pure-tone 
noise  component  from  each  source. 

3.  Determine  the  octave-band  power  level  spectrum  for  each  noise 
source,  and  combine  the  levels  from  all  of  the  sources  as 
described  in  the  body  of  this  report. 

Find  the  corresponding  effective  signal  spectrum  level  Ls  ^ 
from  ’ 

LS.l  -  LS , oct  -  10  lo«  f  +  X1-5  <"> 

where  t  represents  the  combined  octave-band  level  (as 

S ,  oc  t  ^ 

found  in  Step  3)  at  the  center  frequency  f. 

5.  Compare  the  power  levels  Lw  of  all  pure  tones  (from  Step  2) 
and  the  combined  effective  signal  spectrum  level  Lg  ^  (from 
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Step  *0  with  the  pure-tone  detection  level  spectrum  (from 
Step  1).  Find  the  greatest  differences  Lw-Ld  and  Ls 
and  the  corresponding  frequencies . 

Calculate  the  uncorrected  detection  range  Ru  from 
20  log  Ru(ft)  =  (L-Ld)max 

where  (L-L^)  represents  the  greater  of  the  two  maximum 

differences,  L  -L  ,  and  Lc  , -L, . 

*  w  d  S ,  1  d 
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Ml.  SUPPLEMENT  ARY  NOTES 


ITS ABSTRACT 


I  t  2  •  SPONSORING  MILITARY  ACTIVITY 


Air  Force  Flight  Dynamics  Lab. 
Wright-Patterson  Air  Force  Base 
Ohio 


*The  concepts  which  underlie  the  detection  of  aircraft  by  the  human  ear 
are  described.  A  scheme  is  delineated  for  predicting  the  range  at  which 
a  given  aircraft  will  first  be  heard  by  an  Jfaverage'^listener ;  this 
scheme  is  also  applicable  to  comparing  the  aural  detectabilities  of  al¬ 
ternate  aircraft  configurations  and  to  identifying  those  components  of 
the  noise  of  a  given  configuration  which  bear  prime  responsibility  for 
the  aircraft's  detectability.  Means  are  presented  for  predicting  the 
noise  due  to  all  sources  likely  to  be  significant  for  light  aircraft, 
the  attenuation  of  acoustic  signals  propagating  from  an  aircraft  to  a 
listener  on  the  ground,  and  the  ability  of  a  listener  to  detect  an 
acoustic  signal  in  the  presence  of  background  nolse._  ^  4  .  >ij 


ne^Crom  a 


This  report  gathers  together  for  the  first  tim^pQrom  a  diffuse  literature  of  ✓ 
acoustics,  engineering  mechanics  and  aeronautical  engineering ^kfiowlcdge  and  data 
on  methods  of  predicting  the  source  intensities,  propagation  and  aural  detection 
of  sound  radiated  by  military  aircraft  of  the  utility  class.  This  report  will  be 
used  for  the  design  of  quiet  aircraft  and  discloses  -frh£>de tailed  methods  used>  by 
the  Air  Force  in  evaluating  such  aircraft  from  the  acoustical  standpoint,  w. 
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